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‘lo fly faster 
than time itself... 


Trisonic jets that would land you in Los Angeles 2 hours earlier than your take-off 


time in New York are feasible says a late Federal study. Douglas officials say 


they could be operational by the early 1970’s. 


A 2100 mph civilian jet transport that would fly 
13 miles high, cross the continent in one hour and 
thirty minutes*, and use present jet runways is on 
the drawing boards at Douglas. 

Such an airplane is needed —says a recent 
Federal Aviation Agency study made with White 
House approval —to maintain U.S. leadership in 
commercial aviation. This is important because 
the export value of aircraft and parts in 1960 was 
$1.4 billion or 5.2% of total U.S. exports! 

The study also notes that substantial govern- 
ment assistance would be needed to underwrite the 
$500 to $550 million estimated development costs. 

Douglas believes that the estimated market of 


200 to 300 Mach 8 aircraft would more than repay 
these development costs. 

They are backing this belief with continuing 
studies based on 15 years experience with missiles, 
supersonic and hypersonic aircraft...to bring the 
trisonic civilian jet transport to reality at the 
earliest possible date. 


DOUGLAS 


MISSILE AND SPACE SYSTEMS « MILITARY AIRCRAFT * DC-8 
JETLINERS * RESEARCH AND DEVELOPMENT PROJECTS « 
GROUND SUPPORT EQUIPMENT « AIRCOMB® * ASW DEVICES 
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ARMA SYSTEMS 


Vibration exciters, precision centrifuges— 
these are among Arma environmental test 
facilities used to subject weapons systems to 
performances far beyond actual operation 
—to assure reliability, accuracy and long 
service life. 


predictable by test... proved by performance 


In the receiving areas, in the materials test labora- 
tory, on the production lines, in assembly and at the 
final test station ... wherever you look in the Arma 
plant you’ll see quality control and inspection equip- 
ment in use. Here rests the responsibility for main- 
taining the reliability of Arma systems. 

And, Arma systems are famous for reliability: All- 
inertial guidance for the Atlas ICBM, a proved system 
of intercontinental range accuracy ...the defense 


system for the B-52... shipboard fire control and 
navigation systems...anti-submarine warfare 
systems. 

Arma’s chief contribution to our national defense 
effort is the ability to deliver reliable systems—on 
schedule and within the budget. All research and 
development, engineering and production are based 
on that philosophy. Arma, a division of American 
Bosch Arma Corporation, Garden City, New York. 


AMERICAN BOSCH ARMA CORPORATION 
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An unsolved problem is a nagging challenge to him. The word “‘impossible”’ is an impertinence. 
Are you cut from the same cloth as this man? Then come to Northrop where you can work in 
the fringe of the future on such projects as space guidance and astronertial navigation systems, 
aerospace deceleration and landing systems, magnetogasdynamics, in-space rescue, repair and 
refueling techniques, laminar flow control, automatic test equipment and world-wide communications 
systems. 
More than 70 such programs are now on the boards at Northrop, with many challenging problems 
still to be solved, and new areas of activity constantly opening 
up for creative research. Whatever your specialty, drop us a line re (, - 
now if you want to find out more about the Northrop challenge. NORTHROP CORPORATION, BEVERLY HILLS. CALIFORNIA 
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The antimissile missile test program has been under way at 


White Sands Missile Range for some time, and further testing will be conducted from Point Mugu, 


Calif., where White Sands limitations on range and altitude do not exist. 


Finally, at Kwajalein 


Atoll, Zeus missiles will be launched, controlled, and directed to intercept ICBM target missiles 


fired from Vandenberg AFB, Calif. 


pany. 


Prime contractor for Nike Zeus is the Western Electric Com- 


Bell Telephone Laboratories is in charge of system research and design activities under 


direction of the Army Rocket and Guided Missile Agency. 
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ali, isa matter of 


Whether the “saddle shape,"’ 
matter of space travél, however, i is the subject . 
propulsion system, currently being studied at 
become the power. source for space vehicles. 


its design incorporates a nuclear reactor only 
1850°K. This is transmitted to banks of therm 
electrical energy for the ion. beam motor which 

without any moving parts, minimizing the po. 


Lockheed’s investigation of propulsion covers a by of potential systems. They include: plasma, 
ionic, nuclear, unique concepts in chemical systems involving high-energy solid and liquid propelients, 
combined solid-liquid chemical systems. The fundamentals of magnetohydrodynamics, as they might 
eventually apply to propulsion systems, are also being examined. Just as thoroughly, Lockheed probes 
all missile and space disciplines in depth. The extensive facilities of the research and development 
laboratories — together with the opportunity of working with men who are acknowledged leaders in 

their fields —make association with Lockheed truly rewarding and satisfying. 


Lockheed Missiles & Space Company in Sunnyvale and Palo Alto, on the beautiful San Francisco 

Peninsula, is an exciting and challenging place to. work. For further information, write Research and 
Development Staff, Department M-24B, 962 West Et Camino Real, Sunnyvale, California, U-S. citizenship 
or existing Departmentiof industrial security required, An Equal, Opportunity Employer. 
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ny shape ot ing conjecture. The 
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pckhead Missiles & Space Company, ini ht well 
ae pose ‘We diameter, generating heat at a temperature of 
i enerators, converting the heat directly into 
$s Cesium vapor as a fuel. The entire system is designed 
ibility of failure. { 
LOCKAHEE D MIssites SPACE COMPANY 
DIVISION OF LOCKMEED AIRCRAFT CORPORATION 
h Air Force AGENA Satellite in the DISCOVERER and mioas 
VENT and such NASA projects as OGO, OAO, ECHO, and NIMBUS. 
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GARRETT-U.S. AIR FORCE ‘SPUR’ 


--Oone answer to America’s future power needs in space 


E liminates 


giant size 
and weight 
barriers to 


large power 


systems 


Long lead time is essential to the 
development of large nuclear 
space power systems. Present meth- 
ods of power generation would require 
an impractical heat rejection surface 
nearly the size of a football field for a 
power output of one megawatt— power 
which will be needed for critical space 
missions already in the planning stage. 

Garreit’s AiResearch Divisions have 


THE 


Aerospace Engineering 


The joint U.S. Air Force-AEC Project SPUR is a 
research and development program being conducted 
by an industry team of Garrett, Aerojet-General 
Nucleonics (reactor) and Westinghouse, Lima 
(generator) to provide a power source to produce 
300-1000K W electrical power in space for one year or 
more. Power is obtained by conversion of nuclear 
fission energy to mechanical (shaft) power with a 
potassium vapor turbine. Equipment includes: reac- 
tor, primary and secondary loop pumps, boiler, turbo- 
generator and condenser-radiator. 


now completed the initial SPUR design 
studies and proved the project’s feasi- 
bility to supply continuous accessory 
power and low thrust electrical propul- 
sion in space for long periods of time. 

Cutting projected 1 MW power sys- 
tems to 1/10th the size and 1/5th the 
weight of present power systems under 
development will be possible because 
of SPUR’s capability to operate at 
higher temperatures, thereby sharply 

reducing the required radiator area. 


CORPO: 


Garrett has been working with the 
Air Force and the Atomic Energy 
Commission on SPUR as the prime 
contractor for more than one year and 
has more than five years of experience 
in space nuclear power development. 
Also an industry leader in high speed 
rotating machinery, heat transfer 
equipment, metallurgy and accessory 
power systems, the company is devel: 
oping design solutions for SPUR in 
these critical component system areas. 


AiResearch Manufacturing Divisions 


Los Angeles 45, California * Phoenix, Arizona 
Systems and Components for; AIRCRAFT, MISSILE, SPACECRAFT, ELECTRONIC, NUCLEAR AND INDUSTRIAL APPLICATIONS 
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The Next (X) Years—? 


ie IN AUGUST, as an integral part of the celebra- 
tion of the first half century of Naval Aviation, the 
U.S. Navy and the IAS held a notable joint exercise 
in San Diego. This was not the first time that the 
Institute collaborated with our Senior Service, 
but the occasion was unique in many ways. It was 
marked by the successful flying demonstration by 
Test Pilot Don Germeraad of the IAS-built replica 
of the Navy’s first airplane—the Curtiss A-1 (see 
IAS News, Sept. 1961). It included the dedication 
of an historic airfield to the memory of a great Naval 
hero—Admiral Halsey. But, more important from 
the Institute point of view, the 1,000-odd attendees at 
the three-day technical sessions heard scientists, 
engineers, economists, and policy planners discuss 
the blueprint for the future of the Navy’s Air Arm. 


The program was originally announced as a look 
at the ‘“‘next 50 years of Naval Aviation.’’ Almost 
toa man, however, participants backed away from 
such long-range prognostication on the grounds that 
the present state of affairs is so cloudy, and require- 
ments and capabilities are changing so rapidly, 
that the effective range of the crystal gazer is far 
short of a half century. 
ashot at the next decade 


Most were willing to take 
not beyond. It becomes 
too dangerous to base technical decisions upon ex- 


trapolations of extrapolations. 


The IAS Council is now facing up to a similar 
problem. During the past several years a great deal 
of time and energy have been spent by your Officers 
and Council in providing the necessary flexibility 
in programs and in operations to meet the require- 
ments of the IAS membership in this period of ex- 
ploding technology. Although the results have been 
good, much of the work was, of necessity, done on an 
ad hoc basis. Events moved so rapidly that time to 
sit down to think connectedly about the long-range 
luture was hard to come by. Long-Range Planning 
Committees were often concerned more with next 
year's budget than with questions as to where the 
Institute might be, or what it ought to look like a 
decade hence. 


EDITORIAL 


Now that some of the immediate problems have 
been dealt with, however, the attention of the Coun- 
cil has turned to the more distant future. At the 
last meeting (in Los Angeles, June 12) it was sug- 
guested that the next meeting be strictly a ‘‘think”’ 
session—that Officers and Council sit down together 
for a couple of days in some relatively quiet spot and 
plot the course for IAS for the foreseeable future. 
This struck a responsive chord, and plans are now 
well advanced for a week-end meeting of the Council 
and its Advisory Groups early in November, prob- 
ably in the Palo Alto area. 


It was early agreed (for the same reasons which 
influenced many prognostications at San Diego) 
that a “10-year look”’ was about the practical limit. 
Accordingly, the planning committee headed by 
IAS Vice-President Harold Luskin has come up 
with the Conference theme. JVhe IAS in 1970 
A Plan for Growth and Development. A four-part 
program has been laid down covering: 


(1) Current services—their 1970 character. 
(2) New and expanded services in 1970. 

(3) Implementation 
(4) Implementation 


IAS organization in 1970. 
LAS finances in 1970. 


The end point of the exercise will be the formula- 
tion of a specific program to carry the work of the 
Institute forward through the next decade. We 
want to know where we are going, and what must 
be done to arrive at our selected goals. 


During the past year we have made a number of 
extensive surveys of member opinion regarding IAS 
programs and policies. These have proved ex- 
tremely valuable as groundwork for the up-coming 
exercise. If, however, any member wishes to give 
the Council the benefit of additional thinking as to 
future IAS operations—now is the time. Letters 
addressed to any Officer or Council Member on any 
subject of common interest will be brought to the 
Council’s attention at the November meeting. Let 
us hear from you.—SPJ 
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Control & Guidance 


Piloted Simulation Studies of Re-Entry Guidance 
and Control at Parabolic Velocities 


Rodney C. Wingrove € Robert E. Coate 


Ames Research Center, NASA 


Tue GENERAL PROBLEM of concern here is that of 
designing pilot-operated systems to guide and con- 
trol maneuverable vehicles through a planetary 
atmosphere to a desired touchdown point without 
exceeding specified deceleration limits. The guid- 
ance methods most discussed in the literature use 
either predetermined stored reference trajectories or 
an on-board prediction computer in which future de- 
celerations and end points are repeatedly predicted 
as fast as possible from the present flight-path condi- 
tions. 

While the former ‘‘fixed-trajectory’’ method can 
be quite accurate, it is inherently limited to the con- 
ditions or situations stored on board. A requirement 
for a more general approach able to cope with a 
wider range of emergency conditions, including 
aborts, led to the study that is reported in refer- 
ence 

It is realized that an actual design choice between 
the two methods may depend on many factors out- 
side the scope of a preliminary simulation study; 
these include compatibility with the mid-course 
guidance system, status of component development, 
and system weight and volume. However, it ap- 
peared desirable after completing the study’ to per- 
form a reasonably parallel study of a fixed-trajec- 
tory method for the identical vehicle and conditions. 
It was hoped thereby to permit at least a comparison 
of range and entry-corridor capabilities and of ability 
to avoid exceeding the pilot’s time tolerance to de- 
celeration. 

Accordingly, both of these re-entry control tech- 
niques have been simulated, by use of an electronic 
analog computer. The vehicle is assumed to have a 
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Mr. Wingrove graduated with a B.S. degree 
in aeronautical engineering from the Univer- 
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time he has been in Government service and 
at present is a research scientist with NASA, 
Ames Research Center. Mr. Wingrove is 
the author or coauthor of several NACA- 
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and aerospace medicine technical papers. 
His early work was in the flight testing of 
experimental aircraft; stability and control 
analysis of handling qualities led to piloted- 
analog simulation; and his later work has 
been with fixed base, moving base, and centrifuge simulations for 
pilot control during re-entry. 


Mr. Coate, research scientist at the NASA 
Ames Research Center, received his A.B. 
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from the San Jose State College in June 
1950. Prior to his current employment he 
served as a pilot in the United States Air 
Force. Mr. Coate has specialized in manned- 
flight simulation techniques and has been 
associated with research in the fields of air- 
craft stability and control, and also in re- 
entry energy management systems. 


maximum L/D of 0.5 and to be entering the earth's 
atmosphere at an altitude of 290,000 ft at a velocity 
of 35,000 fps. Previous studies, references 2 through 
4 for example, have considered end-point control 
from near-earth orbital conditions. Re-entries at 
the speed chosen here, which is slightly less than 
parabolic, require that the control system damp the 
phugoid oscillation so that the vehicle does not 
exceed a maximum g limit (arbitrarily set at 10g), 
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» Results of piloted simulation studies are presented and compared for 
two piloted re-entry guidance methods for controlling lifting vehicles 
through a planetary atmosphere to a desired touchdown point. 
lt was found that the usable range with the fixed-trajectory method 
is about one third that with the prediction-control method, while the 
usable entry corridor is about the same for the two guidance systems. 
With either control system, this low L/D vehicle can be piloted in such 


a way that the pilot's time-tolerance-to-acceleration boundaries are not 


exceeded. 
Symbols 

A = acceleration factor, g units 

Cp = drag coefficient, D/[(1/2)pV2S] 

D = drag force, lbs 

g = gravitational acceleration unit, ft/sec? 

h = altitude, ft 

h = rate of climb, fps 

. Cp for desired trim condition 
K = drag ratio, ———————_ —— 
Cp for existing trim condition 

L_ = lift force, lbs 

m = mass of vehicle, slugs 

r = distance from planet center, ft 

R, = downrange value along great circle route in space, 
miles 

R, = crossrange value normal to local great circle route 
in space, miles 

S = surface area, ft? 

= circumferential velocity component normal to 
radius vector, ft/sec 

i = ratio of velocity to orbital velocity, u/+/ gr 

V = total velocity, ft/sec 

Y = side force, lbs 

Z = dimensionless function of 7 determined by Chapman 


equation and appropriate boundary conditions, 
(pu/2) (CpS/m)v/r/B 

8 = atmospheric density decay parameter, ft~! 

Y flight-path angle relative to local horizontal direc- 
tion; positive for climbing flight 


p atmosphere density, slugs/cu ft 
roll angle 
¥ = lateral deflection angle 
Subscripts 
d = respect to destination 
f respect to fixed trajectory 
1 = initial condition 


max = trim for maximum L/D 

min = trim for zero L/D 

0 condition at roll angle 0° 

u circumferential component normal to radius vector 


and also does not skip out of the atmosphere (repre- 
sented by a minimum limit of 0.075g). These L/D 
and g limits result in a permissible entry corridor of 
—2.1° to —6.2° flight-path angle (as determined by 
the methods of references 5, 6, and 7). 

In addition to the general comparison mentioned 
previously, this paper will show how the deceleration 
limits influence the range attainable by the vehicle 
and how varying the entry angle affects the control 
problems. The maximum heating the vehicle and 
pilot can tolerate is, of course, another important 
limit. Although heating was not considered in these 
studies, it can easily be taken into account in either 
control method; it affects the attainable range and 
permissible entry corridor in a manner analogous to 
the demonstrated effects of deceleration. 

(Continued on page 28) 
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Fig. 1. Maneuver capability with prediction-control method. 
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Flight Mechanics 


Supercircular Entry and Recovery With 
Maneuverable Manned Vehicles 


Robert H. Smith, MIAS € 


The Boeing Company, Aero-Space Division 


= PRIMARY OBJECTIVE of any manned space mis- 
sion is to perform some useful function in space at 
reasonable, if not minimum, cost. However, there is 
always an overriding responsibility to return the 
crew safely to their base. Resolution of these con- 
flicting demands without compromise of either is a 
formidable problem. Several approaches to this 
problem are under consideration; each has its par- 
ticular features and shortcomings. This paper, 
which is concerned with atmosphere entry and re- 
covery of manned space vehicles, emphasizes the 
factors of crew comfort and safety, direct return to 
home base, and minimum interference with space 
mission planning. 

Greater crew comfort and security can and should 
be provided by future entry vehicles. Current 
booster payload limitations, which place a premium 
on low vehicle weight, have almost forced 10g ac- 
celeration and ocean ditching into the status of stand- 
ard operating procedures for entry vehicles. Such 
operations will surely be relegated to their former 
emergency-procedure status when payloads can in- 
clude more sophisticated entry vehicles. 

Entry and recovery with a high lift/drag ratio 
glider will offer, along with other advantages, a high 
degree of crew comfort and security. A glider is 
tolerant of deviations from ideal entry conditions 
and allows control of discomfort from acceleration 
loads. After entry, transportation to the home base 
is provided by the long-range capability and ma- 
neuverability of the glider. These features allow ex- 
panded freedom in space mission scheduling. 
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Mr. Smith is engaged in aerodynamics re- 
search on the staff of the Flight Technology 
Department, Boeing Aero-Space Division. 
He works in the hypersonic aerodynamics and 
space system performance fields and partici- 
pated in the Dyna-Soar development during 
proposal stages. Before joining Boeing in 
1953 he was on the research staff of the 
University of Washington Aerodynamics 
Laboratories. Mr. Smith is a graduate of the 
University of Washington, with a B.S. degree 
in aeronautical engineering (1946) and an 
M.S. in aeronautical engineering (1952). 


Mr. Menard received a B.S. in aeronautical 
engineering (1954) and an M.S. in aero- 
nautical engineering (1960) from the Uni- 
versity of Washington. After serving 3 
years in the Air Force Research and Develop- 
ment Command as a pilot and aerial 
weapons test officer, he joined the NASA 
Flight Research Center as an Aeronautical 
Research Engineer where he participated in 
the investigation of the landing performance 
of low L/D vehicles. For the past 18 months 
he has been affiliated with the Flight Tech- 
nology Department, Boeing Aero-Space Divi- 
sion, working on superorbital entry vehicle performance. 


Flight Plan 


Direct single-pass entry will be used with con- 
trolled glide to the landing site. No power impulse 
is used and the glider stays in the atmosphere at all 
The mission is considered 
in three essentially independent phases: 

(1) Entry maneuvers involving a pull-out from 
the initial dive into the atmosphere and a pull-up 
to a suitable glide altitude just above satellite speed. 


times after initial entry. 
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A concept is presented for atmosphere entry at parabolic speed with a hypersonic 
glider. The concept provides for glide, after entry, within the atmosphere to a 
preselected landing site. Three principal advantages of using a high lift/drag ratio 
glider for manned return-from-space missions are established: 
(1) A very wide tolerance of entry angle is available if the crew endures accelera- 
tion loads approaching human limits. 
(2) Entry angle tolerance is adequate even if acceleration loads are held to levels 
which are relatively comfortable for the crew. 
(3) Controlled glide within the atmosphere to a preselected landing site is often 


possible regardless of the space-referenced location of the entry point. 


a (2) Global range control during trans-satellite lift away from the earth and tends to minimize 
on. glide. penetration of the atmosphere. The high drag as- 
= (3) Subsatellite banked glide to the recovery site sociated with maximum-lift attitude is desirable at 
ng with lateral offset control and minor range control this time to reduce speed. Both effects, speed re- 
_ maneuvers. duction and reduced penetration, are beneficial be- 
ics A logical fourth phase, not considered here, is the cause they reduce the peak dynamic pressure and 
- landing maneuver. Fig. 1 is an altitude-velocity heating rate which will be reached near the bottom 
an profile of typical trajectories to be considered. of the pull-out. When resultant acceleration load 


reaches a prescribed level, further increase is arrested 
by gradual reduction of angle of attack. During this 


Entry Maneuvers 


Pull-Out attitude modulation, the resultant of lift and drag 

The entry pull-out begins at 400,000 ft altitude is held constant; hence, resultant acceleration load 
and 35,000 fps with the glider in maximum-lift atti- is constant in spite of rapidly increasing dynamic 
tude. For steep (undershoot) entries, this directs pressure. This mode of (Continued on page 34) 
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Fig. 1 (left). Entry and recovery maneuvers. Fig. 2 (right). _Lift-drag polars. 
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Many honors, including his country's highest, have been bestowed 
upon this year's Gardner Lecturer. It was on May 19, 1942, that he 
received the Congressional Medal of Honor from President Roosevelt 
for having led the U.S. air raid on the Japanese mainland a month and 
a day earlier. Looking on (above left) are Mrs. Doolittle and Lt. 
Gen. H. H. Arnold. At right above is the 1961 Gardner Lecturer in 
his ‘open cockpit days." Gen. Doolittle recalls: “It was taken on 
Sept. 1, 1929, in Cleveland, just before taking off in a Curtiss P-1 to 
put on a stunting exhibition at the National Air Races. The airplane 
..-had a smaller radiator than the conventional P-1. As a result, it 
was somewhat cleaner and dived faster. While pushing under in a 
vertical dive the wings failed in reverse loading.’ Gen. Doolittle 
parachuted and about an hour later gave the scheduled exhibition in 
another airplane. in photo below, Gen. Doolittle receives the 


Daniel Guggenheim Medal the same year. Left to right are C. H. 
Colvin, G. W. Lewis, Brig. Gen. Doolittle, Glenn L. Martin, Joseph 
Hartson, Gen. Arnold, J. C. Hunsaker, and Lester D. Gardner. 


This year's Gardner Lecture was presented on April 28 at the Massa- 
chusetts Institute of Technology, and repeated on September 28 at the 
Smithsonian Institution. The assistance of both these organizations in 
scheduling this enlightening document for publication in AEROSPACE 
ENGINEERING is gratefully acknowledged. The Lecture bears the 
name—in the words of the 1961 Lec- 
turer—"of a great public servant and 
outstanding citizen,” Lester D. Gardner, 
whose long career in aviation included 
founding of the IAS in 1932. His death 
occurred in 1956 and a bequest to 
M.LT. provided funds to support a 
periodic address on the history of aero- 
nautics. The initial lecture was given in 1959 by Grover Loening; 
the second in the following year by Igor Sikorsky. Photographs, 
rather than a written biography, recite some of the highlights in the 
career of the 1961 Lecturer. J. H. “Jimmy” Doolittle, who lived 
through, and in fact made, some of the history about which he writes, 
now is Chairman of the Board of Space Technology Laboratories. 


On The Lecture, The 
Man Whose Name It 
Bears, And The 1961 
Lecturer 


TODAY 
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Early 


The Third 


J. H. Doolittle, HFIAS 


Space Technology Laboratories, Inc. 


- DEVELOPMENT OF 
flight has been a gradual, 
evolutionary process. 
There are, however, cer- 
tain landmarks along the 
route. We may arbi- 
trarily consider that 
there are roughly de- 
fined periods 
these datum points. 

In the very early days 
of flying, a slight breeze 
could cancel or delay a 
flight. 

I recall seeing my first air show in the winter of 
1909-10 at Dominguez Field near Los Angeles. 
In those old days it was customary for a pilot to 
wet his finger in his mouth and hold it up. I/ there 
was enough air movement to cause uneven evapora- 
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Blind Flying 


lester Gardner Lecture 


A Historical Review of 


Early Experiments in Instrument Flying 


tion and make one side cooler than the other, then 
there was too much wind to fly. 

Soon a modest wind was not a deterrent to flight 
and a good breeze was enjoyed because it shortened 
take-off and landing distance. 

Next came the period when a pilot was happy 
only as long as he could see the visible horizon. 
He might fly in or through clouds, but he wanted 
clear air and a visible horizon when he came back 
through them. 

This talk will deal with the next era, at the end of 
which a visible horizon was not required and during 
which it became possible to fly, and even land, with- 
out seeing outside an instrumented cockpit. 

Today aircraft fly safely and reliably in all but the 
most inclement weather, and I look forward to the 


Many personalities mentioned in Gen. Doolittle’s 1961 Gardner 
Lecture were present at both the initial presentation at M.I.T., and for 
its repeat at the Smithsonian Institution. These photos were taken at 
the M.L.T. event which took place on April 2. The first (opposite page, 


Instruments 


not-very-distant future when the airplane—or its 
successor—will fly absolutely regardless of weather 
and will then become not only the fastest but the 
most reliable form of transportation. 

Whenever pilots got together for meditation and 
discussion in the immediate post World War I era, 
the subject of “‘flying instinct’ or flying by the 
“seat of the pants’ was apt to come up. 

Views were divided. Some pilots believed that 
they could fly indefinitely without reference to the 
visible horizon—because they had. Others agreed 
that they had, but claimed it was not pilot skill but 
the inherent stability of their aircraft which made it 
possible. 

In the early twenties the navigational instruments 
most commonly used were (Continued on page 56) 


bottom of Col. 2) shows (from left to right) Jerome Lederer, Elmer 
Sperry, Jr., and C. Stark Draper. Center: Harry Guggenheim, 
James R. Killian, and Luis de Florez. Right: C. A. Frische, J. C. Hun- 
saker, and Vice Adm. Carl F. Espe. 
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Meteorology 


High-Altitude Wind Observations 


as RESPONSE of a space vehicle during its 
flight depends not only on wind velocity, but also on 
the rate of change of velocity with respect to time. 
This rate of change of velocity is essentially deter- 
mined by wind shear, which is the change of wind 
velocity with respect to altitude. Wind shear values 
have a significant influence on vehicle response 
even when associated with relatively small increases 
in wind velocity.! Optimum achievement of con- 
trol response depends upon a broad knowledge of 
wind velocity, wind shear, and the gusts to be en- 
countered. Current information concerning these 
factors is less adequate than is desired, especially 
for altitudes above 20 km where few observations 
have been made, and the accuracy of the data ob- 
tained doubtful. 
High-altitude wind 
Canaveral, Florida, are increasing in number, alti- 


measurements over Cape 
tude, and accuracy; however, wind data continue 
to be relatively scarce and expensive to obtain in 
comparison to routine rawinsonde measurements. 
During the period from April 1959 through May 
1960, 37 sounding rockets were launched successfully 
from Cape Canaveral. The objective of these firings 
was to gather information on high-altitude wind 
flow characteristics for use in missile and space 
vehicle design studies, and in the evaluation of test 
results obtained from ballistic missile test firings. 
The data* recorded are presented here, and, al- 
though the sample lacks sufficient size to have a de- 
finable statistical significance, it gives an indication 
of the general features of the atmospheric circulation 

* Data furnished through the cooperation of the U.S. Army 
Signal Research and Development Laboratories, the Air Force 


Cambridge Research Center, and the Air Force Missile Test 
Center. 
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George C. Marshall Space Flight Center, NASA 


at high altitudes over Cape Canaveral (Atlantic 
Missile Range). 


observations to 
several geographical areas has been conducted under 
the direction of the Meteorological Rocket Net- 
This effort, through the cooperation of the 
military services, the Weather Bureau, and the 
NASA, will provide valuable data in higher altitude 
regions. 


An extension of rocketsonde 


work.” 


Mr. Smith attended Oklahoma City University 
while working for the U.S. Weather Bureau, 
prior to World War Il. During the war he 
served 4 years with the Army Air Force as 
an Instructor in Meteorology and as Per- 
sonal Equipment Officer. He returned to the 
Weather Bureau after the war, and in 1957 
transferred to the U.S. Army Ballistic Missile 
Agency, Huntsville, Ala. He has been a 
Meteorologist in the Aeroballistics Division 
of the Marshall Space Flight Center, NASA, 
since the creation of that agency in 1960. 
Mr. Smith is a member of the American 
Meteorological Society. 


Mr. Vaughan received his B.S. degree in 
mathematics and physical sciences from the 
University of Florida in 1951, and did 
graduate work in meteorology at Florida 
State University while an officer with the 
USAF. In 1955, he joined the Air Force 
Armament Center, Eglin AFB, Fla., as an 
R&D Meteorologist. He then transferred to 
the ABMA, and since 1960 has been asso- 
ciated with the Aeroballistics Division, Mar- 
shall Space Flight Center, NASA. Active in 
the study and interpretation of atmospheric 
environmental data for development of mis- 


sile and space vehicle design criteria and test performance evaluc- 
tions, Mr. Vaughan is Chief, Aerophysics and Astrophysics Branch, Aero- 
He is a member of the American Meteorologi- 


ballistics Div., MSFC. 
cal Society and the American Geophysical Union. 
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Above Cape Canaveral 


A summary of 13 months’ wind data from rocketsonde observations above 20-km 


altitude over Atlantic Missile Range at Cape Canaveral. 
Wind speed, direction, shear, and components from April 12, 1959, to May 25, 1960, 


are analyzed on a seasonal, as well as annual, basis and the results are presented 


graphically. 
90 
= 80 
INOICATES WIND FROM NW AT IS METERS / SEC 
2 ™ WIND FROM W AT SO METERS / SEC 
> 
< 70 t = 
\ ay > u 
JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC JAN 
Fig. 1. High-altitude wind arrowgram by rocketsonde observation, April 12, 1959, to May 25, 1960, Cape Canaveral. 


Origin of Data and Methods of Computation 

The rockets used to obtain the data were of the 
unguided, solid-propellant type which eject a target 
hear apogee a few minutes after firing.* The targets 
which consisted of either chaff, parachutes, or bal- 
loons were tracked by radar. The resulting tracking 
data were resolved by electronic computer to provide 
wind speed, direction, shear, components, and rate-of- 
fall of the target. The radar provided a record of 


the time, azimuth angle, elevation angle, and slant 
range at one-sec intervals. These data were 
converted to rectangular coordinates, and the result- 
ing data were smoothed over 101 one-sec data 
points, using a second-degree polynomial. Allow- 
ance was made for the curvature of the earth. A 
three-point linear smoothing was then applied using 
30-sec time-spaced data points, and wind variables 
were determined at 60-sec (Continued on page 53) 
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Structures & Materials 


A Unified Engineering Theory of 
High Stress Level Fatigue 


Sitaram Rao Valluri, AFIAS 


California Institute of Technology 


siaus THE LAST hundred years of research in 
fatigue, several theories concerning fatigue failure 
were suggested which held sway for some time and 
were subsequently discarded in the light of more re- 
cent information found to be inconsistent with the 
available theories.* Needless to say, any theories 
proposed now are necessarily subject to the same 
fate eventually. However, during the short period 
of its life, a theory will have served its purpose if it 
can explain consistently most of the available experi- 
mental information, can logically predict some re- 
sults as yet unforeseen, and can be checked experi- 
mentally to establish its region of validity. It is 
with this end in view that the following theory of 
metal fatigue at high stress levels is proposed. 
Limitations of a theory are almost always im- 
posed by the explicit and implicit assumptions made 
in the development of the theory. In fatigue a test 
specimen or a structure may fail in the first loading 
or withstand several million cycles of loading before 
failing—or may not even fail at all. To be com- 


This paper has been based upon research supported by the 
Aeronautical Research Laboratories, Air Force Research 
Division, ARDC, USAF, at the Wright-Patterson AFB, Ohio. 
The theory proposed here is discussed more extensively in a 
series of reports which will be published by the Laboratories 
and which are at present available as GALCIT Solid Me- 
chanics Reports (references a, b, c, and d). 

It is a pleasure to acknowledge the sponsorship and financial 
assistance of this agency, and the many helpful discussions 
with Profs. Y. C. Fung, E. E. Sechler, and M. L. Williams of 
GALCIT in the preparation of this paper. Thanks are also 
due to R. H. Christensen of Douglas Aircraft Company, Inc., 
for kindly supplying Figs. 2, 3, and 4. 

* Symbols list appears on pages 68 and 69. 
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The author has been working in the Solid 
Mechanics group of C.I.T.'s Guggenheim 
Aeronautical Lab. since his graduation in 
1954; at present he is Senior Research 
Fellow there. He received the B.Sc. in 
mechanical and electrical engineering from 
Banares University in India, the D.I.I.Sc. from 
the Indian Institute of Science in Bangalore, 
and the M.S. and Ph.D. in aeronautics from 
C.1.T. Interested in metal fatigue and, more 
recently, in the problems of fast propagating 
cracks and the phenomenon of spalling, Dr. 
Valluri's primary interest has been to see if 
the knowledge compiled by physical metallurgists in the dislocation 
theory of metals can be used by engineers for a better understanding 
of the engineering behavior of metals. He is an Associate Fellow 
of the IAS and a member of the American Society of Metals, Society of 
Experimental Stress Analysis, and Sigma Xi. He is also a consultant 
to Douglas Aircraft Company. 


plete, a theory of fatigue, apart from explaining the 
dependence of fatigue upon the innumerable vari- 
ables, must encompass the wide range of stresses over 
which failure may take place. Existing basic ex- 
perimental observations show that failure under two 
fundamentally different deformational processes has 
been generally classified as fatigue. Under such 
circumstances, it becomes somewhat illusory to 
think in terms of a theory of fatigue that covers the 
whole range of stress and qualitatively, if not quan- 
titatively, explains the dependence of fatigue on the 
several variables. If this idea is accepted, the next 
problem is one of defining a suitable mechanism, 
determining its range of validity, and discussing the 
dependence of the several variables. If one theory 
can be set apart from others, it will be because of 
the simplicity of the chosen model, the range of its 
validity, and the completeness with which it de- 
scribes the dependence of the several variables. 

In fatigue the number of items a theory has to 
explain are staggeringly large if we insist that the 
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“An analytical theory of fatigue for high stress levels 
has been formulated based on a fairly simple model 
derived from a synthesis of the macroscopic elasto- 
plastic fracture theory and concepts derived from the 
dislocation theory of metals. The theory (1) predicts 
the shape of the stress vs. number of cycles curve in 
fatigue; (2) afler determining two of the unknown 
constants from test results, makes quantitative predic- 
lions regarding the residual static strength, changes re- 
sulting from changes of stress ratio (and mean stress), 
cumulative damage due either lo multiple step or to 
random loading; and (3) as a formal exlension, gives 
a quantitative theory of acoustic fatigue.” 


explanation of all the ideas must come from the 
same basic model with at most a few (very few!) 
assumptions. Some of the various aspects which 
should be explained are shape of the stress vs. 
number of cycles curve in fatigue; existence or non- 
existence of an endurance limit; and effects of cold 
working, grain size, frequency, temperature, over- 
stressing, understressing, and applying single over- 
loads. One should also be able to discuss residual 
static strength; effect of changing test stress ratio 
or mean stress; cumulative damage (specifically, 
response to stress magnitude and response to the 
order of application of stress); response to random 
loading; and problems of acoustic fatigue. In addi- 
tion, the theory should yield quantitative results 
with at most a few constants to be determined ex- 
perimentally. If not immediately amenable, it 
should at least have the elements required to inter- 
pret the statistical distribution so characteristic of 
fatigue. Finally, it must be simple conceptually if 
not analytically. 

It is nearly trite to comment that no theory has 
yet been devised which can do this over the whole 
stress spectrum in fatigue, especially for all the dif- 
ferent materials. It therefore becomes necessary to 
restrict the scope of the theory. In order to do this, 
we shall have recourse to a brief survey of the avail- 
able literature to ascertain basic facts that are suffi- 
ciently reliable, propose a hypothesis that is intui- 
tively plausible, construct a quantitative theory, and 
find out to what extent the result so obtained can 
discuss the various aspects of fatigue in a systematic 
and consistent manner. (Continued on page 68) 
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Fig. 1. Operation of a Frank-Read dislocation source leading 
crack formation. 


Time Sec 


Fig. 2. Crack growth vs. number of load cycles in a damaging stress 
range. 


Fig. 3. Schematic representation of fatigue crack growth. 
Note repetitive and alternate periods of micro and macro crack- 
ing which result in “jumps” or “spurts” during progressive frac- 
ture. (Numbers of cycles of loading at various stages are 
from an actual test.) 
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Fig. 1. Geometry of nose-wheel steering problem. 


Nearcy ALL currently operational aircraft utilize 
landing gear systems with nose-wheel steering capa- 
bilities. _ Nose-wheel steering minimizes ground 
handling (towing) requirements and facilitates park- 
ing, taxiing, alignment on runway, and ground roll 
during landing and take-off. Control of nose-wheel 
steering is handled in most aircraft by the pilot 
through a separate steering wheel, and in a few 
through a direct linkage with the rudder bar. In 
either instance the pilot retains complete directional 
control of the aircraft throughout its so-called 
ground phase. 

The essential thought of this paper is that a human 
operator (the pilot) must act as an element in a 
closed-loop servosystem. The functions of this 
human element are as follows: 

(1) Selection of a reference (point, line, angle, 
curve, etc.) from which to determine the aircraft's 
misalignment or deviation—i.e., error. 

(2) Detection of such errors. 

(3) Establishment of an error criterion. 

(4) Response to the detected error (via rudder 
bar or steering wheel) in such a manner as to main- 
tain the aircraft within the established criterion 
limits throughout the ground phase. 

These are complex functions for a servoelement 
of any type to perform. With this consideration in 
mind, it is a natural consequence to seek ways in 
which the ground steering task of the pilot could be 
made safer and simpler. 


Analysis 


In order to provide insight into the way in which 
the human functions as a servoelement in a nose- 
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Joseph G. Wohl 
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wheel steering system, the following simplified 
analysis of steering geometry is presented. 

Fig. 1 illustrates the geometry of the nose-wheel 
steering problem during take-off. Certain simplify- 
ing assumptions are used for convenience: 


6 < 10° 
sin 6 = 6 
y < 10° 
sin y = y 


v = a constant 


(Effects of suspensions, cross wind, wheels, tires, 
skid or slip, inertia, etc., are neglected.) 

Under these assumptions, the following simple and 
fundamental equations can be derived from the rela- 
tionships illustrated in Fig. 1: 


R = W/6é (1) 
y =v/R (2) 
t= of (3) 
where 
R = turn radius, ft 
W = wheel base, ft 
= 


nose-wheel angle, rad 

Y = direction of motion of aircraft with respect 
to runway centerline, rad 

v = ground speed, ft/sec 

& = dx/di 

= dy/dt 


Substituting Eq. (1) in Eq. (2) results in 
vy = v0/W 4) 


Taking the derivative of Eq. (3) with respect to 
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Dynamics of Nose-Wheel Steering 


Aircraft ground steering dynamics are analyzed from a human factors standpoint, 


and a velocity-modulated steering system is proposed to reduce required 


variations in human steering sensitivity during the ground phase. 


time (holding v constant) results in 
w (5) 
Substituting Eqs. (3) and (4) in Eq. (5) results in 
(6) 


Eq. (6) relates lateral acceleration with respect to 
the runway centerline to vehicle speed, wheel base, 
and nose-wheel angle, and represents a simplified 
model of the system dynamics at constant speed 
without the human operator. 

Although several possible avenues of approach 
may be taken in analyzing human operation in the 
steering loop, a good starting place is the perceived 
information upon which the pilot might base his con- 
tinuous judgment concerning steering. 

The human perceptual problem is indicated more 
clearly in Fig. 2, which shows that there are at 
least six perceptual variables (the three shown plus 
their rates of change) with which the man might 
have to deal on an instantaneous and continuing 
basis. 

From Fig. 2 we can infer that, for small angles 
< 10°), 


= Dy (7) 
Differentiating this once, we obtain 
x= Dy+ yD (8) 


Solving this expression for ¥ and substituting for % 
from Eq. (3), we get 


¥ = — D)/D] (9) 
Setting Eqs. (4) and (9) equal, we obtain 


Mr. Wohl is Associate Program Director 
whose responsibilities include the human 
factors program of the Fleet Ballistic Missile 
System. He received a B.S. degree in phys- 
ics and math from the University of Wis- 
consin in 1949 and has done graduate work 
at several universities. A Lecturer in man- 
machine systems engineering at Columbia 
University, he worked in the engineering 
psychology field at the Naval Research 
Laboratory and at ONR from 1948 to 1954, 
and for the ensuing 3 years was Senior 
Engineer and Group Leader at Sperry Gyro- 
scope Company, performing duties involving systems engineering and 
preliminary design of drone and missile systems. He joined Dunlap 
and Associates in 1957, and holds memberships in several professional 
societies. 


v0/W = — D)/D) (10) 
Solving Eq. (10) for 0, we obtain 
= ¥(W/2) — D)/D) (11) 


If as in many aircraft and ground vehicles the 
steering system is reversible,* then the nose-wheel 
displacement @ is proportional to the force F applied 
by a human operator to the rudder bar or steering 
wheel—viz., 


6= kF (12) 
Eq. (11) can then be rewritten as 
F = (y/k) (W/v) — D)/D] (13) 
(Continued on page 90) 


* In the zrreversible case—-e.g., certain aircraft flight control 
systems—the displacement @ may instead be proportional to 
the angle a through which the steering wheel or rudder bar 
is operated, with the disadvantage that there is no force 
feedback to the pilot. For a discussion of why the force- 
displacement concept is more desirable, see reference 3. 
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Aerodynamics & Fluid Mechanics 


Application of Laminar Flow Control 


to Transport Aircraft 


Welko E. Gasich, AFIAS, Northrop Corporation 


Mocu HAS BEEN WRITTEN about the aerodynamic, 
propulsive, and structural aspects of laminar flow 
control. This paper discusses some of the areas 
which have received relatively less emphasis, and 
which are now commanding greater attention among 
aircraft designers and airline operators. In order to 
determine the applicability of laminar flow control 
to transport aircraft, one must compare the operat- 
ing costs of normal aircraft having turbulent bound- 
ary layers with aircraft employing laminar flow 
control on wings and empennage. This considera- 
tion leads to the determination of manufacturing 
costs, maintenance costs, and operational considera- 
tions. 

In order to demonstrate the practicality of lami- 
nar flow control, the United States Air Force has 
awarded Northrop Corporation a contract to con- 
struct and fly two demonstration aircraft. A brief 
description of this airplane configuration and the 
program objectives will be discussed. 


LFC Definition 


The definition of laminar flow is illustrated in 
Fig. 1. When the boundary layer air, or thin layer 
of air adjacent to the wing surface, moves in prac- 
tically parallel air layers that do not intermix with 
each other, the air is said to be laminar and is of a 
low drag nature. 

Turbulent flow is characterized by oscillating or 
fluctuating air movements. These disturbances or 
energy interchanges within the turbulent boundary 


The author would like to extend his appreciation to S. H. 
Brown and M. G. Huben of the Advanced Design Staff and to 
Dr. W. Pfenninger, Chief, Boundary Layer Research, for 
their contributions to this paper. 
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layer are responsible for the relatively high tur- 
bulent friction drag. 

Laminar flow control (LFC) prevents the forma- 
tion of turbulent flow by pulling the slow-moving, 
turbulent-prone air adjacent to the wing surface 
through a large number of fine slots cut spanwise 
in the outer skin, and passing it through ducts inside 
the aircraft structure, then pumping it overboard. 

A schematic diagram of a typical LFC system is 
shown in Fig. 2. The wake drag of the laminar 
surfaces and the equivalent suction energy required 
for the LFC process is balanced by a portion of the 
main engine thrust, the remainder of the main en- 
gine thrust being used to overcome induced drag 
and the drag of the fuselage and engine nacelles. 

The equivalent suction energy C, that is required 
is accounted for in the skin friction drag coefficients 
C, quoted tor LFC surfaces, giving a true picture 
of the net benefit resulting from LFC. The suction 
compressor can be driven by a free turbine, direct 


Mr. Gasich received his A.B. (mechanical 
engineering) and his M.S. (aeronautics) at 
Stanford University, and his professional de- 
gree of Aeronautical Engineering at Cali- 
fornia Institute of Technology. He joined 
Northrop in 1953 as Chief, Preliminary De- 
sign, and was promoted to his present posi- 
tion of Director, Advanced Systems Section, 
in 1957. Preparatory to his work at 
Northrop were associations with the Ames 
Aeronautical Lab. (NACA) as flight test and 
research engineer, Douglas Aircraft as 
structures engineer specializing in cero- 
elastics, and RAND Corp. where he was responsible for the aircraft 
design group performing design synthesis of various system studies. 
In his current position, Mr. Gasich directs the systems analysis and ad- 
vanced design of new aerospace systems undertaken by the Norair 
Division. In addition to his affiliation with the IAS, he is a member of 
the Society of Automotive Engineers and Sigma Xi. 
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shaft, bleed and burn turbine or by a separate turbo- 
shaft engine.! 


Research Results 


As described in the preceding section, laminar flow 
control appears elementary in concept and, although 
its advantages were understood in the late 1930’s, 
it has commanded great attention only recently for 
the following reasons: 

(1) It was easier for many years to make dramatic 
reductions in drag with developments such as canti- 
lever wings, retractable landing gears, etc. These 
developments were a necessary part of the evolution 
to provide a smooth, streamlined airplane com- 
patible with LFC. 

(2) LFC required a comprehensive research pro- 
gram to bridge the gap between its theory and ap- 
plication. 

Yet, LFC appeared to be the last significant step 
remaining in the history of friction drag reduction 
and thus posed a challenge. To solve the problem, 
a joint Air Force-Northrop effort was initiated in 
1949 under the direction of Dr. Werner Pfenninger. 
The resulting 12-year program has covered extensive 
theoretical and experimental research, laboratory 
testing, wind-tunnel programs, and flight demon- 
strations to substantiate fully the friction drag 
figures quoted for LFC. 

Some program results concerning drag coefficient 
vs. Reynolds number are shown in Fig. 3, with drag 
values of a typical production airfoil and a highly 
polished airfoil? shown for comparison. The F-94 
airplane referenced in Fig. 3 was equipped with a 13 
percent glove airfoil over a portion of the F-94 
wing. Suction was applied to the upper glove sur- 
face through 69 slots. Full chord laminar flow was 


¢ Demonstrated drag reduction 


LFC is feasible from the following standpoints: 


¢ Practicality of construction and maintenance 


¢ Possibility of lower operating costs in future subsonic transports 


Fig. 1. LFC definition. 


maintained up to 36 million Reynolds number 
the test limit of the airplane. Over 200 flights were 
made with this airplane.’ 

The swept wing curve of Fig. 3 indicates the feasi- 
bility of LFC on a modern jet transport airplane. 
This curve is a result of tests that were conducted in 
the University of Michigan 5 by 7 ft tunnel and 
the Ames 12-ft tunnel on a wing of 30° sweep, with a 
modified 66-012 airfoil. As compared with a typical 
turbulent swept wing, the drag of the laminar wing 
is five to six times lower.* 

A necessary part of achieving the low friction 
drags associated with LFC is a smooth airfoil sur- 
face. The smoothness required for the swept wing 
of a typical high subsonic airplane is shown as a 
function of altitude® ° in Fig. 4. 

To lend some physical feeling to the numbers 
quoted in Fig. 4, a few comparisons can be made. 
The roughness of 0.005 in. is approximately equiva- 
lent to fine sandpaper. This requirement is easily 
bettered by conventional sheet aluminum. Skin 
joints should be designed such that possible misalign- 
ments result in “up steps” instead of ‘‘down steps” 
to capitalize on maximum allowable step height. 
The ‘‘up step” allowable dimension of 0.017 in. at 
40,000-ft altitude is approximately equivalent to 
the thickness of a playing card. This tolerance can 
be met in most areas with current manufacturing 
methods, although special attention is required 
where the leading and trail- (Continued on page 44) 
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Vehicle Design 


Some Notes on Aerospace Vehicle 
Design Trends 


Adolph Burstein, AFIAS 


The author is a member of the IAS Aerospace Tech- 
nology Panel on Vehicle Design, and Assistant 
Chief Engineer, Research and Development, General 
Dynamics /Convair. 


i GENERAL ACCEPTANCE of jet transportation by 
the airlines and the public has been most gratifying 
and lends emphasis to several design trends. 


Jet Transports—Present and Future 


Subsonic Jet 


There is a fast-growing use of jets on relatively 
short routes involving operations from limited fields. 
This leads to more extensive incorporation of high- 
lift and leading-edge devices and more sophisti- 
cated trailing-edge flaps. 

The turbofan engines answer the need for both 
higher unit thrusts and greater fuel economy. 
Higher thrust engines coupled with lower take-off 
speeds made possible as a result of higher lift 
capabilities of the wing increase one-engine-out 
control difficulty. To solve this problem, several 
rearrangements of power plants from the current 
individual wing-mounted pod-pylon systems have 
been either used or proposed. The two most 
prominent approaches are tail-mounted or twin- 
pod, wing-mounted engines. 

The greater thrust adds emphasis to the general 
desirability of increasing directional control and 
stability. Powered rudder systems were introduced 
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—as seen hy the panels 


because they offered greater control effectiveness 
than servotab actuated systems, while permitting 
low control forces. Another major advantage of the 


irreversible control system is that the rudder 


position is unaffected by aircraft yaw angle. 
The added complexity of powered controls had 


to be accepted to provide the desired effectiveness. 


Supersonic Transport 


Supersonic transport studies have indicated that a 
Mach 3 airliner carrying in the neighborhood of 130 
passengers and designed for a range of some 3,500 
miles could be operated economically in the 1970 
time period. Since fewer aircraft will be required 
than in the subsonic jet fleet, the high development 
costs would make a privately financed venture 
highly unlikely. 

Normal extensions of current technology would be 
required to produce a Mach 3 transport. A few 
illustrative problems and current conclusions are 
presented. Dry turbojets or turbofan engines appear 
to be most promising for operation in the vicinity 
of Mach 38. Sonic boom may be entering the 
power-plant selection considerations. Sonic boom 
is a major operational problem which may require 
| be restricted to 
Turbofan cycle would be 


that acceleration beyond M = 
altitudes above 40,000 ft. 
favored for such operation. 

Obviously important in supersonic transport 
design is L/D. Two items of interest will be men- 
tioned here. At supersonic speeds, the volume of the 
fuselage has much more influence on the total drag of 
the vehicle than in the subsonic case. To illustrate: 
reducing the fuselage volume standards from that 
used in subsonic jets to the Caravelle or DC-7 
class will produce a direct operating cost reduction 
in the neighborhood of 15 percent. This certainly 
calls for a serious discussion of the luxury that 
needs to be provided for the supersonic transport 


passengers. Another drag item that we wish to 
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single out for comment is trim drag. Because of the 
center of pressure shift between subsonic and 
supersonic flight, methods must be found to min- 
imize the penalty of trim drag during the supersonic 
cruise. Canard control surfaces offer interesting 
possibilities of providing efficient methods of achiev- 
ing trim and desirable stability variation between 
subsonic and supersonic regimes. 

Efficient structural design and selection of ma- 
terials is important to minimize weight and cost 
penalties caused by the high speed. Both titanium 
and steel have to be considered. Cabin insulation, 
as well as wing tank insulation, are important prob- 
lems. Much concern has been expressed about 
fuel temperature rise at the end of a Mach 3 mission. 

Some operational problems should be mentioned. 
It is expected that supersonic transports will require 
field lengths similar to the current subsonic jets. 
Therefore, no special problem exists there. How- 
ever, sonic boom will require major consideration 
and may affect airframe design as already mentioned. 
The major problem is the lack of supersonic operating 
experience. The only vehicle that has significant 
supersonic endurance, at the present time, is the 
B-58. This is compared to the thousands of hours 
that were accumulated on large, high-speed, sub- 
sonic, military jets prior to inauguration of pas- 
senger jet service. There is, and must be, serious 
planning to take full advantage of available means to 
gather maximum experience in preparation for 
supersonic passenger traffic. It is desirable that a 
program be instituted which will permit the develop- 
ment of safe operating procedures and the regulations 
to ensure such procedures. 


Exoatmospheric Transports— The Aerospaceplane 


This discussion is not complete without a brief 
look beyond the earth’s atmosphere where it is 
possible to visualize a number of missions where 
men and material must be placed in earth orbits 
simply and economically. Such transports must be 
manned and therefore recoverable. They should 
approach the operational flexibility and safety of the 
transports discussed earlier. The various inves- 
tigations show that this can be achieved in one 
stage by taking full advantage of the air-breathing 
engines to as high a speed as efficiency considerations 
permit, switching at that point to rockets. Right 
here we have a broad field for ingenuity and re- 
search in the area of air-breathing engine cycles. 
But propulsion is by no means the only problem 
area. In selecting the configuration for the aero- 
Spaceplane, the designer is faced with a set of 
formidable requirements, some of them conflicting. 
Aerodynamically, we must have an efficient, control- 
lable configuration throughout the speed range 
from low subsonic to hypersonic re-entry velocities. 
We must take advantage of the pressure field under- 


neath the wing or body for our air inlets; at the 
same time, we have to consider the increased local 
heating that will occur around the intersection of 
two surfaces formed by the intersection of a cowl and 
wing or body. 

Structural heating and insulation problems are 
well known. Since liquid hydrogen is apt to be the 
fuel for this vehicle, we shall have a structure 
exposed to cryogenic temperatures on the inside and 
to over 2,000°F on the outside. (Higher tempera- 
tures would be experienced at the leading edges.) 

Human factors, communication, and navigation, 
all present additional challenges. 

A successful vehicle will indeed call for a most 
thorough and imaginative synthesis of the aerospace 
technologies to achieve the optimum combination. 

Two-stage recoverable versions are interesting be- 
cause of much greater payload potential and there- 
fore place fewer demands on technology advances re- 
quired to produce a useful vehicle. Of course, defi- 
nite sacrifice in operational flexibility would be in- 
volved in going from single to two-stage design. 


Flight Mechanics 


Reliability—One of Today’s Dominant 
Problems 


C. B. Westbrook, MIAS 


Mr. Westbrook is a member of the IAS Aerospace 
Technology Panel on Flight Mechanics, and Chief 
of the Aero Space Mechanics Branch of the Flight 
Control Laboratory, USAF, at Wright-Patterson 
AFB, Ohio. 


= DOUBT, one of the dominant problems of 
the design of future aerospace vehicles is that of 
reliability. It is a subject, recently become fashion- 
able, which, like the weather, many people are 
talking about. Also like the weather, a much smaller 
number are doing anything very constructive about 
it. 

One of the problems in this respect is that not 
everyone is talking the same language or about the 
same missions and vehicles. For example, one 
person may be thinking of the reliability problem of 
operating the B-58 in service with a reasonable 
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utilization. Another may be thinking of the prob- 
lem of reliable firing of an Atlas or Titan. A 
third may be concerned with reliable operation of an 
unmanned space vehicle for a year or more. There 
is a vast difference between the extent of the problem 
and the methods and procedures used to combat 
the problem for these various missions. 

Some groups appear to believe that reliability 
can be “tested” into a system. While gradual 
and sometimes significant improvements can be 
made in this way by elimination of “‘bugs,”’ it is 
certainly not the answer to the general problem. 

Other groups appear to pin their faith on building 
the various components to the very highest quality 
and/or selecting high quality components from 
production. This approach is necessary and de- 
sirable in many cases where key elements must be 
in the system and the best is none too good. 

It is disconcerting to note, at times, the emphasis 
placed on quality and reliability of certain classes 
of components when the overall gain in system 
reliability may not warrant it. For example, 
with regard to the flight control subsystem, relatively 
little effort has been placed on the electromechanical 
parts of the system as contrasted with the electronic 
portions, and yet a chain is only as strong as its 
weakest link. 

The key factor in proper understanding and then 
design is the ability to analyze a proposed system 
early in the proposal stage. With this knowledge 
available, an engineer can identify the critical areas 
and plan what must be done. His first action no 
doubt would be to review his design to see if these 
critical points can be avoided by their outright 
elimination, redesign of the system to a simpler 
one, or substitution of other kinds of equipment. 
If a critical item is found to be indispensable, a 
proper course of action can be determined—develop- 
ment of a “‘super’’ component, use of redundancy, 
etc. Many times, unfortunately, flight control 
system designers accept requests for additions, frills, 
and niceties to the system as a challenge to their 
ability to perform the operation, with little regard 
to its effect on system reliability. Designers must 
learn to refuse such requests resolutely, based on 
concrete knowledge of what it will do to the system. 

The mathematics of reliability analysis is relatively 
simple and straightforward. As with any analysis, 
however, the result is only as good as the data put 
into it. Efforts are being made to obtain these data 
particularly on electronic components. Industry 
has recognized the need enough to put aside its 
normal proprietary interest and share information 
under the IDEP program. In certain areas, how- 
ever, our knowledge is weak or totally inadequate. 
For example, data and knowledge of how systems 
operate when turned on and off are practically non- 

(Continued on page 94) 
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Relating Industrial Research to 
Company Goals 


Lewis A. Rodert, FIAS 


Mr. Rodert is a member of the IAS Aerospace Tech- 
nology Panel on Management, and Research 
Assistant to the Chief Scientist of Lockheed-California 
Co. 


— THE long-range planning problems the 
aerospace industry faces, none is more pressing 
than that involving research. Industry managers 
require better knowledge of how to return a profit 
on research investments. Research scientists and 
engineers want to know how to make recognized 
achievements in the industrial environment. 

Ambitions, however, are not limited to the narrow 
point of view. Managers want their scientists to 
receive recognition for their achievements, and 
scientists want to see their ideas turned into profit. 
Yet, frustration beclouds many who hold to these 
objectives. The aerospace industry has problems to 
solve in harnessing the research function to com- 
pany goals. Greater attention to long-range 
planning of industrial research arises directly from 
apprehensions that the research effort might fail to 
advance national scientific strength and promote 
industrial growth. The impact on the nation’s 
strength and the dollar magnitude of industrial 
research! eliminate the possibility of ignoring the 
problem. Its complication®* precludes solutions 
found quickly by casual processes conducted on a 
part-time basis. 

Relating industrial research with company goals 
involves first identifying company goals encompas- 
sing (1) economic objectives, (2) product area plans, 
(3) resources to support the accomplishment of 
objectives, and (4) the environment in which the 
company is to compete.‘ A detailed enumeration of 
these items comprises one of the principal parts of 
relating research to company goals. 

The other side of the balance sheet is concerned 
with how research may be described. Principal 
headings include programing, facilities, organization, 
and finance. An important decision in programing 
involves the selection of technical areas in which the 
company plans to be strong. Furthermore, the 
decision must be made that research will be con- 
ducted either to produce product line innovation or 
to preserve the continuity of the existing product 
lines and maintain technical ability. From thiese 
decisions many subordinate topics must be decided, 

(Continued on page 94) 
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Adleentures 
in Engineering 


James D. 
Mason 


a SENSITIVITY of germanium transistors 
was discovered in the development of a guidance 
system for a major ballistic missile program. This 
guidance system uses approximately 1,000 ger- 
manium transistors. Cooling air is provided near 
0°C prior to turn-on. Transistor malfunction 
occurred in approximately 0.1 percent of a given 
transistor population immediately after turn-on. 
The malfunction usually disappeared after a short 
warmup period, and was caused by a drastic ex- 
cursion of Icpo over a narrow temperature range. 
This problem was assigned to the writer and a task 
force of engineers was provided to help effect the 
solution. 

The transistors had been purchased to a specifica- 
tion which included 100 percent leak testing and 
elevated temperature test. 

When transistor temperature sensitivity appeared, 
an extensive program of failure analysis and test was 
indicated. As an interim step, all transistors of this 
type were subjected to a series of temperature cycles 
from —35 to +85°C with initial and post tempera- 
ture cycling readouts of Icpo and Inpo. Any tran- 
sistor exhibiting a change in excess of 2uA was re- 
jected. This screening procedure was established 
deductively in an attempt to improve stocks and to 
require vendors to institute more stringent manu- 
facturing controls. 

Test results indicated that moisture contamina- 
tion was the cause of difficulty. Vendors’ manu- 
facturing and assembly techniques indicated that 
moisture content was within limits when manu- 
factured. Study indicated that transistor hermetic 
seals were adequate to exclude moisture from ex- 
ternal sources for at least a period of time which 
exceeded the age of the transistors under study. 
Further investigation indicated that the internal 
structure of the transistor was outgassing under 


Temperature Sensitivity in 
Germanium Transistors 


The author is Supervisor, Piece Parts Reliability Standards Group, Product 
Reliability Department, Arma Division, American Bosch Arma Corpora- 
tion. He received his B.E.E. degree from Syracuse University in 1949, 
joining Crouse-Hinds Co. as Development Engineer and Laboratory Super- 
visor. On becoming associated with Arma in 1957, he was assigned 
responsibilities in the Engineering Standards Section as a Semiconductor 
Specialist, shortly thereafter being assigned to his present position. Mr. 
Mason is a member of the IRE and the Standards Engineers Society. 


storage and use, increasing the moisture content of 
the atmosphere within the transistor housing. 

Dew point testing was investigated as a possible 
screening mechanism. (The moisture level in the 
internal atmosphere of a transistor should be such 
as to exhibit a dew point below — 65°C.) 

Dew point testing was discarded as a screening or 
weeding mechanism for several reasons: 

(a) Individual handling was required. 

(b) A subjective decision was required. 

(c) Characteristic anomalies, when present, were 
displayed very briefly. 

(d) Test method caused operator fatigue. 

(e) Precautions were required to avoid false indica- 
tions from moisture condensation on the exterior 
of the transistor housing. 

It became obvious that further work had to 
be done to evaluate, within a reasonable time, 
the efficiency of the temperature cycling screening 
mentioned earlier. Furthermore, immediate atten- 
tion had to be given to recent developments in tran- 
sistor construction. 

The most noteworthy change in transistor con- 
struction at the time of investigation was the advent 
of the ‘‘gettered”’ transistor. A “getter” in a tran- 
sistor consists of the addition of a molecular sieve 
material to the interior of the transistor. The mo- 
lecular sieve mechanically traps moisture molecules 
and holds them within its structure over a wide 
temperature range. A test series was established 
using a temperature plateau test. This test con- 
sisted of exposing groups of transistors to tempera- 
tures in 10°C steps from —55 to +25°C with elec- 
trical measurements at each step or “plateau.” 
Groups tested included: 

(a) Transistors accepted via temperature cycling 
tests. 

(Continued on page 95) 
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Piloted Simulation Studies (Continued from page 11) 


Fig. 2. 
in simulation of prediction-control method. 


Fixed cockpit and analog computer used 


Methods of Control 


Prediction Control 


Basic concept: A brief description of 
the prediction method of end-point 
guidance and control is presented here 
for clarity. (Full details of this system 
are given in reference 1.) The aero- 
dynamic range capability of the vehicle 
is predicted from existing conditions 
along the flight path by solving equa- 
tions of motion for the constant trim 
conditions for maximum manuever. 
In this study the three constant trim 
conditions used were: 


(L/D)maz = 0.5 and bank angle equal 
to zero (¢ = O°) for 
maximum range 

(L/D)min = O and = 0° for minimum 
range 

(L/D)mazr = 0.5 and g = 45° for maxi- 
mum crossrange 


This area of aerodynamic range 
(shaded in Fig. 1) is solved repetitively, 
as ‘‘fast” as possible, for changing flight 
conditions by an airborne computer; 
the range and maximum manuever 
capability are thus known continuously 
as the vehicle parameters vary. 

One of the major problems in a con- 
tinuous prediction system is to choose 


INITIAL CONDITIONS 


r == 
VISUAL CONTROL 
cues INPUTS | VEHICLE 6-DEGREE 


MOTION EQUATIONS 


Fig. 3. Block diagram of analog simulation of 
prediction-control method. 
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the form of the equations of motion 
used in the airborne prediction com- 
puter so that a satisfactory compromise 
is achieved between the conflicting re- 
quirements for a reasonable amount of 
computing equipment, accuracy, speed, 
and easily obtainable input information. 
Chapman! has derived a mathematical 
model which defines the important 
variables for re-entry trajectories. A 
comparison is made below of the quanti- 
ties needed to describe a trajectory for 
the standard equations of motion and 
those needed for the Chapman model. 


Equation Airborne Vehicle 
Measurements Parameters 
Standard V, yh mCpS, L/D 
Chapman V, vy, Aw L/D 
BOUNDARY FOR SKIP OUT OF ATMOSPHERE 
@ROLL ANGLE DESTINATION 
BOX 
4s" xf 
60” 
ROLL INDICATOR i TRIM INDICATOR 
BOUNDARY FOR LIMIT OF 10g DECELERATION 


Fig. 4. Pilot's guidance and control display. 


The Chapman model combines the 
vehicle’s mass-to-drag ratio m/CpS and 
altitude measurement into a decelera- 
tion measurement A, which gives an 
effective ‘“‘density altitude.”” This den- 
sity altitude is then related to an ex- 
ponential altitude density variation by 
the Chapman equation 

There are three factors which make 
the Chapman equation satisfactory 
for use in the airborne prediction com- 
puter. The first factor is that the com- 
putations are simplified by the use of 
acceleration (density-altitude) inputs. 
The second factor is that the relation 
between density and geometric altitude 
will usually be unknown during an 
actual re-entry; therefore, use of the 
effective density-altitude variation will 
provide a way of compensating for 
nonstandard atmosphere variation. 
The third factor is that the Chapman 
model does not use geometric altitude. 
White® shows large errors may be pres- 
ent in the measurement of geometric alti- 
tude during re-entry if an inertial plat- 
form is used for position information. 
Computations which use acceleration 
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measurement instead of geometric altj- 
tude will be less affected by inherent iner- 
tial platformerror. It is recognized that 
the continuous-prediction method must 
still use an inertial platform to measure 
longitudinal and lateral range. How- 
ever, a cursory examination indicates 
this information would be available 
with the required accuracy. 

Simulation details: The prediction 
method of end-point guidance and con- 
trol was simulated by use of an elec- 
tronic analog computer and a fixed 
cockpit (Fig. 2). The elements of the 
analog simulation important to the 
navigation study are shown in Fig, 3. 
The analog computer functioned to (a) 
accept control inputs about all three 
axes from the pilot and then solve the 
resulting six-degree-of-freedom motions 
of the vehicle, (b) compute the repeti- 
tive prediction equation, and (c) gen- 
erate the pilot display information. 
Complete details of this simulation are 
given in reference 1, but for clarity a 
short discussion of the simulation will 
be given here. 

The vehicle equations of motion de- 
scribe six-degree-of-freedom dynamics 
along a great circle route (with small 
lateral displacements) over a nonrotat- 
ing curved earth. The flight regime 
simulated was from 100,000 to 300,000 
ft variation in altitude using the 1956 
ARDC standard atmosphere, and from 
rest to 35,000 fps variation in velocity. 
The vehicle simulated had proportional 
reaction controls about all three axes, 
with constant gain damping augmenta- 
tion included so that the short-period 
damping would be near-critical at the 
highest dynamic pressures. This damp- 
ing along with very little cross coupling 
in the aerodynamic derivatives made 
the vehicle easy to fly in the short- 
period mode so that the pilot’s function 
was mainly one of guidance control. 
The vehicle’s aerodynamic character- 
istics, the re-entry altitude and velocity, 
and the allowable deceleration excur- 
sions imply®~’ safe re-entry flight-path 
angles of —2.1° to —6.2°. 

The analog computer used in this 
simulation was not equipped with high- 
speed computer components. Because 


300,000 SKIP TO THE EDGE OF THE 
ATMOSPHERE FOLLOWED BY 
us TRIM OF MAXIMUM L/D 
200,000 

100,000 
FLOWN TO 10g LIMIT FOLLOWED 
a BY TRIM AT L/D=0 
0 1 1 

= 500 

S o 

8 

8 500 


' 
° 500 1000 1500 2000 2500 
DOWNRANGE FROM ENTRANCE, MILES 


Fig. 5. Extremes of range possible using pre 
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of this, the repetitive solution of the 
prediction equation could not be con- 
sidered fast. The integration rate was 
dy/dt = 0.2/see to which was added a 
|-see reset time. This means the solu- 
tion from parabolic velocity was 8 
sec, from satellite velocity 6 sec, and 
from one-half satellite velocity 3.5 sec. 
It can be seen that the repetition rate 
became more rapid as the vehicle neared 
its destination. One item of interest 
left unanswered by this simulation is 
to determine the effects of a truly fast 
solution time, particularly in the region 
in which pilot control was most sen- 
sitive. (This sensitive control region 
will be described later in this paper.) 

The guidance display used in this 
simulation is shown in Fig. 4. The 
footprint”’ or range envelope is a trans- 
parent overlay on the face of a 5-in. 
oscilloscope. The destination is a mov- 
ing dot on the pilot’s display oscilloscope 
presented in relation to the predicted 
range capabilities of the vehicle in the 
following manner: 


yoscilloscope input = Rzq — Rtmin 
Rimes Romina 


xoscilloscope input = Rug 


Rumaz 


The correspondence then between 
the oscilloscope presentation for the 
pilot (Fig. 4) and the aerodynamic 
range capability of the vehicle (Fig. 1) 
can readily be seen. In fact, the oscil- 
loscope display appears to the pilot as 
though he were looking through a win- 
dow in the vehicle as it approaches its 
destination. The pilot’s field of view 
through the window represents the 
maneuver capability of the vehicle. 
Then it is good procedure for the pilot 
to control so that his destination is in 
the center of his maneuver capability— 
that is, to fly the destination dot to the 
middle of the oscilloscope. 

In this simulation it was decided to 
impose a maximum deceleration limit 
of 10g on the vehicle’s trajectory as the 
arbitrary physiological acceleration limit 
on the pilot. It was also decided to 
impose a minimum deceleration limit 
of 0.075g to keep the vehicle from 
skipping out of the atmosphere. The 
effect of these g limits on the range ca- 
pability of a vehicle entering at parabolic 
velocity is illustrated in Fig. 1. The 
repetitive prediction equations of range 
capability also predict the decelerations 
the vehicle will experience for fixed 
L/D. If the predicted value of de- 
celeration for L/D = 0 exceeded the 
10g limit, then an empirical relation- 
ship taken from references 5 and 8 was 
used to indicate on the oscilloscope 
(Fig. 4) the minimum L/D that could 
be held without exceeding the 10g limit. 
If the predicted value of deceleration for 
L/D = 0.5 was less than 0.075g (equiv- 
alent to a skip-out of the atmosphere), 


then the next repetitive maximum 
range solution included a value of 
A(L/D) proportional to the amount of 
overshoot, and this new L/D maximum 
value is displayed on the pilot’s display 
(Fig. 4). Thus, the repetitive calcula- 
tion is continually presenting to the 
pilot an L/D that if flown would cause 
the vehicle to reach the 10g decelera- 
tion limit, and another L/D that if flown 
would cause the vehicle to skip just 
to the edge of the atmosphere. The 
vehicle range capabilities could con- 
ceivably be extended beyond those of 
this simulation by allowing the vehicle 
to skip out of the atmosphere for added 
range. 

Pilot technique and opinion: The 
method of pilot control for maximum 
or minimum range is illustrated in 
Fig. 5. For maximum range the pilot 
will trim to the L/D shown to be the 
boundary for skip-out of the atmosphere. 
The parabolic entry velocity will even- 
tually diminish, and as it does the skip- 
out boundary will change until the 
pilot finds he is permitted to fly an 
L/D of 0.5. Conversely, to obtain 
minimum range, the pilot would trim 
to the L/D indicated by the boundary 
for the 10g deceleration limit and, as the 
vehicle velocity diminishes, the bound- 
ary will change until eventually the 
pilot can fly a trim of L/D = 0. Fig. 
5 also shows the range area attainable 
by this vehicle when it is flown by this 
guidance method. 

The pilot control technique used to 
fly to a desired end point is as follows: 
During the first part of the trajectory 
when the vehicle is at parabolic velocity 
and just entering the atmosphere, the 
pilot trims at those conditions which 
the guidance display indicates are neces- 
sary to reach the destination. If the 
destination mark on the guidance dis- 
play should appear below the maximum 
deceleration boundary, the pilot just 
trims to that boundary. After the 
initial dive into the atmosphere and 
the point of maximum deceleration is 
passed, the pilot overcontrols slightly 
the destination to the center of the 
vehicle’s range capability. 

The pilots considered this system 
satisfactory. The pilot function needed 


here was to close the control loop be- 
tween the bank and trim indicators and | Graphic Reproduction Sales Division | 
navigation display (Fig. 4), which in- | Se Corer | 
volved keeping up with the changing | Please send me a copy of the booklet, | 
trim conditions during the re-entry | “Etching, Chemical Milling and Plating | 
trajectory. | with Kodak Metal-Etch Resist’’ (P-36). | 

Pilots had no difficulty coping with | ome | 
emergency or abort re-entry conditions lc | 
while using the continuous-prediction en” 
guidance and control system. Fig. 6 | Street 
shows the piloted range capability for | City 
two typical abort trajectories. This | Zone State 
ability of the guidance system to pre- | 
dict range and deceleration limits for | 
any entry conditions makes it an almost | 
universal control method. Eee 
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NEW KODAK 


METAL-ETCH 
RESIST 


opens new chemical-milling 
applications 


Speeds up, simplifies deep etch 
weight reduction and parts manu- 
facture. Reproduces fine-line detail 
as in plating, dial and name-plate 
making. This new photographic 
process ends time-consuming hand- 
work, results in high accuracy. 
Kodak Metal-Etch Resist withstands 
acids, alkalies, electrolytic fluids 

. adheres well to aluminum, 
titanium, magnesium, stainless and 
other alloy steels. High stability 
and strict uniformity simplify vol- 
ume production. Send today for a 
detailed 16-page brochure that 
gives all the facts. 


Text for this advertisement was set photographically. 


No statement or suggestion is to be considered a recommenda 
tion or inducement of any use, manufacture or sale that may 
infringe on any patents now or hereafter in existence 
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800,000 


600,000 
ALTITUDE, FT 


REENTRY 
400,000 V=30,500 FPS 


Y¥=-2.75 DEG 


REENTRY 
V*22,600 FPS 
200,000 DEG 


CROSSRANGE, 
MILES 


1000 2000 3000 4000 5000 
DOWNRANGE DISTANCE FROM 
LAUNCH POINT, MILES 


Fig. 6. Maneuver capabilities for typical abort 
re-entries. 
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Fig. 7. Downrange and corridor limits for 
nearly parabolic re-entry. 
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Block diagram of a piloted fixed- 
trajectory control method. 
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Fig. 9. Block diagram of analog simulation of 
fixed-trajectory control method. 
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Sensitive control regions: Flying those 
trajectories near the maximum decelera- 
tion boundary as indicated on the guid- 
ance display was not difficult. The 
procedure was to hold trim conditions 
just above the maximum deceleration 
boundary until the peak deceleration 
had passed and then to hold the trim 
conditions necessary to reach the des- 
tination. Flying to those destinations 
which required a skip to the upper limit 
of the atmosphere (300,000 ft for this 
simulation) to obtain the necessary 
range resulted in a sensitive control 
problem. For a maximum range tra- 
jectory (Fig. 5 is an example), little 
control effectiveness was available near 
the top of the skip (at maximum 
altitude), and the trajectory control 
exercised during the first dive into the 
atmosphere strongly influenced the 
rest of the trajectory. Fig. 7 is pre- 
sented to show those values of range and 
entrance flight-path angle (h = 290,000 
ft, V = 35,000 fps, and (L/Dmaz = 
0.5 for all runs) for which pilot control 
was most sensitive—that is, the shaded 
area indicates entry conditions for 
which small changes in L/D control at 


35,000 
30000 ONTROLLED TRAJECTORY 


25000 - FIXED TRAJECTORY ~ 


FPS 


20000 


15,000 


VELOCITY, 


10,000 
5000 


2500 2000 ~ 1500 000 500. O 
RANGE TO DESTINATION, MILES 


Fig. 10. Typical control to a fixed trajectory. 


the beginning of re-entry results in very 
large changes in range. 

Fig. 7 is a good index of two major 
flying qualities problems of flight within 
the atmosphere at near parabolic veloc- 
ities. The region of entry conditions 
acceptable to the pilot is squeezed be- 
tween the sensitive control region and 
the limit of deceleration the pilot can 
tolerate with current restraint harness 
and control techniques. Research on 
both of these pilot-control problems 
is described in references 10 and 11. 
Further, Fig. 7 shows not only sensitive 
control regions, but also the range of 
readily controllable entry angles as 
related to maximum and minimum 
range obtainable by this guidance 
method. 


Fixed-Trajectory Control 


Method of analysis: Fixed-trajectory 
control is conceived to be either a 
piloted technique (Fig. 8), which is 
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analogous to the system for the cop 
tinuous-prediction method, or possibly 
an automatic control method (Fig. 9) 
in which vehicle flight-path parameters 
are compared to the parameters of some 
precomputed design trajectory. Jy 
either case, the pilot or automatig 
controller simply initiates vehicle cop 
trol inputs to null any error betweg 
the two. For purposes of direct cog 
parison it would have been desirabjg 
to study the system as represented # 
Fig. 8; however, to expedite the study, 
the Chapman equations’ were agaj 
used as the re-entry equations. 
system as represented in Fig. 9 
then used in a closed-loop analysi 
since the independent variable in thes 
equations is velocity and not time, the 
human pilot could not be introduced, 
The particular comparisons in this 
paper are believed to be unaffected singe 
the function of a human pilot in closing 
the control loop was closely approx 
imated by the analytical method. 

Fig. 9 shows the complete block dia 
gram of the closed-loop method of 
analyzing this guidance system. The 
primary objective of a fixed-trajectory 
guidance system is to find a combination 
of predetermined trajectory and feed- 
back control that will minimize the 
maximum range error at the end poifi 
from large variations in initial com 
ditions. Typical control to a fixed 
trajectory is shown in Fig. 10. The 
outer loop or range control of Fig, § 
will cause the vehicle trajectory to cot 
verge to the design trajectory (Fig. 10). 
The inner-loop controls of Fig. 9, 
which are the deceleration and the rate 
of-climb feedback, primarily affect the 
overshoots or oscillations about thé 
fixed trajectory. The analysis of thi 
system is made with the same limi 
tions on the vehicle and its trajectorie® 
that were imposed on the predictidif 
control method; (L/D)maz = 
deceleration limits of 10g maximum ang 
0.075g minimum, and entry condition® 
of h = 290,000 ft and V = 35,000 fps 
These conditions imply a safe re-entty 
corridor —2.1° to —6.2° flight-patlt 
angle. 

Simulation details: The use of the 
Chapman equation as the model for the 
vehicle trajectory equation of motiol 
allows the important feedback quantities 
for range control to be recognized, and 
its relative ease of mechanization on all 
analog computer permits a straight 
forward method of investigation of fixed 
trajectory control. The control 
range by the L/D input is essentially 
a third-order control system. The 
inner loop (rate of climb and deceleré 
tion feedback) which contribute damp 
ing by the control system can be see# 
in Fig. 9 to come from the Z’ and Z 
terms. It is interesting to note the 
feedback loop from Z’ and Z within the 
Chapman equation. The upper lop 
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He’s solving 
a real estate 
problem 


This AMF engineer’s job is deter- 
mining how best to move big missiles 
off shore for launching. Should they 
be floated out horizontally, flooded 
to an upright position, and then 
launched? Or, would it be more fea- 
sible to barge them out? Might they 
be moved to or assembled on “Texas 
Towers,” or would a causeway or 
simply land-fill be the answer? 


Behind the project is our shrinking 
real estate at launching sites, plus the 
‘ hazards inherent in launching Sat- 
urn-sized missiles (and the coming, 
nuclear-powered missiles) near other 
installations. Off-shore launching 
; may be the answer. 


Feasibility studies of all types are 
{ an AMF specialty. What kind of 
remotely controlled machinery is 
required to service nuclear-powered 
aircraft? What kind of habitation 
could be built on (or under) the sur- 
face of the moon? What sort of 
machines (manned and unmanned) 
could survey the moon’s surface 
without, for example, falling into a 
: fissure? What is the best way to 
assemble a space station? All these 
are problems AMF engineers are 
presently investigating. 


lt cil If your problem is the first of its 
ee. kind, AMF will not, of course, have 
met it before. But AMF’s long expe- 
rience in accepting totally unique 
- challenges gives it an advantage 
. enjoyed by few other organizations 


concerned with ground support, 

launchability and space environ- 

eS ment. To get further information 

write American Machine & Foundry 

‘ Company, 261 Madison Avenue, 

New York 16, N. Y. 


Bineering and.manufacturing AMF has ingenuity you can use... 
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PREDICTION 
CONTROL METHOD 


2 FIXED TRAJECTORY 
CONTROL METHOD 
500 1000 1500 2000 2500 
DOWNRANGE FROM ENTRANCE, MILES 
Fig. 11. Extremes of range possible using 
fixed-trajectory control method. 
w 3000; 
UNDERSHOOT 
= OVERSHOOT LIMITS 
ys 6g 10g 15g 
& 2000 + 
= 1500} ee 
w 1000F 
z 
<a 
= 
= 


| 2 3 6 5 6 7 8 
ENTRANCE FLIGHT PATH ANGLE, DEG 


Fig. 12. Effect of maximum deceleration 
limit on range andentry corridor; fixed-tra- 
jectory method. 


(1—*~)/Z is a stabilizing quantity 
for a less than one, but for @ greater 
than one this loop is destabilizing. The 
lower loop (4Z’—Z) within the Chap- 
man equation is a stabilizing loop that 
is essentially a rate-of-climb feedback 
loop that tends inherently to damp 
phugoid oscillation about a design 
trajectory. 

The simulation of this control sys- 
tem has been done on an electronic 
analog computer in which % was the 
independent variable. It was assumed 
that the response times of control in- 
puts would not greatly affect the over- 
all trajectory motion. The constant- 
gain inner-loop quantities were taken 
as those measurements which may be 
available during a re-entry. 

Crossrange control was simulated 
by a second-order system. The equation 
for crossrange was taken from Slye’? 


These desired L/D and Y/D values 
are related to the vehicle control inputs, 
roll angle y, and trim angle of attack 
(L/D), in the following manner: 


L/D = [(L/D), cos ¢ 
Y/D = [(L /D) sin 


This coupling was usually neglected 
in the past for high L/D vehicles where 
the desired L/D programs angle of 
attack and Y/D programs roll angle 
(¢ is usually less than 45 This study 
of a low L/D vehicle will consider the 
complete coupling. 

Feedback control analysis: The flight- 
path conditions of a re-entry from the 
center of the safe entry corridor were 
chosen as the initial conditions for the 
design trajectory. The determination 
of such a trajectory for parabolic entry 
presents a problem. In studies of 
satellite re-entries, the design trajec- 
tory has usually been chosen as an 
equilibrium glide for a constant L/D in 
the middle of the vehicle’s lifting ca- 
pability. However, if constant L/D is 
used at parabolic velocities, the tra- 
jectory is unstabie, and, therefore, an 
equilibrium glide condition cannot be 
used to define the design trajectory in 
this case. The trajectory chosen for 
parabolic re-entry was that which results 
from L/D control proportional to rate 
of climb. This simple rate-of-climb 
feedback loop is a very effective means 
of stabilizing the design trajectory. 
The gain of the rate-of-climb feedback 
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‘ Fig. 13. Comparison of range and entry corri- 
and was expressed in terms of the head- dor limits. 
ing angle y and the Chapman Z func- 
tion. The present analysis considers 
only the control necessary to null the 
crossrange error to zero. tall 
w 
The equations that result from the 3 PREDICTION 
analysis of longitudinal and lateral wl 1500- 
a 
range control are of the following form: é 
z 
= 1000- 
L/D = Kh — hy) + KA — Ay) ne 
‘ FIXED TRAJECTORY 
inner-loop Zz CONTROL METHOD 
feedback 
<= 
= 
+ K,(R., — outer-loop 
5 10 20 
range control MAXIMUM DECELERATION, g's 
Y/D = K,R crossrange 
wa 8 Fig. 14. Comparison of effects of maximum de- 
control celeration limit on range. 
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Fig. 15. Comparison of typical controlled 
trajectories with nearly maximum range. 
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Fig. 16. Pilot's time tolerance to acceleration, 


was chosen as that constant gain which 
commanded an L/D near the middle 
of the vehicle’s lifting capability (L/D « 
0.3) during the suborbital velocity 
portion of the trajectory. 

The control of errors in range from 
the design trajectory was tried using 
only the rate-of-climb feedback for 
damping and range-to-go feedback ior 
primary control. The results were ur- 
satisfactory since the best control that 
could be expected gave end-point errors 
of +25 miles for initial errors of +100 
miles. 

Therefore, a third control loop, ae- 
celeration feedback, was added to the 
system to make the range control ac- 
ceptable. The addition of this contra 
loop allowed higher gains on the range 
control feedback for closer range con- 
trol. The total system then gave end- 
point errors of +5 miles for entry range 
errors from +100 to —500 miles. The 
acceptable control equation developed 
in this preliminary study is 


L/D+*5 = + KA — A;) + 
—0.5 

Ki(R;, 

where K, = 0.001/ft/sec, Ky = —0.33/s 


and K; = 0.006/mile. 

This system could be improved by 
programing optimum variations of K 
K», and K; as functions of velocity i- 
stead of allowing these gains to bt 
constant. For instance, if Ks; is pre 
gramed to zero until just below satellite 
velocity, then the available range c 
pability would be about 2,00\) miles 
compared with 600 miles for 
stant gain system in this report 
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Guidance capabilities: The ground 
yea attainable using the fixed-tra- 
. method is shown in Fig. 11. 
This range capability is for control 
about a design trajectory having an 
initial entry angle of —4.5° and should 
e compared with the prediction con- 
trol range capability from Fig. 5. 

The effect on maximum longitudinal 
range and permissible entry corridor 
‘changing the maximum deceleration 


jimit of the fixed-trajectory method of 
olled ntrol is shown by Fig. 12. In general, 2g 


— 90 


ubling the maximum g limit doubles 
ALTITUDE. . 


the longitudinal range and the entry 
rridor for the deceleration range con- 
sidered in this study. 


_80 


Comparison of Prediction Control and 
Fixed-Trajectory Control COMPREHENSIVE 
The maximum and minimum range 
attainable by the prediction control and INSTRUMENTATION 
the fixed-trajectory control methods are 
mpared in Fig. 13 for various entry 
neles. The figure shows that either 


8 10 
ntrol method can be used to keep from 
ation @ exceeding the peak deceleration limits | HOWELL INsTrRuMENTS are cali- 
075g and 10g for the same interval | brated for use with all transducers for 
fentrv anoles. whic ras —_9° 
fentry angles ; which w is about —2° to the measurement of the phenom on 
-(\° The main difference between in att aj f d 
hich § the two methods is that the fixed- countered in a phases of aircraft an 
iddle F trajectory method of control provides = missile development, test and operation. 
Dx nly about one third of the range at- Each instrument is completely self-con- 
ocity able by the prediction-control sys 0° (miniaturized, 
j 
{rom fhe effect of changing the maximum ized, — driven, hermetically sealed) 
ising F deceleration limit on range capability - with Zener reference, power supply, am- 
lor has been shown for the fixed-trajectory plifier, servo motor, cold junction com- 
c ntrol method in Fig. 12. A com- | pensation as needed and the 144-inch 
- Ul § parison of the effect of the maximum | lide-wi d hed t to 1i . 
that deceleration parameter on the range 
tots § capability for the two control methods TURBINE PRESSURE e.m.f. for exact, counter-type digital read- 
10 


s shown in Fig. 14. The percentage 
change in range capability obtained 
, ac y varying the peak deceleration limit 


out. Needle pointers are provided for 
quick reference. Accuracies are within 
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| 
) the F for the prediction-control method is | oe Available in 3”-dia. MS33639 and 
1 ac slight compared to the variation in | 2”-dia. MS33598. 
nttol § range capability for the fixed-trajectory Styled Series BH183 & 185 — The In- 


ethoc strument with the Tape-Slidewire —these 


COB: Fig. 15 exemplifies the main reason for ; 

end- F these Area a in range capability. instruments are produced by the makers 

range It shows the acceleration-velocity pro- of the JETCAL® Analyzer, the only jet 
The } files for two maximum-range  trajec- engine tester used throughout the world! 

lope F tories for which the initial conditions 


Full information is available for the 


the vehicle and flight path were . 
asking. 


) + — identical. One was guided by prediction 
control and the other by the fixed- 
| “alectory method. The main dif- 
ference between the two methods of 
control is that the fixed-trajectory con- 
trol technique inherently damps_ to 
d by | the design trajectory, whereas the pre- 
K, “etion-control method damps_ only 
ty in- | “© the end point. Therefore, the pre- 
bef “Cuon-control technique spends more 
prof “me at the lower deceleration levels 
elite | ‘Fig. 15), thereby stretching the glide, 
ca | “an does the fixed-trajectory method. 

‘miles Data presented in Fig. 15 correspond 
con toa maximum range extension ata 3479 WEST VICKERY BOULEVARD ° FORT WORTH 7, TEXAS 
Steep initial entry angle condition (close 
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to the limit shown in Fig. 13) for which 
it would be difficult to avoid the peak 
deceleration limit. 

Fig. 16 is a graph, taken from refer- 
ence 11, showing the pilot’s time toler- 
ance to acceleration as presently con- 
ceived for accelerations in the eyeballs- 
out or eyeballs-in direction. The 
dotted line on this graph indicates the 
length of time the pilot is at or above 
the corresponding g level for the two 
parabolic re-entries shown in Fig. 15 
that is, Fig. 16 shows that for either 
guidance method the pilot will be above 
5g for at most one-half minute and above 
8g for at most one-third minute. For 
deceleration levels less than 5g, the 
time tolerance exceeds the total re- 
entry time of the vehicle, and thus is 
of no significance. Thus by either guid- 
ance method the pilot can maneuver 
so as not to exceed his time-tolerance- 
to-acceleration boundaries, and he is 
well within the region where he can 
perform a useful control function. 


Concluding Remarks 


The results of this study of two basic 
methods of terminal guidance and con- 
trol during re-entry from a lunar mis- 
sion show that either method can utilize 
nearly all the permissible re-entry cor- 
ridor fixed g limits, but the range ca- 
pability of the “‘continuous-prediction’’ 


Supercircular Entry and Recovery With Maneuverable Manned Vehicles 


guidance and control method is nearly 
three times the range attuinable by the 
“fixed-trajectory”’ control method. 
Finally, it has been found for entries 
from parabolic orbits that a vehicle can 
be maneuvered by terminal 
guidance technique in such a manner 
that the present criteria for 


either 


pilot’s time 


tolerance to acceleration are not ex- 
ceeded. 
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entry pull-out for gliders was analyzed 
by Grant! and has been discussed by 
others.?~4 

Peak dynamic pressure occurs just 
before the bottom of the pull-out, 
but modulation will generally continue. 
Beyond peak dynamic pressure, angle 
of attack is gradually increased to 
return the glider to maximum lift 
attitude. Throughout this process, lift 
and drag have been maintained at the 
highest level allowed by the prescribed 
acceleration limit, so penetration has 
been minimized and energy dissipation 
(speed reduction) has been maximized. 
This means, generally, that both heating 
rate and time of exposure have been 
held to a minimum so that the total 
heat load up to this point is as small 
as possible for the prescribed conditions. 

The lowest acceleration limit which 
can be prescribed for a given entry 
angle is that which requires modulation 
to proceed from maximum-lift attitude 
at entry to nearly minimum-drag 
attitude at peak dynamic pressure. 
(Heat load in this case is rather large.) 
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selection of a 
requires less 


On the other hand, 
high acceleration limit 
reduction of lift and drag so total 
heat load is reduced. Attitude modula- 
tion allows very steep entries with a 
high L/D glider without exceeding 
10g. But the more significant advan- 
tage of attitude modulation is that 
resultant acceleration loads can be 
limited to more comfortable values at 
moderate entry angles. In either case, 
a heat load design penalty is incurred by 
extensive attitude modulation. A pro- 
spective glider design will involve a 
three-way compromise between design 
heat load, maximum resultant accelera- 
tion, and maximum entry angle or 
undershoot limit. 

When operating at less than the 
maximum entry angle, it is important 
to recognize the different implications 


of heat-load limit and acceleration 
limit. Because of the concern for 
crew comfort, resultant acceleration 


load should always be reduced as far 
below the limit value as_ possible. 
But once a heat load limit is established 
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and provision for it made in the design 
there is no incentive to hold heat loads 
to lower values. For entry angles 
less than the steepest provided for in 
the design, acceleration relief through 
attitude modulation should be used 
to the maximum extent allowed by the 
design heat load limit. 


Pull-Up 


At the bottom of the pull-out from a 
steep entry, the speed, dynamic pressure, 
and heating rate are still very high 
A pull-up maneuver is desired which 
will avoid skip-out, continue contro 
of acceleration, and avoid further 
accumulation of large heat load. Skip 
out is prevented by banking, as required, 
to regulate the vertical component 0 
lift. (In this paper skip-out is dis 
tinguished from overshoot;  skip-ow! 
is defined as an exit from the atmosphere 
caused, in part, by inappropriate use 0! 
lift.) Resultant acceleration load 
controlled by continuing the angle 0 
attack modulation started during the 
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Fig. 3. Manned entry limits. 


pull-out. Maximum drag _ consistent 
with the prescribed acceleration limit 
is maintained, as before, while the 
independent banking control regulates 
the flight path angle. The high drag 
tends to minimize heat load by keeping 
exposure time short. A considerable 
latitude in choice of trajectory is 
generally available. 

The pull-up is controlled to allow 
termination of entry maneuvers above 
satellite speed at an altitude where 
glide with radiation cooling is possible. 
4n inverted pull-over maneuver is 
required to establish level flight. For 
this study the terminal point is at 
27,000 fps and 250,000 ft altitude. 
In the event that speed loss due to 
drag is too rapid to allow completion 
of the maneuver by bank angle control 
alone, the drag may be reduced by 
lowering angle of attack to allow more 
time. Greater heat load results from 
the additional exposure time. 

It is recognized that the pull-up 
maneuver, which is generally completed 
within 3 min after entry, presents a 
difficult control system problem. This 
problem is being studied, but is beyond 
the scope of this paper. Analog com- 
puter studies have verified that the 
maneuver can be successfully controlled 
with the system proposed by Morth and 
opeyer.° 

The small heading change from the 
use of banked flight during pull-up is an 
incidental assist to subsequent banked- 
glide maneuvers to reach the recovery 
site. 


Overshoot Limit 


Very shallow entries at supercircular 
speed can result in an overshoot which 
will make direct single-pass entry 
impossible. This is the result of in- 
sufficient penetration of the atmosphere; 
the glider will re-exit, and there follows 
a multipass series of entries over an 
extended period of time. At the over- 
shoot limit, a pull-in maneuver is re- 
quired in which maximum lift is used, 
but the glider is inverted so that the 
lift is directed toward the earth. This 
maneuver increases penetration to avoid 


overshoot. If there is any likelihood of 
overshoot, it is best to enter inverted. 
Directing the lift toward the earth 
steepens the flight. path slightly, but 
roll to the upright attitude can be 
made, if necessary, after single-pass 
entry is assured and before atmosphere 
penetration is excessive. 

Shallow entries bordering on the 
overshoot limit are completed by an 
inverted glide, within the heat-load 
region, to the terminal point at 27,000 
fps and 250,000 ft altitude. It is 
recognized that constant attitude equi- 
librium glide above satellite speed is not 
stable. The rate of divergence from 
the speed-altitude schedule is low® which 
implies that control would present no 


problem. Entry at the overshoot limit 
involves low resultant acceleration, 
but considerable heat load accrues 


during the supercircular glide within 
the heat load region—i.e., below the 
temperature placard in Fig. 1. 


Entry Corridor 


Entry at an angle between the steep 
(undershoot) limit and the shallow 
(overshoot) limit is said to be within 
the entry corridor. The notion of a 
corridor’ arises from a consideration of 
the virtual or no-atmosphere perigees 
of the undershoot and overshoot limit 
entries. The difference between perigee 
altitudes of these two hypothetical 
paths is a corridor width more easily 
visualized as an aiming slot for the 
returning space vehicle. Entry corridor 
expressed as an allowable range of 
entry angles may be converted to a 


corridor width by calculating the 
appropriate Keplerian perigees and 
observing the difference in perigee 
altitude. 


An entry corridor is required to 
accommodate atmosphere irregularities 
and inaccuracies in approach guidance 
and control. A wide entry corridor 
may be desired to reduce the require- 
ment for flight path deflection by im- 
pulse prior to entry. Corridor width is 
a major problem for vehicles with small 
maneuverability. However, when entry 
angle corridors of several degrees are 
available, there is freedom to make 
design trades between acceleration 
limits, heat load, and entry corridor. 


Glide Maneuvers 


Range and offset control are confined 
to the heating rate region where radi- 
ation cooling allows freedom for ex- 
tended-time maneuvers. This restric- 
tion governs the maximum speed for 
beginning glide maneuvers. In_ this 
study, glide begins at 27,000 fps—about 
the highest speed allowed by the 
temperature placard for maximum L/D 
glide. The corresponding altitude is 
250,000 ft. Equilibrium glide, which 
is essentially horizontal flight with no 
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STRAIGHT WALL 
TANTALUM 
CAPACITOR 
CAN'T LEAK 


Meets MIL C 3965-B, Style CL-64, CL-65. 


A new space-saving approach to the 
design of wet tantalum capacitors 
ends mounting problems encountered 
with flanged types and yet will not leak. 
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ITT’s compact, sintered slug tantaium 
capacitor features a wedge-shaped 
seal held under compression by an 
epoxy retainer ring formulated for 
thermal characteristics inverse to 
those of silver. Ordinary, straight- 
wall capacitors leak along the lead 
when elastomer compression is re- 
duced as the silver can expands. Not 
so with the new ITT design! 


This new, compact capacitor conforms 
to specifications MIL C 3965-B, Style 
CL-64, CL-65 and provides both the 
compactness and rugged reliability 
required in missile, airborne and 
mobile equipment. For details, write 
today requesting Bulletin No. 610. 
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Fig. 4 (left). 


vertical acceleration, is used from the 
end of entry maneuvers to arrival over 
the recovery site. During this phase, 
continuous response to range-to-go guid- 
ance is possible because aerodynamic 
control is maintained by always gliding 
within the atmosphere. The goal in 
the glide maneuver phase is to obtain 
range control which will allow landing 
anywhere in one circumnavigation while 
at the same time maintaining large 
lateral offset capability for all downrange 


distances within the range control 
capability. The result will be a belt of 


landing site coverage which is twice as 
wide as the offset available and which 
completely encircles the earth. 


Trans-Satellite Glide 


Control of range is the primary 
objective for glide from 27,000 to 25,000 
fps. The very low dynamic pressures 
available in this region permit circum- 
navigating glide at maximum L/D with 
a transition through zero-lift attitude 
at satellite speed. Glide at a high-drag 
attitude along a temperature-placard 
trajectory for maximum allowable dy- 
namic pressure allows very short-range 
flights. The small speed change, long- 
range requirement, and high flight speed 
in this phase make it impossible to 
deflect the flight path laterally more 
than a few degrees; consequently, 
this small added turn capability is 
neglected in favor of a simpler concept. 


Subsatellite Glide 


In the subsatellite glide phase begin- 
ning at 25,000 fps, substantial lateral 
offset can be achieved. Near-maximum 
lateral offset is obtained by a simple 
two-step maneuver consisting of glide 
at 45-deg bank angle until about an 
S0-deg heading change is obtained, 
then an unbanked glide is made to the 
maximum offset point. It is known 
that slight improvement can be obtained 
by following a more complex schedule 
of gradually diminishing bank angle. 

The downrange distance at which 
maximum lateral offset is obtained is 
somewhat larger than half of the maxi- 
mum no-turn range. Offset is the 
primary objective in subsatellite glide. 
Therefore, to avoid deviation from 
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Acceleration magnitude and duration (glider with max L/D = 2; entry angle = 9 deg.). 
max L/D = 2). 


Fig. 6 (right). 


maximum offset capability, range con- 
trol in this phase is limited to small 
corrections. 


Glider Characteristics 


Desired glider characteristics implied 
by the preceding discussion may be 
summarized as (1) high lift and high 
drag with the options of high L/D and 
low resultant aerodynamic force, (2) 
trim and control capability allowing the 
use of these options and all available 
conditions in-between, and (3) capability 
of withstanding the heating environ- 
ments. 

This general study is not based on a 
particular glider design but representa- 
tive glider characteristics have been 
assumed. The assumptions are believed 
to be realistic in relation to currently 
recognized design possibilities. 


Lift and Drag 


Example glider lift-drag polars used 
for this study appear in Fig. 2. These 
are known as Newtonian polars because 
the lift and drag-due-to-lift are com- 
ponents of flat-plate normal force given 
by a modified Newtonian relation. 
(See equation in Appendix.) The mod- 
ification, consisting of adjustment of the 
classic Newtonian coefficient from 2.0 
to 1.82, gives more realistic values at 
high angles of attack. Minimum drag 
coefficient accounts for skin friction 
and other parasitic drag items, and it is 
selected to give the appropriate maxi- 
mum L/D. Aside from being con- 
venient for preliminary studies, lift- 
drag polars of this type furnish good 
approximations of hypersonic test re- 
sults above the angle of attack for 
maximum L/D. 

A maximum lift-coefficient of 0.7 
with maximum L/D from 1.5 to 2.0 
appears, on the basis of applicable test 
data and design experience, to be 
obtainable in the trans-satellite glide 
regime. Present indications are that 
higher L/D will not be easily obtained 
with gliders of reasonable size. Aero- 
dynamic uncertainty is very great in the 
entry pull-out regime. But even if 
lower maximum L/D occurs in the 
entry phase, entry can still be performed 
satisfactorily. Maximum L/D of 2.0 
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Fig. 5 (center). Entry heat load (glider with 


Entry design trades (glider with max L/D = 2). 


or more is expected in the subsatellite 
glide phase. 


Trim and Control 


The lift-drag polars in Fig. 2 ter. 
minate just beyond maximum Tift 
Angles of attack beyond maximum-lift 
attitude are not considered in this study, 
Higher angles of attack give higher drag 
coefficients but lower lift coefficients, 
Some entry analyses’ have indicated 
it is better to use angles of attack 
up to 90 deg for the principal purpose oi 
obtaining low heat load. An opposing 
consideration is that trim and contro 
difficulty increases at high angle of 
attack. At a 90-deg angle of attack 
an entirely different sort of control 
system and control surface arrangement 
would be required. Control dynamics 
and response characteristics deteriorate 
badly at extreme angles of attack. 
It is concluded that the heat load 
reduction obtained by flight above 
maximum-lift attitude does not justily 
the effort required to obtain satisfactory 
flying qualities. Furthermore, it & 
doubtful whether heat load reduction 
can be obtained without undesirable 
acceleration load or duration penalties 


Wing Loading 
} 


Glider weight will generally be 
lowest if the wing loading is the max 
mum allowed by material temperature 
limits. This study is based on assumed 
temperature limits of (1) 4,000°F at 
the 6-in. radius stagnation point and 
(2) 2,500°F on the lower surface i 
turbulent flow at an effective distance 0! 
10 ft downstream of the leading edge 
These limits define temperature plac 
ards which form the boundary betweel 
the heat load region and the heat ratt 
region in rig. 1. Aside from the entry 
maneuver, which necessarily plunges 
into the heat load region, the tem 
perature placards are most 
approached during subsatellite banked 
glide. The critical points are at 22,00 
fps for the stagnation point and 19,00 
for the lower surface. At these desigt 
points, the wing loading is adjusted § 
that the glide altitude for maximull 
L/D attitude and 45-deg bank wi 
not be below the 
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One engineer and an EAI PACE TR-10 transistorized desktop analog computer can be the equal of several men 


limited to conventional design tools. The TR-10 multiplies the design capabilities of the most able engineer. One 
TR-10, equipped with a removable patch panel, will serve the needs of a roomful of engineers. This new patch panel 
permits pre-patching of programs away from the computer. It plugs into the computer instantly — valuable 
computation time is conserved. ™@ Versatility of the TR-10 is unequalled. Unique simplicity makes it excellent for 
students or engineers unfamiliar with analog techniques. Accurate to a tenth of one percent, the TR-10 performs 
admirably for the experienced user in 95% of routine engineering problems. And with High-Speed Repetitive 
Operation added, or with several units slaved to operate as one, the TR-10 meets the demands of many advanced 


applications. ™ For technical data on this almost infinitely versatile computer, write for Bulletin No. AC 934. 


FAI ELECTRONIC ASSOCIATES, INC. Long Branch, New Jersey 


Leader in Analogics Analog/Digital Computers Data Reduction Process Control Instruments Computation Service 
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condition restricts the wing loading 
to a maximum of 42 psf for glide at 
maximum L/D of 2.0 and 45-deg bank. 
If maximum L/D is lower, it will occur 
at a higher lift coefficient, and glide 
altitude will be higher unless wing load- 
ing is raised. For simplicity, wing 
loading is kept constant, and gliders 
with maximum L/D below 2 will operate 
with peak temperatures slightly below 
the limits. 


Heat Protection 


Survival during flight in the heat 
load region dictates the use of some type 
of protective system. No particular 
system will be specified in this paper, 
but, for illustrative purposes, an abla- 
tion system for entry maneuvers in 
conjunction with radiation cooling for 
glide maneuvers will be assumed. 

Certain properties of the heat load 
have been observed which affect selec- 
tion of a protection system. The heat 
load occurs as a pulse with duration as 
low as 90 sec. Driving enthalpies are 
large because of the high flight speeds. 
Hot gas radiation input has been 
ignored in this study; although this 
heat load component probably is not 
negligible in all cases considered, it is 
certainly not dominant. Some consider- 
ation of this factor will probably be 
required in the design of a heat-load 
protection system. Because the rela- 
tively thin ablation coatings over much 
of the glider will have limited insulation 
effectiveness, surface temperatures will 
probably have to be below the material 
limit temperatures—i.e., ablation tem- 
perature less than placard temperature. 
Convective heat loads presented in 
this study have been calculated for 
surface temperatures equal to the as- 
sumed material limit temperatures. 


Entry Performance 


Entry performance was calculated 
for the various entry maneuvers de- 
scribed in the Flight Plan section. 
Of primary interest is the heat load 
accumulation and acceleration loading 
prior to trans-satellite glide. 

The performance of two hypersonic 
gliders, with maximum L/D of 1.0 and 
2.0 and lift-drag characteristics described 
previously, was investigated. The wing 
loading of 42 psf was selected on the 
basis of critical temperatures on the 
lower surface as explained previously. 


Basis for Calculations 


Entry performance calculations for 
the initial entry pull-out and portions of 
the pull-up were programed for solution 
on the IBM 7090 high-speed digital 
computer. The complete equations of 
motion for a point mass with aerody- 
namic lift and drag were used.“ A 
spherical nonrotating earth and a 1959 
ARDC model atmosphere were assumed. 
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The initial velocity and altitude were 
35,000 fps and 400,000 ft, respectively. 

The basic trajectory program has 
provisions for holding the resultant 
acceleration constant at a_ selected 
limit value by automatically reducing 
the angle of attack as the dynamic 
pressure increases above its correspond- 
ing limit value. This procedure main- 
tains lift as high as possible, and the 
prescribed resultant acceleration is 
held constant. 

The heating rate and heat load 
accumulation per unit area were cal- 
culated in a subroutine to the basic 
trajectory program. A modified Fay- 
Riddell heating rate equation was used 
for the nose stagnation point calculation. 
The reference temperature method for 
turbulent flow with adjustments for 
angle of attack variation was used for 
the lower surface reference point. 
These heating calculation procedures 
are given in reference 11. The radius of 
curvature at the nose stagnation point 
is selected as 1/2 ft, and the reference 
point for turbulent heating on the lower 
surface is 10 ft back from the leading 
edge. This choice of lower surface 
reference point gives, within 10 percent, 
the average heating rate for the entire 
undersurface of flat bottom delta wings 
with root chords from 35 to 90 ft. 
These heating data may be scaled to 
other reference lengths by applying 
the well-known rules that stagnation 
point heating varies inversely as the 
square root of the nose radius, and that 
lower surface turbulent heating varies 
inversely as the fifth root of the distance 
from the leading edge 


Undershoot-Overshoot Corridors 


The results of the investigation 
of initial entry pull-out for 35,000 fps 
entry are shown in Fig. 3. The maxi- 
mum permissible entry angle for a 
prescribed limit in resultant acceleration 
for the two gliders is shown by the 
upper two curves. At the maximum 
permissible entry angle, entry is made 
at maximum lift and attitude is modu- 
lated to nearly zero lift and minimum 
drag at peak dynamic pressure. These 
data confirm the conclusions of Grant.! 

Results are also shown for the L/D 
1.0 glider entering at fixed attitude. 
The attitude is that for maximum L/D 
(W/SC;z is 108 psf) which gives approxi- 
mately minimum resultant acceleration 


for a fixed attitude entry. The mini- 
mum entry angle (overshoot limit) 
for these gliders is also shown. Both 


gliders have the same maximum lift 
coefficient and wing loading; therefore, 
both have the same overshoot limit. 
Data for the ballistic vehicle, a 
peculiar feature of which requires expla- 
nation, is also shown in Fig. 3. When 
acceleration level is below about 8g, 
minimum entry angle is governed by 
acceleration rather than overshoot. 
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Hildebrand? and Styer!? have described 
a second peak acceleration which occurs 
below satellite speed for ballistic entries 
near the overshoot limit. It is this 
second peak acceleration which limits 
minimum entry angle in this case 

The entry corridors for the 10g limit 
converted to virtual perigee differences 
as suggested by Chapman are 95 
nautical miles for the ballistic vehicle. 
65 for the fixed attitude L/D 19 
glider, and 185 for the attitude modu. 
lated L/D 2.0 glider. These corridors 
narrow slightly as entry velocity in- 
creases, and, as Chapman® has shown, 
all 10g ballistic corridors vanish at 
46,300 fps. 

The ability of an attitude modulated 
glider to limit the g level during steep 
entries depends upon the average lift 
throughout the pull-out. Average lift 
determines the depth of atmosphere 
penetration and, therefore, the peak 
dynamic pressure. If drag is low, 
lift can be larger for a given resultant 
acceleration; therefore low drag, high 
L/D gliders can control g level more 
effectively than aerodynamically in- 
ferior vehicles. 

It is doubtful that the very large entry 
angle corridors shown will be required 
for manned entries; therefore, a trade 
for lower g’s will be made. 


Acceleration Magnitude and Duration 


The attitude modulation mode of 
entry allows a limit to be placed on 
the magnitude of the g load, but this 
load must be sustained for a period of 
time. The added duration leads to 
concern that the g-time combinations 
might exceed human tolerances. Fig. 4 
shows representative time histories of 
the resultant acceleration for entries at 
maximum lift unmodulated and modu- 
lated to lower g values. Also presented 
are human g-time tolerances for normal- 
to-spine and parallel-to-spine loadings 
as determined from human centrifuge 
tests and reported by Creer etal.? These 
data show that reduction of g level by 
attitude modulation improves the com- 
fort level despite the extended duration 
of the reduced g level. The comfort 
level can be seen to increase rapidly as 
the resultant acceleration level 1s 
lowered. Comfort level is indicated by 
the amount the allowed duration exceeds 
the actual duration at a given accelera- 
tion level. 

For the entries shown, the normal-to- 
spine tolerance is not exceeded even for 
the unmodulated entry. However, the 
parallel-to-spine tolerance is exceeded 
at greater than 8g resultant acceleration 
When the pilot is seated upright im 
the glider, resultant acceleration is 
essentially parallel-to-spine between 
and 55-deg angle of attack; thus the 
lower human tolerance curve in Fig. 4 
will apply. 


Heat 
It 
lowe 
Fig 
accu 
out 
glide 
heat 
entr 
tole 
und 
The 
for t 
Ir 
per 
but 
invo 
per 
it in 
of tk 
is as 
Whe 
fron 
to n 
pull 
r 
con 
lift ¢ 
obte 
The 
stee 
heat 
sho 
ave 
moc 
leve 
entr 
volv 
atm 
1S 
ave 
for 
The 
the 
rate 
of a 
A 
ior 
this 
glid 
glid 
esse 
hea 
diff 
ma) 
perc 
mo 
duc 
is 
of 
Ent 
T 
and 
reg: 
low 
IS 
wei; 
wei 
uni 


bed Heat Loads 


- The stagnation point and turbulent 
es lower surface heat loads presented on 
is Fig. 5 are the summation of heat loads 
nits F cumulated during both the initial pull- 
: out and the subsequent pull-up to the 
mit glide initiation point. The figure shows 
Qe heat load per square foot plotted against 
try angle tolerance. Entry angle 
the Increased technical responsibilities in the field 
and overshoot entry Emits. of range measurements have required the cre- 
The overshoot limit is fixed at 4.4 deg f i 3 
lors 9 for this glider. ation of new positions at the Lincoln Labora- 
tory. We invite inquiries from senior members 
* Ff but only a small area of the glider is of the scientific community interested in par- 
involved. The lower surface heat load th “ bl th 
per unit area is more significant because ticipating wr US SOLVING proolems of 
it involves the entire wing lower surface. n 1 N 
a. greatest urgency in the defense of the nation. 
ts of the total heat for shallow entries and 
ere 


oak isas little as 10 percent for steep entries. RADIO PHYSICS 


When comparing these data with data 


low, from other general studies, it isimportant and ASTRONOMY 


to note whether heat load beyond the 
-out has been included. 

RE-ENTRY PHYSICS 
The maximum lift lines shown denote 

constant attitude entries at maximum 


lift and indicate the minimum heat load PEN ET RATIO N Al DS 
cred obtainable for a given entry angle. DEVELOPMENT 


The heat load decreases as the entry is 
steepened because the higher average 
heating rate is more than offset by the TARG ET IDENTIFICATION 


shorter entry time caused by the higher 


on average deceleration. When attitude RESEARCH 


modulation is employed to limit the g 


jon | (evel, the heat load increases as the SYSTEMS: 
entry is steepened. Modulation in- Space Surveillance 
in volves deeper penetration into the Strategic Communications 


atmosphere because the lift coefficient 
is reduced. The resulting increased 


7 average heating rate is the main reason NEW RADAR TECHN IQU ES 
<a for the higher heat load in this case. 
ona The lower surface heat load also reflects 


ot the secondary variation in heating SYSTEM ANALYSIS 


al rate resulting from changes in angle 
ak of attack during modulation. COMMUNICATIONS: 
dings Although the data presented are 
ifuge for a glider with maximum L/D of 2.0, . y gy y 
hese this discussion applies also to lower L/D 


ol gliders. Analysis of heat loads for a INFORMATION PROCESSING 


cole glider with an L/D of 1.0 showed 


ation | Ssentially the same results except that SOLID STATE 


heat loads were somewhat lower. The 
nfort 
ly as difference is negligible for fixed attitude Physics, Chemistry, and Metallurgy 
13 maximum lift entries, but up to 15 
ceeds modulated attitude entries. This re- 


slera | ‘duction occurs because exposure time work will be sent to you upon request. 
is shortened by the higher drag level 


of lower L/D gliders. 


al-to- All qualified applicants will receive consideration for employment without 
on for Entry Design Trades regard to race, creed, color or national origin. 
r, the 
eeded The heat load, entry angle corridor, 
ation and can Research and Development 
in regardec 
iS ome: nak LINCOLN LABORATORY 
on 15- is indicative of heat protection system 


Massachusetts Institute of Technology 
BOX 23 
LEXINGTON 73, MASSACHUSETTS 


is the | Weight. An estimate of ablation system 
Fig. 4 weight is obtained when heat load per 
unit area is divided by the wing loading 


October 1961 + Aerospace Engineering 39 


T 
| 

| 
| 
F 


000 400 6000 


GREAT-CIRCLE DISTANCE (Nautical Miles) 


Fig. 7 (left). Subsatellite glide offset. 


and effective heat of ablation. The 
result is the fraction of gross weight 
to be used for lower surface ablation ma- 
terial. The entry angle tolerance is se- 
lected on the basis of pre-entry guidance 
accuracy and approach path control 
considerations. Resultant acceleration 
is the crew comfort index. The de- 
signer’s choice among these factors is 
shown in Fig. 6. 

A representative glider design might 
have a design heat load limit for the 
lower surface of 15,000 Btu/sq ft. 
With 42 psf wing loading, an ablation 
weight of 6 percent of the total vehicle 
weight is required if an effective heat of 
ablation of 6,000 Btu/lb is assumed. 
This design heating limit provides for a 
maximum entry angle tolerance of 6 
deg for maximum acceleration of 10g. 
With better control of entry angle so 
that the tolerance can be reduced to 3 
deg (40 nm corridor), the resultant 
acceleration will not exceed 4g. For 
entries within the entry angle tolerance, 
lower heat loads could result, but it is 
preferable to lower the acceleration 
level. 


Glide Performance 


An evaluation of glide performance 
will reveal when and under what con- 
ditions it is possible to glide to a pre- 
selected landing site after entry has 
been achieved. Glide maneuvers be- 
ginning at 27,000 fps allow global range 
control and substantial lateral offset. 
These capabilities are expressed as a 
belt of landing site coverage which 
encircles the earth. The half-width 
of this belt is equal to the maximum 
lateral offset obtained during sub- 
satellite glide. The freedom to land 
anywhere along this belt is obtained by 
control of range during trans-satellite 


glide. This evaluation, rather than 
determining specific flight paths to 


certain sites, will establish the limits 
of capability to reach preselected land- 
ing sites. 
Trans-Satellite Glide 

A wide variety of operating procedures 


is available for use in the trans-satellite 
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Fig. 8 (right). 


SLIDE RANC 


Glide range (from 27,000 fps). 


speed regime (taken here to be from 
27,000 to 25,000 fps Because range 
control, rather than range per se, is 
the objective in trans-satellite glide, 
procedures which allow the vehicle to 
leave the atmosphere are avoided. By 
gliding within the atmosphere, sufficient 
range control is obtained and continuous 
response to range-to-go guidance is 
possible. In general, control of range 
is obtained by regulating the rate 
at which energy is dissipated. In this 
case, range control involves an inter- 
dependent control of flight path and 
glider attitude. Low dynamic pressure 
flight paths and low-drag attitudes give 
long range. Conversely, high dynamic 
pressure and high-drag attitude give 
short range. The interdependence be- 
tween flight path and glider attitude 
arises from the requirement for vertical 
force equilibrium. The flight path 
must be essentially horizontal to re- 
main in the atmosphere therefore, 
lift is required to maintain vertical 
forces in balance. 


Long-Range Glide 


Maximum irans-satellite glide range 
is obtained in a three-step maneuver 
following entry termination at 27,000 


fps and 250,000 ft altitude. The 
altitude-velocity profile of the three 
steps is included in Fig The first 
step is flown inverted at maximum 


L/D attitude along the equilibrium- 
glide speed-altitude schedule for the 
corresponding lift coefficient and 42 psf 
wing loading. This step is terminated 
just above 300,000 ft altitude where 
the dynamic pressure has decreased to 
1 psf and speed is about 100 fps above 
satellite speed. Through step 2 the 
dynamic pressure is held constant. 
The transition from constant lift coeffi- 
cient glide to constant dynamic pressure 
glide requires a slight change in flight 
path angle. It is assumed that this 
transition maneuver will cause no range 
penalty. Vertical force equilibrium is 
maintained during step 2 by gradua'ly 
reducing angle of attack until zero-lift 
attitude is reached at satellite speed. 
Then the glider is rolled upright and 


October 1961 


angle of attack is gradually increased g 
speed reduces. Maximum L/D at 
titude is regained about 100 fps beloy 
satellite speed, and the third step, g 
constant attitude glide to 25,000 fps 
begins. The ranges in these three 
steps for a glider with maximum L/) 
of 2 are 11,600, 11,620, and 9,920 nm 


respectively. (See equations in Ap 
pendix.) The total glide range j 


33,140 nm for this trans-satellite maney 
ver. Range from the entry point adds 
an additional 450-2,300 nm. Maxi 
mum trans-satellite glide range for L/] 
from | to 2 is included in Fig. 8. 


Short- Range Glide 


Very short range results from opera- 
tion at the highest allowable drag level 
A dynamic pressure of 30 psf is availabk 
just above the temperature placard 
altitude, and a drag coefficient of about 
1.0 occurs at maximum lift attitude 
With wing loading of 42 psf, this com 
bination yields a deceleration of about 
0.7g. Speed will dissipate from 27,000 
to 25,000 fps in about 90 sec; the dis 
tance flown will be under 400 nm. Dur- 
ing this maneuver the vertical lift 
component is regulated to maintain 
the desired trajectory by varying 
bank angle from a nominal value « 
90 deg. Because of the short duratio 
this maneuver may require excessivel 
rapid control response. To allow more 
time, the minimum trans-satellite glide 
range is assumed to be 900 rather than 
100 nm. 


Trans-Satellite Turn Performance 


Lateral offset maneuvers are ex 
cluded from the trans-satellite glide 
phase as indicated in the Flight Plan dis- 
cussion. Analysis of the trans-satellite 
turn performance with L/D of 2.0 con 
firmed that the operational complexity 
introduced is not yet justified by the 
available performance gain. To obtain 
an offset at a large downrange distance 
by continuous turning, it is necessary 
to alternate bank. The bank angle 
schedule required is different for each 
downrange distance. The theoretical 
maximum contribution to lateral offset 
using 90-deg bank turns between 
27,000 and 25,000 fps is 465 nm. This 
is reduced to 306 nm when bank angles 
are limited to those allowed by equilib- 
rium requirements and the temperature 
placard. Both of these values are ac 
companied by less than 4,000 nm range 
The lateral offset contribution at 4 


point one circumnavigation farther 
downrange is 114 nm. Because of the 


once-around range requirement, only 
the smaller contribution can be credited 
to widening the coverage belt. These 
contributions are not negligible, but 
they are small when compared with 
subsatellite lateral offset contributions 
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Sybsatellite Glide 


Fig. 7 shows banked glide performance 
for maximum lateral offset with turns 
beginning at 25,000 fps. This maxi- 
mum lateral offset will establish the 
width of the landing site coverage belt; 
the downrange distance at which 
maximum offset occurs is to be included 
in the glide range. 
lateral offset for glide 
25,000 fps in a _ nonrotating 
atmosphere is approximately propor- 
tional to (L/D)1:7 for L/D from 0.5 to 
20, This variation has been suggested 
fore? and has been verified by more 
exact analysis.’ Downrange distance to 
mé lateral offset varies from 
shout 60 to 70 percent of maximum 
no-turn range as L/D varies from 2.0 
to 1.0. The data in Fig. 7 are based on 
approximate relations and a 
nominal maximum offset of 2,100 nm 
at an L/D of 2.0. 


Maximum 
from 


th 


The maximum no-turn glide range 
fom 25,000 fps is included in the 
maximum range shown in Fig. 8. 


This range is given by the constant- 
ttitude equilibrium-glide range equa- 
tion presented in the Appendix.) The 
minimum no-turn range from 25,000 
fps is assumed to be 2,400 nm for any 
L/D between 1.0 and 2.0. Even shorter 
ange can be obtained by high angle of 
ittack banked maneuvers similar to 
those suggested for minimum trans- 
satellite glide range. Note that Fig. 8 
shows only glide range; range covered 
luring the initial entry maneuver (450 
to 2,300 nm) is not included. 


Coverage Belts 


Sufficient glide range data have been 
assembled to establish that a coverage 
belt extending around the earth is pos- 
sible. The belt begins at the minimum 
range at which maximum offset can 
be obtained, and ends at the maximum 
range at which maximum offset can be 
obtained. When the difference between 
these limits equals or exceeds the 
circumference of the earth, the desired 
belt exists. As indicated in Fig. 8, 
this is true when maximum L/D is 
greater than 1.52. 

The half-width of the belt is the 
maximum lateral offset given in Fig. 7. 
This varies from 1,325 to 2,100 nm for 
maximun L/D from 1.52 to 2.0. 
Landing site coverage belts for these 
two cases are shown over a Mercator 
map in Fig. 9. On a globe, the belt 
centerline or baseline is a great circle and 
the edges are parallel minor circles. 
The belts are shown at an arbitrarily 
chosen inclination of 45 deg to the 
equator. The map shows more than 
360° longitude in order that the over- 
lapping ends of the inner belt (L/D = 
1.52) could be shown stretched out full 
length. The wider belt (L/D = 2.0) 


is longer and when stretched out ex- 
tends outside the map borders. 

The wider landing site coverage belt 
encloses the preselected landing site 
indicated in Fig. 9 for which Omaha 
has been chosen. Thus, landing there 
is possible by turning left with a glider 
having maximum L/D of 2.0. The 
belts are fixed in space, and earth 
rotation will shift Omaha to the right 
as time elapses. About one hour later, 
landing at Omaha will be possible for 
a glider with maximum L/D of 1.52. 
About 13 hours later, Omaha will be 
outside both belts on the right. This 
example illustrates the effect of time 
on accessibility of a preselected landing 
site. The baseline inclination and land- 
ing site were arbitrarily selected. Choice 
of another inclination and/or landing 
site latitude would alter the situation. 

Because the belts circle the earth, the 
entry point could be located anywhere 
along the space-fixed baseline without 
changing coverages. A new entry point 
location would require a different glide 
program to reach a given site within the 
belt, but the definition of the belt 
ensures that the site can be reached. 


Space Mission Planning and 
Conditions at Entry 


Round-trip space missions must be 
scheduled to terminate with entry 
conditions which are suitable for the 
vehicle used. 

The approach trajectory leading to the 
atmosphere entry point can be specified 
by the following entry point conditions: 
(1) entry speed, (2) entry altitude, 
(3) entry angle, (4) space-referenced 
(or inertial) entry longitude, (5) ap- 
proach plane inclination to the equator, 
and (6) time at entry. 

Desired entry conditions may be 
specified in advance and the entire 
space mission scheduled to suit the 
resulting restrictions. Because inevit- 
able space trajectory dispersions and 
perturbations will cause variations from 
intended entry conditions, provisions 
will be made for enroute path corrections. 
Some space mission objectives may be 
impossible without major impulsive 
maneuvers to obtain specified entry 
conditions. Thus, space mission sched- 
uling is either hampered or penalized 
in weight by advance specification of 
entry conditions, but these difficulties 
cannot be avoided where entry vehicles 
with small maneuverability are used. 
The relaxation of entry condition restric- 
tions by the use of more sophisticated 
entry vehicles will allow greater freedom 
in space mission scheduling 

Consider, for example, moon-earth 
trajectories when entry angle is the 
only entry condition specified. Entry 
speed is essentially constant for practical 
trajectories. Time of arrival at the 
entry point will depend almost directly 
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work in the fields of the future at NAA 


DYNAMICS 
ENGINEERS 


Immediate openings for Struc- 
tural Dynamics, Vibration, and 
Flutter Engineers to work on 
advanced present and future 
manned aerodynamic vehicles. 
Graduate Engineers in Aero- 
nautical, Mechanical or Ap- 
plied Mechanics needed for: 
Vibration Prediction 
Structural Vibration Analysis 
Dynamic Response 
Calculation 
Oscillatory Aerodynamics 
Specialists or Senior Engi- 
neers have excellent opportu- 
nities in the fields of vibration, 
dynamic analysis, structural 
response, and flutter. Openings 
available at all levels of expe- 
rience and responsibility. 
Experience in at least one of 
these fields, plus an advanced 
degree in AE, ME, or Applied 
Mechanics required. 


Write: Mr. H.K. Bowman, Engi- 
neering Personnel, North 
American Aviation, Inc., Los 
Angeles 45, California. 

All qualified applicants will receive consid- 
eration for employment without regard to 
race, creed, color, or national origin. 
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on time of departure from the moon for on the two lower curves at 45-deg with L/D of 2.0 on an equatorial ap- less 
a given type of trajectory—i.e., energy approach plane inclination with 11- and proach plane. Another feature of Fig. crev 
level, direct or retrograde, etc. Entry 14-hour coverage. The curves show 10 should be noted. For high inclina- real 
point inertial longitude will vary with the variation of coverage for other ap- tions, to the right of the break in the assu 
calendar time of departure at about proach plane inclinations. A_substi- three lower curves, the coverage occurs bea 
13° per day. It will also vary with tute for lateral offset is enlargement in two equal periods. In this case, syst 
entry angle and type of trajectory. of the site to a circular region and ac- the coverage belts go over the poles glid 
Approach plane inclination is the most ceptance of a landing anywhere in and are intercepted twice a day by the site 
uncertain item because it varies greatly the circle. The increase in coverage landing site as it rotates with the earth. pat 
with minor variations in trajectory near for circular recovery areas centered at Many possibilities are suggested by retu 
the moon. It appears that relaxation 40° latitude is shown by the upper pair these limited data but the point is wit] 
of restrictions on approach plane in- of curves in Fig. 10. A_ 3,300-nm that a high L/D glider can land at some T 
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month. Fig. 10 shows that equatorial entry allows great freedom in_ scheduling exp 
The time of day determines the (zero approach plane inclination) with the space mission. witl 
geographic longitude of the entry point point landing sites is excluded for BAe. veh 
for a given inertial longitude, and time —_ both gliders. If the site were moved Other Possibilities | and 
of the month determines inertial longi- south 300 nm (to 35° latitude), 12-hour The utilitarian advantages of highly } L/l 
tude of the entry point. This relation- _ coverage would be available to the maneuverable gliders for entry and i 
ship and its consequences were implied higher L/D glider. However, there recovery are believed to be unique. spa 
in the discussion of coverage belts on the will be a slight sacrifice in coverage at It may, however, be possible to obtain refi 
Mercator map. high inclinations. Similar results for comparable entry and recovery per- pen 
: the glider with L/D of 1.52 would formance with aerodynamically less 
Landing Site Coverage require a landing site at 22° latitude. sophisticated vehicles by including Ap 
In the discussion of space-fixed Another possibility exists for allowing impulsive maneuvers of various kinds. 
coverage belts on the Mercator map equatorial return. The approach plane One especially interesting possibility is Ney 
(Fig. 9), it was shown that landing inclination can be adjusted, in effect, injection into a near-earth waiting orbit 7 
site coverage is independent of entry by heading changes achieved before the after a direct entry has reduced speed Ne 
point location along the _ baseline. glide phase. For steeper entry angles to near-satellite speed. Then de-orbit . 
This means that inertial longitude of and modulation to low acceleration could be arranged to allow widespread 
the entry point does not affect ability to limits, there is opportunity for signifi- choice of landing sites. In this type of 
glide to a preselected site. It was also cant heading change during the pull-up. recovery operation with a less maneu- 
observed that landing at the selected Banked flight is used to control the verable vehicle, entry angle would 
site was possible during a period of flight path angle, and heading changes necessarily be more closely controlled whi 
several hours. Thus, entry could take up to 7 or 8 deg occur incidentally. not only to get comparable comfort w 
place any time during a period of the Shallow entries yield smaller heading levels for the crew, but also to ensure onl 
same duration. It was noted that changes which approach zero for over- capability of achieving the desired ‘ 
results would change for a different shoot entries. Note that these heading waiting orbit. 2.0 
landing site latitude and/or baseline changes are more useful than those Many other schemes will be advanced stu 
inclination. (Baseline and approach which might be obtained during trans- for refining the comfort and utility of ter: 
plane inclinations are nominally the satellite glide because they occur before, entry vehicles by the time their ap- toll 
same.) Fig. 10 shows how the com- rather than during, long-range glide. plication is required. In choosing ton 


bined effects of these factors determine 
the prospects for landing at a 40° 
latitude site. Landing site coverage is 
expressed as the number of hours per 
day that the landing site can be reached 
from the entry point. 

The two cases shown on the map 
(Fig. 9) are included in Fig. 10 as points 
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For an adjustment of approach plane 
inclination due to heading change 
during pull-up, the horizontal scale of 
Fig. 10 is displaced to the right or 
left an amount equal to the heading 
change. A 7-deg heading change ad- 
justment, when available, will give 
about 6-hour coverage for the glider 


October 1961 


among the various possibilities, low 
system weight will always be a major 
factor of merit. Because of the current 
booster situation, low system weight 
is the dominant factor now amounting 
to a GO-NO-GO feasibility criterion. 
But in the future the opportunity 
will be available to trade weight for 


sin 


najor 
rrent 
eight 
nting 
rion. 

unity 
for 


utility. Then there will be a demand 
for entry vehicles with greater maneuver 
capability. This capability will allow 
entry with low acceleration loads on 
the ‘crew and recovery with greater 
flexibility in mission planning. 


Conclusions 


The use of a hypersonic glider for 
direct entry into the atmosphere at 
escape speed and subsequent controlled 
fight to a preselected landing site has 
been considered and found promising. 

Entry corridors on the order of 60 
miles wide can be maintained with 
acceleration levels controlled to 6g or 
less by attitude modulation. Thus, 
crew comfort levels consistent with 
really reliable crew performance are 
assured. This entry performance can 
be achieved with feasible heat protection 
system weights. After entry, controlled 
glide to a previously selected landing 
site is possible regardless of the approach 
path location or orientation, while 
return from space can occur any time 
within a period of several hours. 

These are the predominant perform- 
ance advantages of a glider used for 
entry and recovery; they are the 
product of the basic maneuverability 
of a glider. The remaining question is 
whether these advantages justify the 
weight penalty that may reasonably be 
expected when a glider is compared 
with less aerodynamically sophisticated 
vehicles. This concept of refined entry 
and recovery maneuvers using a high 
L/D glider is submitted with the ex- 
pectation that our rapidly expanding 
space technology will soon require such 
refinements and a reasonable weight 
penalty will be justified. 


Appendix 


Newtonian Lift-Drag Polars 


The lift and drag coefficients for a 
Newtonian polar are 


C, = Nsin? a cos a 
Cp = Nsin® a + min Cp 


where V = Newtonian coefficient, and 
w = angle of attack (positive values 
only), 

The classic Newtonian coefficient is 
2.0, but it is modified to 1.82 for this 
study as explained under Glider Charac- 
teristics. These additional relations 
follow from the foregoing basic New- 
tonian polar equation. 


sin ax {K[1 —K?+(l — 34 
— K?— (ij — 2K2)"2] —K 


sin? ao COS ao 
max L/D = 
K + sin? ap 


max C, = 0.385N 


where N = modified Newtonian coeffi- 
cient, K = min Cp/N, and a = 
angle of attack at maximum L/D. 

Maximum Cy, occurs at a 54%/, deg 
angle of attack in all cases. 


Range Equations 


All glide operations are performed 
with essentially horizontal flight—i.e., 
the flight path angle is very small 
(less than 1/10 deg for L/D of 2) and 
changes very slowly with time. This al- 
lows simplification of the equations of 
motion leading to this familiar expression 
for calculating glide range in steps 1 and 
3 of trans-satellite glide 


AR = (r/2) (L/D) cos ¢ 1 (a 
= (r/2) (L cos ¢ In 1_ v2 


where AR = range increment, r = radius 
of earth, ¢ = bank angle, and V = ratio 
of velocity to circular satellite velocity. 
Subscripts 1 and 2 refer to initial and 
final conditions. 

This equation works above and below 
satellite speed; but, because constant 
L/D and bank angle have been assumed, 
it is discontinuous at satellite speed. 

Glide across satellite speed requires 
variation of attitude—either bank angle 
or angle of attack or both—to satisfy 
the equilibrium rules implied by the 
assumption of horizontal flight. 

The following general equation works 
through satellite speed because bank 
angle has been varied to satisfy equilib- 
rium requirements. Because angle of 
attack and dynamic pressure are held 
constant, it is a constant-drag glide. 
This equation is not directly applicable 
to the step 2 trans-satellite glide but will 
be modified for use as explained here: 


AR = (r/2) (W/Cp Sq) (Vi? — V2?) 


where W/S = wing loading, Cp = drag 
coefficient, and g = dynamic pressure. 
In step 2 of the trans-satellite glide, 
the initial and terminal points are on 
constant-attitude equilibrium — glide 
speed-altitude schedules. For this spe- 
cial case, it can be shown that the last 
factor in the foregoing equation is given 
approximately by the following: 
— 2) = 
W/S 


where C,* = life coefficient of glide 
preceding 1 and following 2, and ¢* = 
bank angle of glide preceding 1 and fol- 
lowing 2. 

Step 2 of the trans-satellite glide in- 
volves variation of angle of attack at 
zero (or 180 deg) bank angle in order to 
maintain equilibrium. Angle of attack 
varies between zero lift and maximum 
L/D attitudes. The drag coefficient 
varies between the minimum and the 
value at maximum L/D. Within this 
variation there is an effective constant 
value which yields equivalent results. 
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The effective value of drag coefficient 
was found to be approximately 


effective Cp = min Cp + 


1/3 [(Cp at max L/D) — (min Cp)] 


For modified Newtonian polars these 
drag coefficients can be expressed in 
terms of ao when higher order terms are 
omitted. 


Cp at max L/D = 3 (min Cp) cos? ao 
effective Cp = (min Cp) X 
(2 3 + cos? ag) 


Substituting the alternate form for 
the velocity factor and effective Cp for 
Cp in the range equation above and 
recognizing that C,* is now Cy, for 
maximum L/D and ¢* is either zero or 
180 deg, the equation for range in step 
2 is obtained: 


3 cos? ap 


AR = r(max L/D) 


3 + cos? ao 


where ap = angle of attack for max L/D 
on the Newtonian Polar, and r = radius 
of earth. 

Note that this range increment de- 
pends only on lift-drag characteristics 
and is unaffected by the altitude (or 
dynamic pressure) at which step 2 
crosses satellite speed. This is reason- 
able because a lower altitude (higher 
dynamic pressure) crossover for which 
drag would be larger would also have 
a larger speed difference and, therefore, 
more energy to dissipate. However, 
crossover at lower altitude would reduce 
range in steps 1 and 3. 
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ing edge sections are attached to the 
main structural box. There is ob- 
viously no difficulty in providing span- 
wise or chordwise gaps less than 0.10 
in. wide. 

A typical wing structure that satisfies 
LFC smoothness criteria is shown in 
Fig. 5. The slot spacing on a jet 
transport size wing would average ap- 
proximately 2 in., being closer together 
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near the leading and trailing edges and 
farther apart near the maximum thick- 
ness point of the airfoil. It should be 
noted that the leading and trailing edge 
sections serve as full depth LFC ducts. 
In the area of the main structural box 
where suction requirements are lower, 
ducts are part of the skin panel, and 
integral fuel volume is provided similar 
to a conventional turbulent-type wing. 
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Wing assembly procedures are discussed 
in the Manufacturing Techniques sec. 
tion of this paper. Another paper 
presents some structural aspects of g 
typical LFC wing. The type of wing 
structure shown in Fig. 5 is utilized on 
the demonstration aircraft program 


Demonstration Aircraft Program 


In August 1960, the Norair Division 
of Northrop Corporation was awarded 
an Air Force contract to modify two 
B-66 aircraft to a laminar flow control 
demonstration airplane configuration, 
The B-66 was chosen because of its 
ability to simulate in size and speed the 
operating Reynolds number range of 
transport aircraft, while at the same 
time meeting cost and schedule con- 
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Fig. 2. LFC schematic. 


siderations. The basic objectives of 
the program are: 
(1) To demonstrate the practicality 

of LFC 

(a) Design—To show the feasibility 
of integrating LFC requirements with 
the smoothness criteria, fuel capacity, 
access provisions, and fatigue integrity 
expected in a typical, sound wing struc- 
ture. 

(b) Manufacturing—To prove that 
LFC structure can, in fact, be fabricated 
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Fig. 3. Research results. 


successfully at reasonable cost, utilizing 
present state-of-the-art manufacturing 
techniques. 

(c) Operation and maintenance 
To demonstrate the compatibility of 
LFC with various weather conditions 
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Fig. 4. Smoothness requirements. 


and other operational environments 
normally encountered by a service-type 
airplane, and to establish the LFC 
maintenance requirements. 

(2) To provide data for direct applica- 
tion to future aircraft systems. 

The LFC airplane, illustrated in Fig. 
6, offers such significant changes from 
the original airplane as increased wing 
area, aft-mounted engines, and pumping 
nacelles underneath the wing. 

The wing area was increased from 780 
to 1,250 sq ft to establish characteristics 
similar to typical operational LFC 


Fig. 5. Typical LFC wing section. 


airplanes, such as the optimum cruise 
aerodynamic condition, the laminar area 
percentage of total airplane wetted area, 
take-off and landing characteristics, 
representative amount of fuel in the 
wing, and wing structural similarities. 
The new wing will have, for test pur- 
poses, electrical cyclic deicing equipment 
on a 10-ft section of the leading edge. 
The aft-fuselage engine location was 
selected to broaden the range of LFC 
experimentation on the wing itself. 
A wide spectrum of test conditions in- 
cluding many off-design points has been 
planned to establish firmly the limits of 


Fig. 6. LFC airplane. 
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LFC operational effectiveness, and this 
is most efficiently accomplished with 


fuselage-mounted engines. Wing- 
mounted turbofan engines are com- 
patible with LFC using established 


acoustic criteria and sound suppression 
techniques.’ 

Pumping nacelles are located on the 
wing, as they would be on a typical jet 
transport. A cutaway of the LFC 
demonstration aircraft pumping nacelle 
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Fig. 7. LFC pumping system. 


is shown in Fig. 7. The suction com- 
pressors in this case are driven by bleed 
and bleed-burn turbines operating on 
air from the main propulsion engines. 
In a jet transport with wing-mounted 
engines, the pumping system would be 
located within the engine pylon. 

Status of the LFC demonstration air- 
craft program is highlighted at present 
by the ‘‘finalized’”’ design of LFC slots 
hole sizes, and complete duct system. 
Wing tooling fabrication is in progress, 
and ground tests of the complete air- 
plane will start in October 1962, fol- 
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1963. 


Manufacturing Techniques 


A typical LFC wing production break- 
down is shown in Fig. 8. A notable 
manufacturing feature of this wing is the 
extensive use of proved bonding tech- 
niques to satisfy the LFC contour and 
smoothness requirements while assuring 
production costs per pound equal to 
conventional structures 

The upper and lower skin panels and 
the leading edge and trailing edge sec- 
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Fig. 9. Machining contoured fixture. 


tions are bonded assemblies. These 
in turn, are attached to each other and 
to the ribs and spars with mechanical 
fasteners. 

One of the first steps in manufacturing 
the wing is the construction of a number 
of bonding fixtures approximately 4 by 
20 ft in size to serve as pressure bases 
in the bonding operation. These fix. 
tures match the wing surface to an 
accuracy of 0.002 and are wave-free, 
The fixture material is aluminum tooling 
plate, machined with a numerically con- 
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Fig. 10. Skin bonding. 


STRINGER 


trolled 
Fig. 9. 

Skin panels are then formed, ma- 
chined as required, and put through 
several consecutive bonding operations 
as shown in Fig. 10. All bonding is 
accomplished by use of the previously 
mentioned bonding fixture as a pressure 
base. 

In step 1 (Fig. 10), the plenum cham- 
ber grooves and suction holes are ma- 
chined and drilled in the structural skin 
utilizing an accurately indexed tooling 
fixture. In step 2, the structural skin 
and outer cover, still unslotted, are 


profile mill as illustrated in 
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Fig. 11. Manufacturing cost comparison. 
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wnded together. Next, the stringers 
are bonded to the outer skin assembly; 
fnally, the inner skin assembly, com- 
olete with bonded rib-attach angles, is 
bonded to the stringers. 

The complete skin panel, approxi- 
mately 4 by 20 ft in size, is then slotted 
ysing a radial-type saw attachment, ac- 
curately indexed to the same tooling 
fxture used in machining the plenum 
chamber slots and in drilling the suction 
holes. Slotting is accomplished on the 
upstroke movement of the saw blade, 
thus avoiding burrs on the slot edges in- 
side the plenum chamber. 

A manufacturing cost comparison be- 
tween conventional turbulent structure 
and LFC-type structure is shown in 
Fig. 11. The turbulent structure re- 
ferred to here is the conventional skin- 
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Fig. 12. LFC maintenance. 


stringer-spar type of construction which 
actually becomes more expensive per 
pound than the LFC structure when 
large numbers are produced. This is 
primarily the result of fewer parts in 
the LFC wing, with its large bonded 
assemblies and lower tooling require- 
ments 

The information shown in Fig. 11 is 
the result of (1) extensive manufacturing 
research on LFC-type structures, in- 
cluding tool development and opera- 
tional experience, and (2) firsthand ex- 
perience in the fabrication of conven- 
tional jet transport wing structures. 

Conceivably, the conventional wing 
could be redesigned for bonded con- 
struction, but the significant feature of 
Fig. 11 is the fact that the high accuracy 
requirements of the LFC wing have been 
more than compensated for, from the 
point of view of cost, relative to con- 
ventional construction. 


Maintenance Estimates 


The maintenance of an LFC airplane 
is generally similar to that of a current 
jet transport. Emphasis is placed on 
maintaining clean, smooth LFC sur- 
faces, primarily a wipedown operation 
performed at the flight line. Periodic 
maintenance on such items as waviness 
checks, trouble shooting, repairing or 
replacing suction equipment, etc., re- 
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TURBULENT LFC 


Fig. 13. Maintenance cost per flight hour. 


quires no unusual skills or equipment. 
LFC malfunction is not a safety-of- 
flight item. 

Northrop has accumulated 12 years’ 
experience in the maintenance aspects of 
LFC through laboratory testing, ex- 
posure of test panels to various climates 
over extended periods of time, wind- 
tunnel experiments, and actual opera- 
tional experience on the F-94 test air 
plane. This background is being fur- 
ther expanded through the demonstra 
tion aircraft program now in progress. 

Utilizing this experience, maintenance 
requirements have been established to 
assure satisfactory LFC_ operation. 
These requirements have been applied 
to current jet transport aircraft main- 
tenance procedures as shown in Fig. 12. 

The 3.53 man-hours of maintenance 
per flight hour is in addition to the 
average figure of 15 man-hours for 
current jet transports. This is an in- 
crease of 23 percent resulting from LFC. 
Converting these maintenance man- 
hours and material costs to maintenance 
costs per flight hour (using the 1960 
ATA formula), the results indicate only 
a 9 percent cost increase attributable to 
LFC (Fig. 13). 

This 9 percent is accounted for in the 
direct operating costs of the LFC air- 
planes shown in the Design Studies 
section which follows. 


Design Studies 


Results of the research, analytical, 
and design efforts to date form the basis 
for this section of the paper, which com- 
pares performance and direct operating 
costs of subsonic turbulent and LFC 
jet transports. 

Historically, as shown in Fig. 14, the 
rate of reduction in the skin friction drag 
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Fig. 14. Skin friction coefficient vs. time. 
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coefficient has come to a _ standstill; 
current jet transports have reached the 
lowest drag levels feasible for turbulent 
airplanes. The coefficients shown are 
for the complete airplane 

The limiting condition, as indicated 
by the dotted line of Fig. 14, is the drag 
level of a smooth turbulent flat plate at 
approximately 40 million Reynolds 
number. The fact that a turbulent jet 
transport, with its curvatures and inter- 
sections, manages to approach this 
flat plate limiting condition so closely is 
attributable to clean aerodynamic de- 
sign and to local areas of laminar flow on 
the airplane. 


DRAG 


FRICTION DRAG 
(LFC) 


WING AREA 
Fig. 15. 


Optimum wing area. 


The only method substantially 
changing the trend of friction drag 
coefficients is through LFC. The jet 
transport with laminar wings and em- 
pennage, potentially available in 1965, 
shows a friction drag of 60 percent of 
its turbulent counterpart 


400,000 
a4 
300000 
z 
200,000 
CRITICAL FIELD LENGTH 8000 FT 
50,000 LB PAYLOAD y 
: RUISE= 0.80 
Fig. 16. Parametric comparison. 
One basic result of LFC is a greater 


optimum wing area for the same take-off 
weight. The interrelation of wing area, 
friction drag, and drag due to lift is 
illustrated in Fig. 15 

The drag due to lift, primarily a func- 
tion of span loading, decreases with in- 
creasing wing area, for fixed aspect ratio 
and fixed take-off weight (Fig. 15). 
But friction drag (a direct function of 
area) increases as the wing area in- 
creases; however, this increase is at a 
much slower rate for LFC airplanes be- 
cause of the lower skin friction coeffi- 
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Fig. 17. Derivation of cost data, 


cients associated with laminar surfaces 
The intersections of the drag due to lift 
curve and the friction drag curves repre 
sent the points of maximum lift-to-drag 
ratio for the respective airplanes, with 
the LFC airplane achieving its mayi- 
mum lift-to-drag ratio at a lower drag 
level and larger wing area than the 
turbulent airplane. This effect raises 
the typical lift-to-drag ratio approxi- 
mately 50 percent. 


Parametric Aircraft Analysis; Variable 


Gross Weight 


In examining the various benefits of 
LFC as applied to jet transports, a para- 
metric study was made to determine the 
direct operating cost variation as a 
function of design range for both tur- 
bulent and LFC transport aircraft 
The study utilized a design payload of 
50,000 Ibs. at the design range. A 
maximum payload capability of 70,000 
lbs was allowed for at the original gross 
weight of the design payload-range point 
by off-loading fuel. 

As is typical of most parametric 
analyses, the assumption of scaled en 
gines was included. Turbofan engines 
with a thrust-to-weight ratio of 4.25 and 
a cruise sfc of 0.78 were used. These 
engine characteristics are typical ol 
current values. The wing areas are 
primarily a function of cruise require- 
ments rather than take-off requirements 
although no configuration used in the 
study exceeded an 8,000-ft critical field 
length. The airplanes were similar in 
appearance to current jet transports 
and all airplanes had the same siz 
passenger-type fuselage. The LFC ait- 
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ig. 18. Direct operating costs, parametric 
analysis. 
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BALLISTIC MISSILE AND 
SPACE VEHICLE SYSTEMS 


Edited by H. S. Serrert, Stanford Univ., and KENNETH 
Brown, John Wiley & Sons, Inc. Covers engineering 
design principles of missiles and space vehicles which 
resemble each other in everything but the trajectories 
followed. Attention is focused on engineering rather 
than on biological or space physics problems. It pro- 
vides a thorough survey of the ballistic problem, in- 
cluding vehicle stabilization and re-entry dynamics 
and support system planning. The vital relationships 
among the many complex parts of a missile or space 
system are given great stress. 

The whole is based on a course taught by a team, 
each member of which is now working at high levels 
in important space and missile programs. 


1961. 526 pages. $12.00 


STATISTICAL ANALYSIS AND 
OPTIMIZATION OF SYSTEMS 


By E. L. Peterson, General Electric Co. Analysis, syn- 
thesis, and optimization of systems where statistical 
uncertainty is involved in the process dynamics. 
Gives analytical foundations, describes computation 
procedures, and numerical results of specific examples. 

1961. 190 pages. $9.75 


ANALYSIS OF NONLINEAR 
CONTROL SYSTEMS 


By Dunstan Graunam, Princeton University and DUANE 
McRuer, Systems Technology, Inc. Key word is analy- 
sis. Offers basic mathematical tools for problems of 
stability, accuracy, and response of nonlinear systems. 
Illustrated with practical problems, design data. 
1961. 482 pages. $9.75 


TEMPERATURE MEASUREMENT 


IN ENGINEERING, Volume II 
By H. D. Baker, Columbia Univ., E. A. Ryper, United 
Aircraft, and N. H. Baxer, Columbia Univ. Methods 
for measurements in situations arising in shop, labora- 
tory, and field. 1961. 510 pages. $13.00 
Vol. I, 1953. 179 pages. $4.95 


PROCEEDINGS OF THE FIRST 
SYMPOSIUM ON ENGINEERING 
APPLICATIONS OF RANDOM 
FUNCTION THEORY AND PROBABILITY 


Edited by J. L. Bocpanorr and F. Kozin, Purdue Univ. 
Record of the first symposium devoted to general engi- 
necring applications of stochastic processes. 


1961. In press. 


DESIGN DATA IN 
AERONAUTICS AND ASTRONAUTICS 


Edited by Ricuarp B. Morrison, assisted by Metva J. 


INGLE, both of University of Michigan. Brings together 
information previously unpublished or scattered 
widely throughout the literature. All data are pre- 
sented in forms convenient for analysis and computa- 
tion relating to industrial, research, teaching, and 
military needs. Correlated data and nondimensional 
parameters are given wherever possible, and graphs 
rather than tables are used when it is desirable to show 
trends. 
Special feature—a Mollier Diagram for Equilibrium Air. 
Includes: Atmospheric properties; thermodynamic 
properties of gases; fluid flow properties; aerody- 
namics; aerothermodynamics; performance; ma- 
terials; aerothermochemical data; human tolerances. 
1961. Approx. 608 pages. Prob. $12.50 


THE MECHANICS OF INERTIAL 
POSITION AND HEADING INDICATION 


By W. Markey and J. Hovorxa, M.I.T. A concise, 
rigorous treatment of the synthesis of inertial naviga- 
tion systems, which can be immediately applied to 
practical work in the field. 1961. 92 pages. $3.95 


A GUIDE TO 
FORTRAN PROGRAMMING 


By D. D. McCracken, McCracken Associates, Inc. First 
practical book on FORTRAN, the computer language 
designed primarily for scientific and engineering pro- 
gramming. A self-study book for those who want to 
program problems but don’t care about computer 
mechanics per se. Includes exercises. 

1961. 88 pages. $2.95 


OPTIMUM DESIGN OF 
MECHANICAL ELEMENTS 


By R. C. Jounson, Yale Univ. Introduces a method of 
design which is explicit in nature, can be applied to 
innumerable problems, minimizes cut-and-try, and 
takes into account inherently unavoidable limitations 
confronting the design engineer. 


1961. 535 pages. $11.50 


INTRODUCTION TO 
SPACE DYNAMICS 


By W. T. THomson, UCLA. Stresses analytical tech- 
niques in the areas of coordinates for space motion; 
orbits in central force field; gyroscopic theory applied 
to instruments, missiles, satellites; generalized theory 
of dynamics from a variational approach; optimiza- 
tion techniques applied to rockets and trajectories. 
1961. Approx. 352 pages. Prob. $11.50 


Send now for your on-approval copies 


JOHN WILEY & SONS, INC. 440 Park Avenue South, New York 16, N.Y. 
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planes had LFC on the wing and em- 
pennage. 

The results are shown in Fig. 16, 
where each curve represents a series of 
airplanes, with each airplane optimized 
for maximum range for the weight 
shown. A 63 percent range improve- 
ment for the LFC airplane is shown at 
260,000 Ibs, a take-off weight typical of 
current jet transports. In comparing 
both airplanes at the same range of 
4,000 nautical miles, the LFC airplane 
shows a take-off gross weight 19 percent 
lower than that of the turbulent air- 
plane. 

Direct operating costs as a function of 
range were calculated by 1960 ATA 
methods for each of five different gross 
weights for both the turbulent and LFC 
airplanes. These direct operating costs 
were plotted separately for the turbu- 
lent and LFC airplanes as represented 
by the numbered curves of Fig. 17. 

Assuming that the curves of Fig. 17 
typify turbulent aircraft, the cost curves 
for the LFC airplanes have the same 
general characteristics but are of a lower 
cost magnitude. In both cases, turbu- 
lent and LFC, envelope curves were 
drawn through the design payload points 
as shown by the solid line, and through 
the maximum payload points as indi- 
cated by the dashed line. It is evident 
that minimum direct operating costs are 
achieved at less than design range and 


Fig. 19. 


Planform comparison. 


at the maximum payload; but in the 
cost comparisons that follow, the design 
payload values are used to maintain a 
consistent relationship with the design 
ranges quoted. Cost reduction for 
LFC is relatively constant for both de- 
sign and maximum payloads. 

The envelope cost comparison, de- 
rived by the method of Fig. 17, for the 


PAYLOAD [OOO LB 


Fig. 20. Range comparison fixed weight. 
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turbulent and LFC aircraft at the design 
payload condition is shown in Fig. 18. 
The LFC aircraft direct operating costs 
vary from 14 to 26 percent lower than 
the turbulent aircraft over ranges of 
2,000 to 6,000 nautical miles 

It is also apparent (Fig. 18) that the 
cost reduction through an LFC airplane 
is substantial when compared to a 
turbulent airplane of the same weight 
for example, a 13 percent reduction at a 
gross weight of 300,000 lbs 

As impressive as these cost reductions 
are through LFC, it must be appreciated 
that these curves show optimized air- 
planes for each specific range, on a 
parametric basis, and therefore do not 
necessarily represent the practical opera 
tional case due to the unavailability of 
fractionally sized propulsion units. 

In summary, it has been shown that 
when comparing an optimized LFC air 
plane to an optimized turbulent airplane 


6 
2 
< 12 
8 & 
10 
= 
PAYLOAD 
(MAXUM) 
DESIGN TAKEOFF WT 267000 LB 
2 CRITICAL FIELD LENGTH 6000 FT 
Mcauise = 0.80 
RANGE X 1000 NMI 
Fig. 21. Direct operating costs fixed weight. 


for the same design payload-range con- 
dition, the LFC airplane provides a 
significant reduction in direct operating 
costs. 


Specific Aircraft Analysis; Fixed Gross 
Weight 


Cargo Transport 


The parametric analysis showed that 
the optimum LFC airplane would be 
somewhat smaller in weight and would 
exhibit superior operating cost char- 
acteristics as compared to the optimum 
turbulent design. Consider now an 
existing airplane modified by incor- 
porating LFC on a new wing and 
empennage. Acontemporary turbulent 
design that meets the USAF SOR-182 
requirements is considered here because 
of its familiarity to the industry. This 
design uses programed turbofan engines 
and meets the 6,000-ft critical field 
length requirement. The payload- 
range requirements are 50,000 lbs for 
4,000 nm and 20,000 Ibs for 5,500 nm. 

The existing and modified airplanes 
are shown in the planform comparison 
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PRYLOAD x 1000 LB 
poc $/nm 


RANGE ~1000 NMI 


694 
RANGE ~1000 Nm 


Fig. 22. 


LFC cost reduction. 


of Fig. 19. The fuselages, wing geom 
etry, and tail volumes are identical for 
each airplane. Also, they each have 
the same cruise speed of 0.80 Mach 
number. 

The wing area on both airplanes is a 
compromise between take-off and maxi 
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Fig. 23. Skin temperature. 


mum range requirements. The take-off 
characteristics require that the wings 
are oversized, such that a 1 to 2 percent 
range reduction occurs. The lower 
drag of the LFC airplane allows the use 
of engines of 30 percent lower thrust 
than the turbulent airplane, but the 
take-off requirements can be main 
tained on the LFC airplane because of 
its larger wing area. The same take-off 
lift coefficients are used in all cases. 
The payload-range characteristics of 
the two airplanes are shown in Fig. 20. 

The LFC airplane carries a 50,000-lb 
payload approximately 30 percent far- 
ther and the 20,000-lb payload 33 per- 
cent farther than the turbulent airplane 
at the same gross weight. 

Fig. 21 compares the 1960 ATA direct 
operating costs as a function of range 
for the two aircraft. The LFC airplane 
is slightly more economical below the 
2,500-nm range point. At the design 
payload-range point of 4,000-nm the 
LFC airplane shows a 23 percent re- 
duction in direct operating costs com- 
pared to the same weight turbulent air- 
plane. The operational advantages of 
the LFC airplane, carrying the same 
payload as the turbulent airplane, are 
shown by the shaded area of Fig. 21. 
For trans-Atlantic or Hawaii-to-Japat 
routes (approximately 3,500 nm), the 
cost reduction is approximately 18 per- 
cent for the LFC. At ranges greater 
than design range, such as a flight from 
Los Angeles to Paris (approximately 
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5,500 nm), the cost reduction with the 
LFC is dramatic—a 43 percent reduc- 
tion as compared to the direct operating 
cost of the turbulent airplane. It is 
also evident from Fig. 21 that the LFC 
airplane carries the design payload 30 
percent farther than the turbulent air- 
plane and at a lower cost. 

In summary, incorporation of LFC 
on an existing airplane can provide low 
operating costs over a broad spectrum 
of ranges. Small savings are available 
below the maximum payload point, and 
increasingly greater savings are avail- 
able above this point. This ability to 
operate over greater ranges at low cost 
with LFC increases operational flexi- 
bility and precludes the necessity of 
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Fig. 24. Drag breakdown. 


purchasing several sizes of aircraft to 
obtain low operating costs at the various 
ranges between 2,000 and 6,000 nm. 


Passenger Transport 


Consider now the economic advis- 
ability of LFC application to a con- 
temporary turbofan passenger transport. 
Assume, in this case, that LFC is in- 
stalled on the wing alone. From the 
payload range curve (Fig. 22), it is seen 
that with a nominal payload of 30,000 
lbs the range is extended about 33 
percent. This enables the airplane to 
operate on several routes that were 
previously impossible even with zero 
payload. 

Fig. 22 also compares the direct 
operating costs in dollars per nautical 
mile as a function of range for the turbu- 
lent airplane and its LFC counterpart 
with LFC on the wing alone. For a 
payload of 30,000 lbs and 10 hours per 
day utilization, a cost reduction of 
approximately 6.5 percent over a wide 
spectrum of operating ranges results for 
the LFC installation. This 6.5 percent 
reduction in operating costs represents a 
saving of $2.1 million over a 10-year 
period. If it is assumed that the modi- 
fication cost for a new wing and LFC 
suction system is about $1.0 million, 
the earning power of the passenger 
transport is increased by $1.1 million. 

In summary, it can be seen that 
through the application of LFC to the 


wing alone on a passenger-type trans- 
port, the operational flexibility is in- 
creased, the direct operating cost is re- 
duced, and over a wide spectrum of 
operating ranges the earning power of 
the airplane is significantly increased 
during its operational lifetime. 


Future Trends 


LFC will influence the future trends of 
subsonic and supersonic jet transport 
designs to a much greater degree than 
generally envisioned today. Just as 
the introduction of cantilever mono- 
planes made it worthwhile to cowl 
engines and retract landing gears, LFC 
will cause evolutionary changes in the 
design of airplanes to augment the 
basic advantages of laminar flow. 
These changes could substantially im- 
prove range and operating costs. 


Subsonic 


Initially, subsonic LFC jet transports 
are expected to resemble current designs 
for traditional reasons but will grad- 
ually change in appearance while be- 
coming more efficient because (1) in- 
duced drag will form a larger part of the 
airplane drag and thus pose a challenge 
to develop wings of greater span while 
retaining present structural-weight effi- 
ciencies; (2) further aerodynamic 
cleanup will be worthwhile (smoother 
windshields, fewer bumps and fairings, 
lower interference drag, and more highly 
streamlined bodies); and (3) further 
reductions in friction drag will be pos- 
sible through application of LFC to the 
fuselage and other turbulent areas. 

The investigation of laminar flow on 
bodies, including bodies of revolution, is 
currently the subject of a research 
program at Northrop under Air Force 
contract. This program includes wind- 
tunnel testing at high Reynolds num- 
bers and will help determine the prac- 
ticality of LFC on bodies. 

If LFC is applied to the major portion 
of the fuselage and engine nacelles of 
the jet transport shown in Fig. 16 (which 
had a design payload turbulent range of 
4,000 nm and range of 6,500 nm with 
LFC on the wings and empennage), the 
range would be approximately 7,400 
miles, or 85 percent greater than the 
original turbulent airplane. 


Supersonic 


The feasibility of supersonic LFC 
was demonstrated in 1958 with wind- 
tunnel tests at Mach numbers 2.23 and 
2.77 on a 5 percent biconvex airfoil at a 
relatively low chord Reynolds number 
of 12.5 million.? Additional research is 
now in progress at Northrop under Air 
Force contract to investigate LFC on 
bodies and wings at high Reynolds 
number in wind-tunnel experiments 
from Mach 2.5 to 3.5. 
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Maintaining laminar flow on a super- 
sonic transport has two potential advan- 
tages, lower skin temperatures and re- 
duced friction drag. The cooler skin is 
due to a large reduction in the con- 
vective heat-transfer coefficient re- 
sulting from LFC (essentially a function 
of the lower velocity gradient in the 
laminar boundary layer), coupled with 
the effect of surface radiation. The 
temperature reduction is shown in 
Fig. 23. 

All temperatures in Fig. 23 are 
equilibrium values and are reached in 
approximately 10 min. The example 
chosen for analysis was Mach 3.0 at an 
altitude of 75,000 ft. The calculations 
are based on simple flat plate heat- 
transfer coefficients for turbulent and 
laminar flow. LFC skin temperatures 
are somewhat higher than natural 
laminar flow, asshown. The suction air 
will have a total temperature only 
slightly higher than the external laminar 
surface temperature, so the duct wall 
temperatures will approximate those of 


Fig. 25. 


M = 3.0 transport comparison. 


the external surface. Results of Fig. 23 
may be summarized as follows: 

(1) Application of LFC to the wing 
could reduce the equilibrium surface 
temperature by as much as 100°F 
from the turbulent boundary layer 
value. 

(2) Equilibrium surface temperature 
over a large portion of the wing might 
be maintained below 350°F with LFC 
for M = 3.0 flight. This fact intro- 
duces the feasibility of construction 
materials other than steel. 

(3) Application of LFC to the fuselage 
appears advantageous. In addition to 
further reducing the total vehicle drag, 
reducing the structure equilibrium tem- 
perature, and lowering the boundary 
layer noise, the cabin conditioning heat 
load is substantially reduced. This 
reduction would be about 35 percent of 
the turbulent flow heat load. 

Friction drag continues to be a sub- 
stantial percentage of the total airplane 
drag at Mach 3.0 and probably will in- 
crease in percentage with time as new 
methods develop to reduce wave drag. 
A typical drag breakdown for a Mach 
3.0 transport, representative of airplane 
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types presented at the recent IATA su- 
personic transport conference, is shown 
in Fig. 24. For a fixed range of 3,500 
nm, the LFC wing area remains ap- 
proximately equal to that of the turbu- 
lent wing. This similar wing area pro- 
vides the optimum lower wing loading 
for the LFC airplane as a function of 
less fuel and less propulsion system 
weight. Wave drag remains constant 
with the fixed percentage thickness 
ratio of the airfoil and the fixed wing 
area. Since drag due to lift is reduced 
only slightly, in proportion to the wing 
loading, the LFC reduction in friction 
drag contributes most significantly to 
the 24 percent reduction in the total 
airplane drag. The friction drag re- 
duction shown assumes LFC on the 
wing and tail only. The resulting per- 
formance improvement shown in Fig. 
25 includes a weight penalty for the 
LFC system. 

At a design range of 3,500 nm the 
weight of the Mach 3.0 transport could 
be reduced by approximately 60,000 
Ibs, or the range could be increased 
by 700 nm. The most advantageous 
choice is considered to be the weight 
reduction from several standpoints 
lower cost, present runway compatibil- 
ity, take-off and landing speeds com- 
mensurate with current subsonic jet 
transports, and less sonic boom (lower 
wing loading). 


Studies indicate the desirability of 
continuing the present research to ad- 
vance supersonic LFC technology to the 
same high confidence level now existing 
at subsonic speeds 


Conclusions 


Application of LFC technology is 
desirable for the following reasons: 

(1) LFC is feasible from the stand- 
point of demonstrated drag reduction 
and practicality of construction and 
maintenance. 

(2) LFC provides for the possibility of 
appreciably lower operating costs in the 
future development of subsonic trans 
ports. 

(3) Preliminary examinations show 
that lower skin temperatures and lower 
take-off weights are 
through the application 
supersonic transport 
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High Altitude Wind Observations (Continued from page 17) 
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intervals. The data were interpolated 
at 500-meter intervals of altitude. 
Shear was then computed for each 1,000- 
meter altitude layer and attributed to 
the midpoint of the layer. 

Care should be exercised in inter- 
preting the reduced wind data because 
smoothing methods used in the com- 
puter operations tended to reduce cer- 
tain wind speed and shear values. Also, 
in some cases the use of chaff targets 
gave abnormally high wind speeds and 
shears. The chaff targets dispersed 
wer an area of several square kilom- 
eters, and the radar focused on the 
point of strongest signal return which 
ften shifted from one point to another. 
Such abrupt shifts caused abnormal 
speeds, shears, and directions to appear 
in the reduced data, thereby requiring 
careful editing. 

Simultaneous tracking of the rocket- 
sonde targets by two or more radars has 
made a comparative check of the com- 
puted wind values possible. Similar 
comparisons of the overlapping portions 
of the rawinsonde wind measurements 
show that there exists an uncertainty in 
the rocketsonde wind speed of approxi- 
mately 8 meters/sec or 10 percent of the 
computed wind speed, whichever is 
greater. With careful investigation and 
correction of the more radical data 
points measured, averages should be 
reasonably accurate even though indi- 
vidual data points are questionable. 
This agrees favorably with other com- 
parable error analyses. ** 
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vailing wind by rocketsonde 
observation, April 12, 1959, 
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Presentation 


While data from 37 soundings were 
used, not more than 24 observations 
were available for any single altitude 
level. When fewer than five observa- 
tions were available for a given altitude 
level, the data are shown on the graphs 
by dashed lines. Data were sum- 
marized for the winter months of 
October through March, for the summer 


(hen) ALTITUDE 


| 
| 
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months of April through September, and 
for the entire year. Because of dis- 
sipation of the chaff targets, few ob- 
servations extended below 20-km al- 
titude. No observations were made at 
altitudes above 84.5 km; consequently, 
most of the graphical data are confined 
to the 20- to 66-km altitude range. 


Analysis of Data 


Inspection of the arrowgram (Fig. 1) 
shows that winds from any direction 
can occur at almost any altitude or 
during any season, at least occasionally. 
However, when the measured winds are 
summarized (Fig. 2), it becomes ap- 
parent that easterly winds prevail 
during summer months and westerly 
winds prevail during winter months. 
On the rare occasions when easterly 
winds do occur in winter months or 
westerly winds in summer months, 
they seem to be characterized by low 
velocity at altitudes up to approx- 
imately 70 km. Annual prevailing 
winds—1.e., those directions which occur 
most frequently—fluctuate between 
2asterly and westerly according to the 
data presented here, probably because 
more observations were made during 
summer months when easterly winds 
prevail than in winter when westerly 
winds prevail. If an equal number of 
observations had been made during 
winter and summer months, it is likely 
that the annual prevailing winds would 
be predominantly westerly at almost 
all altitudes. 

A minimum wind speed was observed 
at an altitude of 23 km® (Fig. 3). Above 
this point, mean wind speed increases 
with altitude up to at least 64 km, and 
possibly higher, to another region of 
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‘wind component“ 
maximum wind speeds apparently cen- ALTITUDE, 
tered near 70 km.°’ The maximum OBSERVATIONS 
wind speed observed in this rocketsonde e¢}-— = 
series was 133 meters/sec from the 
southwest at an altitude of 71 km; 
it occurred in winter. The summer ee 2 
maximum was 99 meters/sec from the s 5 
east-southeast at 69-km altitude. Win- 
ter winds were considerably stronger 
than summer winds, except at the sek yr 3 
wind speed minimum which occurs at re 4 
an altitude of approximately 23 km. 4 20 
44} < 2! 
The stronger winter winds were es- 
pecially noticeable at higher altitudes. sak 
These observ ations agree favorably with oesenven 
the findings of other studies.® 20 Aya ANNUAL MEAN 17 
Easterly winds first appear at al- wiehesr oaserver 4 
titudes above 50 km in late April, 
4 0.01 0.02 0.03 0.04 0.05 0.06 
and westerly winds first occur in late SHEAR (v0c"!) 
October. Unfortunately, no observa- Fig. 5. High-altitude wind shear by rocket- 


tions were made in early October, and 
very few were made during the month of 
April. More observations will be re- 
quired to confirm the date at which the 
reversal of prevailing wind occurs; 
however, it appears likely that this re- 
versal occurs within a very few days. 

Meridional wind components are 
generally small, and the annual merid- 
ional wind resultant is nearly zero, 
although southerly winds seem to occur 
slightly more often at altitudes below 
57 km (Fig. 4). Seasonal meridional 
wind components reverse so frequently 
that they are omitted. 

The summer zonal wind components 
are easterly, and the winter zonal wind 
components are westerly and much 
stronger than the summer easterlies 
(Fig. 4). While approximately 60 
percent of the observations were made 
during the summer months and 40 
percent were made during the winter, 
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sonde observation, April 12, 1959, to May 25, 
1960, Cape Canaveral. (Shear Ah = 1,000 
meters.) 


the combined resultant for the annual 
summation is mainly westerly because 
the wind speed is greater for the westerly 
winds. 

The wind shear, the partial derivative 
of the wind vector with respect to 
altitude, averages less than 0.01 inverse 
sec/1,000-meter altitude layer up to 
60 km (Fig. 5). The highest shear in 
this study, 0.0523 inverse sec/1,000- 
meter altitude layer, was measured at an 
altitude of 68 km. Because of errors 
inherent in the method of observation,° 
the values of shear data obtained might 
not be exact, but there is evidence to 
substantiate an increase in the mag 
nitude of wind shear at higher altitudes. 


October 1961 


Conclusion 


The mean rocketsonde wind speg 
observed in the high altitudes aboy 
Cape Canaveral ranges from 5 meters 
sec at an altitude of 20 km to 4 
meters/sec at an altitude of 58 km 
Individual observations range from 
nearly zero to 133 meters/sec. Winds 
from any direction can occur at almost 
any altitude or during any season, but 


summer winds predominantly 
easterly and winter winds are pre 
dominantly westerly and 


stronger 
Meridional wind components are small 
and zonal wind components are westerly 
and strongest in winter. Wind shear 
ranges from nearly zero to more than 
0.05 inverse sec/1,000 meters and in. 
creases in magnitude with increases in 
altitude. At altitudes up to 59 km 
the mean wind shear is less than 0,0) 
inverse sec/1,000-meter layer. 

The data on high-altitude wind 
characteristics presented here are in- 
tended for the use of design and opera- 
tions personnel concerned with missile 
and space vehicle test flights from Cape 
Canaveral (AMR), Florida. 
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Personnel Uoportunines 


This section is for the use of individual members of the Institute 
seeking new connections and eligible organizations offering em- 
ployment to specialists in the aero/space industry. Any member 
or eligible organization may have requirements listed without 
charge by writing to the Secretary of the Institute. 


Wanted 


Engineers and Scientists—Electronics Re- 
search Directorate of the Air Force Cam- 
widge Research Laboratories announces 
professional staff vacancies for Engineer 
Electronics); Physicist (Electronics, Elec- 
tromagnetics, Plasmas, Solid State); Math- 
ematician; Chemist (Physical, Analytical, 
Nuclear). These vacancies are within the 
Federal Civil Service and offer liberal leave 
benefits, payment of expenses incurred by 
employees and their families in moving to 
Bedford, an extensive program of sponsored 
education, and an opportunity to pursue a 
stable program of non-time-oriented research. 
Graduate level scientific personnel inter- 
ested in positions described are invited to 
write or call Francis McWilliams or John 
Cooney, at CRestview 4-6100, Ext. 3633 
or 3758 at Hanscom Field, Mass. 

Professors—The Air Force Institute of 
Technology, Wright-Patterson AFB, Ohio, 
has a vacancy in the Department of Mathe- 
Most of the work is at advanced 
undergraduate and graduate levels. One 
quarter out of four free of teaching duties for 
research or other academic pursuits, plus 
vacation. Working conditions comparable 
) those in leading universities. Employ- 
ment will be effected in accordance with 
Civil Service regulations. Rank and salary 
will depend upon qualifications of applicant 
and may be anywhere between assistant 
professor, GS-11, $7,560, and full professor, 
GS-14, $12,210. Applications should be 
made on Standard Form 57, available at any 
Post Office or by letter to Head of the De- 
partment of Mathematics. 


Aerodynamicists—Cornell Aeronautical 
Laboratory has attractive openings for 
qualified scientists to perform research in the 
following areas: (1) magnetoaerodynamic 
effects on the flow over blunt bodies; (2) 
flow field analyses at hypersonic velocities 
and low densities; (3) aerodynamic char- 
acteristics of simple geometry hypersonic 
configurations; (4) aerodynamic aspects of 
air-data sensing at hypersonic flight condi- 
tions; (5) base heating of clustered rocket 
nozzle configurations; (6) structure of air 
heat transfer due to turbulent boundary 
layers at hypersonic Mach numbers. The 
Laboratory provides a desirable combination 
of the best features of industry as to financial 
and facilities resources coupled with a 
university associated atmosphere conducive 
toR&D activity at the highest level. Salary 


matics. 


Notices of all change of ad- 
dress should be sent to the Circu- 
lation Department, Institute of the 
Aerospace Sciences, Inc., 2 East 
64th Street, New York 21, N.Y., 
at least 30 days prior to change 
of address. 


The number preceding the notice 
represents the Box Number of the 
Institute of the Aerospace Sciences to 


which inquiries should be addressed. 


commensurate with qualifications and ex- 
perience. Send complete résumé to A. A. 
Rentschler, Cornell Aeronautical Labora- 
tory, Inc., 4455 Genesee St., Buffalo 21, N.Y. 

138. Professor of Aerospace Engineer- 
ing—Positions now available in department 
of aerospace engineering of state university. 
Openings are in areas of aerospace structures, 
aerodynamics, and magnetoaerodynamics. 
Rank and salary will depend upon qualifica- 
tions of applicant. Liberal research and 
consulting policies. Furnish résumé with 
reply. 


Available 


137. British Aircraft Designer (AFIAS)— 
Interested in consultation and formation of 
design team in United Kingdom. Would be 
prepared to consider any design consign- 


ment. Work would be completed rapidly 
and at nominal cost. Résumé on applica- 
tion. 


136. Aeronautical Engineer—Thirty 
years’ experience—civil approval of power 
plants, power-plant supervisor with aircraft 
manufacturer, chief of AEDC engine test 
facility; also manager, engineering sales and 
advanced proposals. Desires engineering 
applications and Government contacts for 
advanced propulsion and space systems. 
B.S. in M.E. and LL.B. | Ohio resident. 


134. Project or Engineering Manager— 
Broad background in aeronautical and ballis- 
tic missile design and management. Ex- 
perience in establishing structural design, 
criteria, performance of operational analysis, 
and weapon system analysis and develop- 
ment. Considerable experience as project 
manager handling all phases of design, 
budgeting, program planning, and supvport 
activities for an operational military system. 
Age 44. AFIAS. Reg. Prof. Eng. Desires 
position in program planning, project man- 
agement, engineering management or at 
staff level. 

133. Aeronautical Engineer—Thirteen 
years’ experience in high-speed wind-tunnel 
testing, data reduction, and data analysis. 
Ten years in a supervisory capacity. Ex- 
perience includes research activities, tech- 
nical report preparation, and aerodynamic 
data analysis. Work also accomplished in 
the field of plasma generation in connection 
with arc-heated wind-tunnel development 
and employment as a testing facility. West 
Coast location desired; other locations con- 
sidered. Résumé furnished on request. 
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132. Director—Research, Engineering, or 
Planning—Age 39, M.S. in A.E., Registered 
P.E., AFIAS. Eighteen years’ extensive 
experience in scientific and management 
fields as consultant, director of research, and 
general manager with well-known company in 
aerospace and related fields. Has proved ca- 
pabilities in planning, organizing, and admin- 
istering successful programs—completed on 
schedule and at a profit. Familiar with 
USAF, USN, Army, and many commercial 
companies’ plans. Desires responsible posi- 
tion with dynamic growth company. 

131. Engineer—M.E., MBS. Fifteen 
years’ experience in research and develop- 
ment on electromechanical devices and sys- 
tems. Includes 9 years in design and evalua- 
tion of inertial guidance components and 
systems. Now Manager of Analysis for 
large manufacturer. 


130. Engineer - Designer - Draftsman— 
Age 50. Thirty years’ experience in draft- 
ing, design, and engineering of airframes. 
turbine engines, automotive parts, fans, 
blowers, superchargers, airmotors, valves, 
and fittings. Two years in aircraft engine 
stress analysis. Registered Professional En- 
gineer (M.E., Indiana). Salary open. Ré- 
sumé on request. 


129. Aero/Mechanical Project Engineer 
—Seeking position with a growing, stable 
company. Over 10 years’ diversified back- 
ground as a project, development, and appli- 
cation engineer in aircraft structural and me- 
chanical design; electromechanical, hydrau- 
lic, and pneumatic controls; rotating machin- 
ery including small turbines. Familiar with 
control dynamics, filament wound fiberglass 
structures, and cryogenics. B.S. Aero Engi- 
neering, Post Graduate Industrial Manage- 
ment. Former USAF and civilian pilot. 
Prefer a medium-sized equipment company 
in Northeastern coastal location. All replies 
considered. Résumé on request. 


An Authorized Binding for 
Aerospace Engineering 


Arrangements have been made 
by the Institute of the Aerospace 
Sciences for members and sub- 
scribers to have their magazines 
bound carefully and economically. 
Simply write for details to: 


Publishers’ Authorized 
Bindery Service Ltd. 
430 West Erie St. 

Chicago 10, Ill. 
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Early Blind Flying (Continued from page 15) 


the magnetic compass, the altimeter, and the air- 
speed indicator. 


It was customary on cross-country flights to 
follow railroads, or less frequently, highways, and 
the maps generally used were the standard Rand 
McNally maps of the individual states. These 
state maps were each about the same size so, un- 
fortunately, were usually not to the same scale. 

Please excuse me if in the following I seem to talk 
too much about myself. The personal experiences 
are related to illustrate conditions which affected 
all fliers of the era. 

On September 4-5, 1922, I flew a DH-4 airplane 
in which additional gasoline and oil tanks had been 
installed from Pablo Beach, Fla.—near Jackson- 
ville—to Rockwell Field at San Diego, Calif., with 
one intermediate stop for fuel at Kelly Field, 
San Antonio, Tex. The elapsed time was 22 hours 
and 35 minutes. It was the first time the North 
American Continent had been crossed in less than 
24 hours. 

It was also the first airplane in which I had used a 
bank and turn indicator. To obtain the instrument, 
it was necessary to go the the Army Air Service 
Engineering Base at old McCook Field in Dayton, 
Ohio, and “promote” an experimental model 
through the help of cooperative technical friends. 
This instrument, invented in 1917 by Elmer Sperry, 
Sr., built by Elmer Sperry, Jr., and first flight tested 
by Lawrence Sperry with Elmer, Jr., as passenger, 
was not yet in common usage or generally available. 

Took off just after dark and arrived at Kelly 
Field soon after daylight. Was quickly refueled 
and immediately proceeded to Rockwell, arriving 
there in the late afternoon. Had chosen a moon- 
light night to facilitate night flying, but about 4 
hours out ran into solid overcast and then severe 
thunderstorms. 
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“On September 4-5, 1922, | flew a 
DH-4...It was the first time the North 
American continent had been crossed in 
less than 24 hours." (IAS Photo, gift of 
Albert R. Crocker, from the Sherman Fairchild 
Collection of Aeronautical Photographs) 


For a while the lightning flashes were alimost 
constant and, in the otherwise black night, so in- 
tense as to light up the ground clearly for a con- 
siderable area. Some flashes were so close that 
their familiar ozone odor could be detected, but, 
although it seemed that one could reach out and 
touch them, none struck the plane. 

The air was extremely turbulent and the air- 
plane was violently thrown about its axes as well as 
up and down and, despite its excellent stability 
characteristics, was held on a relatively even keel 
only with great concentration and effort. After 
the lightning died away, the turbulence appeared to 
intensify and there was about an hour in the jet 
black darkness when no ground reference point 
could be seen; it would have been quite impossible 
to maintain proper attitude and course without the 
blessed bank and turn indicator. 

Although I had been flying for almost five years 
‘‘by the seat of my pants’’ and considered that I'd 
achieved some skill at it, this particular flight made 
me a firm believer in proper instrumentation for 
bad weather flying. 

In 1925 I did a thesis here at M.I.T. for a Doctor 
of Science degree in aeronautics. 

At first I had hoped to study, through carefully 
controlled flight tests, the possibility—or impossi- 
bility—of a specially trained pilot orienting himself 
without flight instruments and determining certain 
phenomena, such as wind direction, without reier- 
ence to a visible horizon or point on the ground. 

The original subject matter was not sufficiently 
abstract for the doctorate, so it was changed to 
an analysis of ‘The Effect of the Wind Velocity 
Gradient” employing flight tests, wind-tunnel data, 
and mathematical analyses. 

The thesis starts: 


There has long been an uncertainty in the minds of avia- 
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Forging 
a metal 
spectrum for 
any flight 
environment 


Whether design criteria be Mach 2 or escape 
velocity — the part a wing spar or rocket mo- 
tor nozzle — most of the significant advances 
in forging conventional, exotic and refrac- 
tory materials come from Wyman-Gordon. 
Here, technical breakthroughs in size, com- 
plexity and strength-to-weight ratios of forged 
airframe hardware and propulsion components 
have materially raised performance ceilings 
for every type of vehicle — global bomber to 
deep space probe capsule. Always impatient 
with state-of-the-art concepts, Wyman-Gordon 
engineers and metallurgists constantly explore 
ways to bring hitherto “unforgeable” metals 
and alloys into operational use — broadening 
the spectrum of ultimate-performance ma- 
terials available for advanced vehicle designs. 
Consultation on your project requirements 
can help achieve better correlation of materials 
to service environment. 


CURRENTLY FORGING THESE AERO-SPACE MATERIALS 


Aluminum Beryllium 

Titanium Molybdenum 

Magnesium Tantalum 

Columbium Zirconium 

Hafnium Tungsten 

Inco 718 René 41 

Astroloy Waspaloy 

Stainless steels Vacuum-melt alloys 

Nickel-base alloys lron-base, high- 

High-strength die temperature alloys 

steels 20-27 Nickel steels 

WYMAN - GORDON 

FORGINGS 

of Aluminum Magnesium Steel Titanium ...and Beryllium Molybdenum Columbium and other uncommon materials 
HARVEY ILLINOIS WORCESTER MASSACHUSETTS DETROIT MICHIGAN 
GRAFTON MASSACHUSETTS LOS ANGELES CALIFORNIA PALO ALTO CALIFORNIA FORT WORTH TEXAS 
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tors regarding the effect of the wind on the flying qualities 
of an airplane. 


Some pilots claim that it is much easier to turn into the 
wind than with it, and that at any altitude they can tell 
the wind direction by the feel of the ship in a turn and 
this even though in a dense cloud which would preclude 
the possibility of obtaining their relative motion from any 
stationary object. 


Other pilots maintain that, regardless of the wind velocity 
or the proximity of the ground, there is no difference in 
the feel of the plane when turning into the wind and when 
turning with it. They claim that any apparent dif- 
ference is due wholly to the psychological effect on the 
pilot, resulting from the difference in ground speed in 
the two cases, and if there is any difference in the ships’ 
performance, from a time altitude standpoint, it is be 
cause the pilot handled the controls differently. In other 
words, if the pilot were blindfolded, he could not tell the 
wind direction when turning, and a turn made into the 
wind would be identical with a turn made with the wind. 
This is, of course, considering the turn in relation to the 
medium in which it is being executed and not in relation 
to the curves traced out on the ground. 


There is a similar difference of opinion regarding the ef- 
fect of a strong wind on the rate of climb. Experienced 
pilots are about evenly divided, half feeling that a plane 
climbs better into the wind, and the other half feeling 
that the wind makes absolutely no difference 


Seven of the leading pilots of the day were ques- 
tioned regarding the effect of the wind on flying 
performance. The answers given were far from 
consistent, as might be expected from such a group 
of individualists. There was therefore still consider- 
able confusion—and controversy—among the ex- 
perts. 

The conclusions from the thesis, somewhat over- 
simplified, were that, in airplane flight, 


(1) There is no measurable effect in level flight, at 
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altitude, due to wind direction so long as the wind js 
steady. 


(2) There is no effect on climb due to wind except 
very near the ground, and there the wind velocity 
gradient increases the rate of climb slightly when 
flying into the wind and decreases it slightly when 
flying with the wind. 


(3) A steady wind has no effect on turning except 
very near the ground when the wind velocity gra- 
dient causes a slight tendency to settle when turning 
away from a head wind and a slight tendency to 
climb when turning into it. This is most noticeable 
in strong winds and when flying at a large angle of 
attack or at minimum power. 


SUMMING UP: A steady wind exercises no meas- 
urable effect on airplane performance at altitude— 
except, of course, on ground speed and direction of 
flight. Very near the ground, however, the effect of 
wind velocity gradient can be serious, particularly in 
the case of a heavily loaded airplane. The danger is 
increased by a strong tendency on the part of the 
pilot to pull the nose up or in beyond the most 
efficient angle of attack. This increases any 
tendency to settle and may even cause the airplane 
to stall and spin in. 


In the early and middle 1920’s the Jones-Barany 
revolving chair test was given to all military pilots 
as a part of their periodic physical examination for 
flying. Normally, this test was given with the 
pilot’s eyes open, and the flight surgeon looked 
for variations in times and amount of the rhythmic 
side-to-side movement of the eyes called nystagmus. 

In early 1926, Capt. (later Col.) David A. Myers, 


“The first important expenditures made 
by the Full Flight Laboratory were for two 
modern airplanes.” One, relates Gardner 
Lecturer Doolittle, was a Consolidated NY-2 
military training plane which mounted a 
special landing gear with long oleo action 
instead of pontoons. This is the NY-1. 
(IAS Photo, from the Sherman Fairchild Col- 
lection of Aeronautical Photographs) 
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an outstanding Air Corps flight surgeon, decided to 
augment the routine test by giving an additional 
test consisting of several rotations of the chair with 
the pilot’s eyes closed. After the rate of rotation 
became steady, a normal pilot, with eyes closed, 
could not tell which way he was turning. If the 
rate of rotation was slowed down and stabilized at a 
somewhat lower speed, the pilot thought the ro- 
tation had been stopped and, when the rotation 
actually was stopped, he thought he was turning in 
the opposite direction. 


The explanation is that man normally maintains 
his equilibrium by sight, touch, hearing, muscle, 
and vestibular sense. Touch and hearing are not 
important in flight orientation. By using the three 
remaining senses he can usually ascertain and main- 
tain his position, accurately sense the rate and di- 
rection of his motion, and generally orient himself 
with relation to the earth. Sight is by far the most 
reliable of these three senses, and when sight is lost, 
we must get our sense of balance and motion from 
the muscles and from the fluid movement sensors in 
the vestibular canals. If an individual is merely 
displaced, the fluid motions return to zero very 
rapidly, but if one is rotated it may take from 5 to 
25 seconds for the fluid motions to stop. 
this period an individual can experience a false sense 
of motion called vertigo. This, of course, explains 
why early-day pilots flying in dense fog or clouds 
frequently became completely confused and _ oc- 
casionally spun in and crashed. 

Captain (later Col.) William C. Ocker, an early 
and extremely competent Army Air Corps pilot and 
flight researcher, had long been interested in instru- 
ment flying and in 1918 had tested the then new bank 
and turn indicator. 


In mid-1926 he took Captain Myers’ new “‘blind- 
fold test.” His first reaction was that Captain 
Myers had played a trick on him or, if not, that his 
senses had failed him. After further consideration 
he decided that there was proof positive that no 
normal pilot could consistently fly ‘“‘blind’’ without 
instruments. 


During 


Ocker, who had had considerable experience flying 
with the bank and turn indicator—he frequently 
carried one, a quickly attachable unit complete with 
venturi, in his flight baggage—believed this instru- 
ment could correct the pilot’s faulty senses. He de- 
signed a lightproof “black box’’ which contained a 
bank and turn indicator and a magnetic compass. 
This box was mounted on the front of the Jones- 
Barany chair. The pilot sealed his face against the 
opening in the box and observed the bank and turn 
indicator and compass. With this piece of equip- 
ment he could correctly identify the direction and 
rate of his rotation. After the rotation stopped and 
the compass settled down, he could then determine 
heading. 


Myers and Ocker continued their experiments, 
and the arrangement of black box and revolving 
chair were patented and subsequently used in the 
training of pilots. Later some pilots were to learn 
to fly by instruments alone before they learned to fly 
under normal visual conditions. 


In the late twenties and early thirties, Capt. 
Ocker and Ist Lt. (later Col.) Carl J. Crane collabo- 
rated on the study of instrument flying techniques 
and developed, among other things, a unitary arrange- 
ment of instruments which would give the pilot a 
maximum of useful flight information with a min- 
imum of effort and fatigue. They referred to this as 
a ‘flight integrator.”’ 


In 1932, Maj. Ocker and Lt. Crane wrote Blind 
Flight in Theory and Practice (published by the 
Naylor Printing Company). This book gave an 
excellent analysis of the problems inherent in instru- 
ment flying and their solution. It was for many 
years the standard reference on instrument flying. 


Daniel Guggenheim was one of the great indus- 
trialists, philanthropists, and citizens of the 20th 
Century. He was interested in everything that could 
lead to a fuller life and a better world. One of his 
many great contributions was the Daniel Guggen- 
heim Fund for the Promotion of Aeronautics. It 
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was established for the purpose of promoting the 
advancement of the art, science, and business of 
aviation. It proved to be a very effective medium 
in the accomplishment of that purpose. 

The Fund was established in January 1926 with a 
grant of $2.5 million. (The previous year $500,000 
had been given to New York University for the pur- 
The Fund, 
in cooperation with the Government, was to be 


pose of starting aeronautical education.) 


administered by a board of trustees composed of 
men of “eminence and competence.’ Harry Gug- 
genheim, gifted son of Daniel and a World War 1 
naval aviator, was chosen president of the Fund. 
In the initial stages of its organization, Rear Adm. 
H. I. Cone, an outstanding naval officer, was vice- 
president. Then, until the Fund’s work was 
completed, Capt. (later Vice Adm.) Emory S. Land 
of the Construction Corps of the U.S. Navy held 
this position. The strength of character, sound 
judgment, organizational ability, understanding, and 
capacity for cooperation of Harry Guggenheim and 
Capt. “Jerry” Land were, in large part, the cement 
which held the Fund together, enabled it to function 
efficiently, and permitted its considerable contribu- 
tions. 

The general purposes of the Fund were defined as 
follows: 

(1) To promote aeronautical education both in 
institutions of learning and among the general 
public. 

(2) To assist in the extension of fundamental aero- 
nautical science. 

(3) To assist in the development of commercial 
aircraft and aircraft equipment. 

(4) To further the application of aircraft to 
business, industry, and other economic and social 
activities of the nation. 

The basic concept of the Fund was ‘“‘to maintain a 
simple, inexpensive directing organization depend- 
ing on established outside agencies, whenever pos- 
sible, to carry out the aims of the Fund.” It was 
to be a primer—a spark plug—to stimulate interest 
and promote action. 
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From the first it was understood that flight sufety 
and reliability were important considerations and 
that one phase of the Fund’s work might certainly 
be to study means of assuring safe and reliable flight 
despite weather conditions. With this in mind, q 
special committee of experts was organized to define 
the problem, and a directive was prepared which 
authorized study regarding 

|) Dissipation of fog. 

2) Development of means whereby flying fields 
nay be located from the air regardless of fog. 

3) Development of instruments to show accurately 
the height of airplanes above the ground so as to 
replace barometric instruments now in general use 
showing height above sea level. 

{) Improvement and perfection of instruments 
allowing airplanes to fly properly in fog. 

5) Penetration of fog by light rays. 

fo assist in carrying out pertinent parts of this 
directive, the directors of the Fund decided in 1928 
to establish a Full Flight Laboratory at Mitchel 
Field, Long Island, N.Y., with all necessary facilities 
and equipment. In August 1928, I was borrowed 
from the Army Air Corps to head this laboratory. 
Professor William G. Aeronautics 
Department of M.I.T. joined us as technical assistant 
in February 1929. Bill Brown worked with the 
Fund's Full Flight Laboratory from then until it 
1929. His technical 
knowledge and unbounded enthusiasm provided a 


Brown of the 


was dissolved at the end of 


constant help and inspiration. 
[his was a most interesting period in my life. 
[he Good Book says, “No man can serve two 


masters for . . . he will hold to one and despise the 
other.” I had three—an Army boss, Lt. Col. 
later Maj. Gen.) Conger Pratt, Commanding 


Mitchel Field 
Capt. Land of the Navy; and Mr. Guggenheim, a 
civilian. They were all such fine, understanding, 
and cooperative people that my existence, instead of 
being complex, was uncomplicated and extremely 
pleasant. 

Our first activity in the Full Flight Laboratory 


Officer of where I was stationed; 


The second Full Flight Laboratory aircraft 
was a Navy-type Vought Corsair O2U-1 
“to be used for cross-country practice fly- 
ing.” (IAS Photo, from the Sherman Fair- 
child Collection of Aeronautical Photo- 
graphs) 
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was to study and endeavor to analyze the work pre- 
viously done on blind landing in fog. 

In England tethered balloons had been lined up 
with the landing field and used with some success 
in still air and when the fog was not thick. This 
concept we abandoned at once as unsatisfactory, for 
experience had indicated that the fog layer might be 
very thick and still air could not be depended upon 
at all times when visibility was restricted—for ex- 
ample, in a blizzard. 

In both England and France the lead-in cable idea 
was tried. In this system an electrified cable circled 
the field and led into a landing. It required very 
sensitive sensing equipment in the airplane, and it 
was necessary to make a precision turn into the field 
at low altitude. This turn presented considerable 
difficulty. 

Lieutenant LeRoy Wolf of the Army Air Corps 
also experimented with the electrified cable concept 
at Wright Field. 

The U.S. Navy had some success, with an elec- 
tromagnetic cable guide at the Lakehurst Naval Air 
Station. 

The low-frequency radio range had been developed 
by the Bureau of Standards and the Army and was in 
limited use for aerial navigation. An adaptation of 
this radio range in the form of a radio homing beacon 
seemed to offer the greatest promise for our use. It 
could also be readily tied in with the radio receiver 
and other conventional airborne equipment. 

Actual blind landings had been attempted with 
dragging weights and long tail skids. These either 
gave an indication upon touching the ground or were 
rigged to actuate the aircraft control. 

There was additional background information but 
the material outlined here was fairly representative. 

The first important expenditures made by the Full 
Flight Laboratory were for two modern airplanes. 
One, a Consolidated NY-2 military training plane 
mounting a J-5 engine, was to be used in the instru- 
ment landing experiments and to test instruments, 
equipment, or devices that might be helpful in over- 
coming fog flying problems. It had the large wings 
used by the Navy for pontoon seaplane training, but 
mounted, in place of the pontoons, a specially rein- 
forced landing gear with long oleo action. As a 
training plane it had a very high factor of safety, 
was extremely rugged, and was inherently stable 
about all three axes. (Special flying and landing 
wires permitted the rigging in of additional dihedral.) 

The acceptance test flight on this airplane was 
made November 3, 1928, at the factory near Buffalo, 
and the airplane was flown to Mitchel Field the next 
day. It was flown considerably in November and 
December and then delivered to the Radio Frequency 
Laboratories in Boonton, N.J., to have voice radio 
installed. It was at Boonton for 4 weeks this first 
time and frequently thereafter. 

The second airplane was a Navy-type Vought 
Corsair 02U-1. It mounted a Pratt & Whitney 
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Wasp engine. It was to be used for cross-country 
practice flying and was an excellent airplane for the 
purpose. It was a fast, good-flying airplane, but 
not so rugged and stable as the NY-2. The 02U-1 
was delivered to the Fund on November 21, and the 
to Boston, with Harry 


was made on the Saine 


first cross-country flight 
Guggenheim as passenger 
day 

The Army Air Corps made a hangar available at 
Mitchel Field for the use of the Full Flight Labora- 
tory. It was initially provided for the safe airplane 
competition, another important contribution of the 
Fund. The Army also provided the full-time serv- 
ices of an excellent mechanic to maintain the 
Laboratory’s aircraft. This was Cpl. (later Sgt.) 
Jack Dalton. The continued excellent mechanical 
condition of the two aircraft thereafter was largely 
the result of his competence and devotion. 

Lieutenant (later Brig. Gen.) B. S. Kelsey was 
made available by the Army Air Corps as flight assist- 
ant and safety pilot. When flights were made 
under the hood it was necessary, in the interest of 
safety, to have another pilot in the airplane to look 
out for other aircraft and also to make sure that the 
pilot under the hood did not get into difficulty be- 
cause of instrument or equipment failure. Ben, from 
the start, was a full-fledged member of the team and 
did much of the experimental work. His piloting 
help, criticism of tests carried out, sound technical 
counsel, and ever-pleasant personality contributed 
greatly to the results achieved. 

As the preliminary practice flights progressed, it 
soon became apparent that, even with the very stable 
and sturdy NY-2, the available instruments were 
not adequate. For determining heading when 
maneuvering and when landing, the compass, be- 
cause of the northerly turning error, was entirely 
unsatisfactory and the bank and turn indicator, 
though excellent for its purpose, was more a quali- 
tative than a quantitative measuring instrument. 
Also, at the moment of touch-down in a blind land- 
ing, it was desirable that the wings be level with the 
ground. This was not easy to assure, particularly 
when the wind was gusty. 

An accurate, reliable, and easy-to-read instru- 
ment showing exact direction of heading and precise 
attitude of the aircraft was required, particularly 
for the initial and final stages of blind landings. 

Two German artificial horizon instruments, the 
Anschutz and the Gyrorector, were studied but were 
not deemed entirely satisfactory. 

I sketched a rough picture of the dial for an in- 
strument which I thought would do the job and 
showed it to Elmer Sperry, Sr., a great engineer and 
inventor who had established and headed the Sperry 
Gyroscope Company and who was very interested 
in aviation. It was, in substance, the face of a 
directional gyro superimposed on an artificial hori- 
zon. He advised that a single gyroscopic instru- 
ment could be designed to meet the requirements 
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but recommended, for simplicity of construction, two 
separate instruments. I agreed and he then as- 
signed his ingenious son, Elmer, Jr., to work with us 
and to be responsible for the design and fabrication 
of the two instruments. We could not have had a 
better colleague. Elmer also soon became a member 
of the team and spent as much time at the hangar 
and at the evening ‘‘discussion’’ sessions in our 
quarters on Mitchel Field as the rest of us. Out of 
this, as you know, came the Sperry artificial horizon 
and the Sperry directional gyroscope which still, 
with their descendants, are on the instrument panel 
of every airliner and military airplane today. 

As time progressed, literally hundreds of blind and 
simulated blind landings were made. To make a 
landing, the airplane was put in a glide at 60 mph, 
with some power on, and flown directly into the 
ground. Although this was about 15 mph above 
stalling speed, the landing gear absorbed the shock of 
landing and, if the angle of glide was just right, the 
airplane didn’t even bounce. Actually, after a while, 
it was possible to make consistently perfect landings 
by this method. To assure just the right amount of 
power in the glide, a mark was placed at the proper 
place on the throttle quadrant. 

Excellent cooperation was obtained from every 
organization and individual with whom we worked. 
Among them were the Pioneer Instrument Company, 
the Taylor Instrument Company, the Radio Cor- 
poration of America, and the Bell Telephone Labora- 
tories, who installed the modern radio transmitter 
and provided miniature earphones with molded 
plugs. Along with the other BTL _ personnel, 
Capt. Russell Luff Meredith and Capt. Ray Brooks, 
recently retired from the Army Air Corps, were very 
helpful, as was the Aircraft Radio Corporation of 
Radio Frequency Laboratories which installed the 
excellent radio receiver. Dr. L. M. Hull, the Presi- 
dent, Dr. A. W. Parkes, Dr. F. H. Drake, and all 
members of the staff were most cooperative. Plus 
the Kollsman Instrument Company, the Sperry 
Gyroscope Company, and many other industrial 
concerns, scientists, engineers, inventors who helped 
and just well-wishers. All helped. 

Valuable support was received from the Bureau of 
Standards which designed and installed most of the 
ground and airborne radio navigation equipment. 
Dr. Harry Diamond and his associates of the Bureau 
spent much time with us and could not have been 
more cooperative. 

While aural signals were satisfactory for rough 
aerial navigation, it soon became apparent that a 
visual indicator would be much better for the precise 
directional control required in blind landings. 

To meet this requirement, the Bureau of Stand- 
ards, working with the Airways and Radio Division 
of the Department of Commerce, designed a 2-kw, 
semiportable, two-leg range, which was used as a 
homing beacon, and a fan-type marker beacon. 
The homing range was installed on the west side 


of Mitchel Field. The marker beacon sat on the 
leg of the homing range and was located on the east 
side of the field. 

The indicator for the former, carried in the air- 
plane, was a pair of juxtaposed vibrating reeds. If 
to the right of the course, the right reed vibrated 
through the greater amplitude. If to the left of the 
course, the left reed vibrated more vigorously. If 
on course, both reeds vibrated through the same arc. 
As the station was approached, the amplitude of 
vibration increased. If the amplitude of vibration 
became too great, it could be reduced through use 
of a rheostat. 

A single reed started to vibrate as the fan-type 
marker beacon was approached. It reached maxi- 
mum amplitude, then quickly dropped to zero when 
the airplane was directly overhead, rapidly built up 
to maximum again, and then tapered down. The 
homing range indicator also had a distinct null 
when the airplane was directly over the range station. 

As the tests progressed, the instrumentation and 
equipment were constantly improved until toward 
the end of 1929, during the final stages of the flight 
tests, the following instruments and equipment 
were carried : 


(1) Normal engine instruments 
(2) Magnetic compass 
(3) Earth inductor compass 
) Bank and turn indicator 
(5) Directional gyro 
) Artificial horizon 
(7) Air-speed indicator 
(8) Altimeter 
(9) Rate of climb indicator 
(10) Outside air thermometer 
(11) Vibrating reed homing range indicator 
(12) Vibrating reed marker beacon indicator 


Considerable thought was given to the location or 
arrangement of each instrument in order to facilitate 
reading and reduce pilot fatigue. Fatigue led to 
errors, and piloting errors could not be tolerated in 
instrument landings. 

Small instruments were preferred to the more 
conventional, larger ones because, even though some- 
what harder to read individually, the small instru- 
ments permitted more compact and logical arrange- 
ment and were easier to read and interpret en masse. 
It was soon determined that small instruments could 
be made easier to read by the use of broader hands 
with white or radium paint applied to the outer half 
of their length. 

It was found that the instruments could be read 
more quickly, and over a longer time period with 
less fatigue, if the arrangement of instruments and 
the position and direction of motion of the indicating 
hands was ‘“‘natural.’’ Also, any abnormality or 
improper indication should be detectable automati- 
cally by a quick glance at the instrument panel. This 
did not ‘solve’ instrument flight but did simplify it. 

The airplane, in addition to its flight instruments, 
carried a RFL radio receiver, a BTL radio trans- 
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mitter, two 6-inch Pyle national landing lights, and a 
parachute flare. 

A larger-than-customary Leece-Neville generator 
was installed on the engine to assure an adequate 
and continuous electrical supply. 

The mast-type receiving antenna, which was 
employed to minimize directional effect, required 
considerable development before a tendency to 
vibrate under certain flight conditions could be 
corrected. 

A trailing wire antenna was used for transmission 
in normal flight. This was reeled in, and a fixed 
wire antenna was used for transmitting when land- 
ing. 

The specialized ground equipment consisted of a 
radio receiver and transmitter which provided voice 
communication with the aircraft and a Kollsman 
sensitive altimeter with which the altimeter in the 
airplane was synchronized by radio. 

In addition, there were the visual homing range 
and the visual fan-type marker beacon previously 
mentioned, and, of course, the standard Mitchel 
Field Army-type aural beacon to lead an aircraft to 
the vicinity of the airfield. 

Two problems that were available very much on the 
minds of the Fund directors were collision of aircraft 
in the air and formation of ice on the wings, structure, 
propeller, and in the carburetor of aircraft, when at- 
mospheric conditions were conducive to icing. 

Dr. W. C. Greer and Dr. Merrit Scott at Cornell 
University carried out work, under Fund sponsor- 
ship, on the underlying causes of ice formation on 
aircraft. 

The Fund also sponsored the work done by John 
P. Kilgore of New Haven, Conn., on an electrically 
heated wing blanket for deicing. 

The Full Flight Laboratory did no development 
work on deicing equipment, but many test flights 
were made under icing conditions to determine 
effects and limitations. 

The NY-2 was frequently out of commission while 
new instruments or equipment were being installed. 
These were convenient periods for cross-country 
flying practice under unfavorable weather conditions 
in the O2U-1. This airplane had all necessary flight 
instruments but no blind landing equipment. 

An extreme example of a cross-country bad 
weather flight took place on March 15, 1929. I 
took off from Buffalo in the O2U-1 headed for 
Mitchel Field. It was night and the weather was 
fair and improving at Buffalo but marginal to the 
south and east. This was to be a difficult flight 
but possible and just the sort of thing required to 
establish flight ‘‘limitations.”’ In a pinch I could 
return to Buffalo any time up to the point where 
nearly half of the gasoline supply was used up. 

I well realized that the pilot who flew within his 
limitations would probably live to a ripe old age, 
whereas the pilot who flew: beyond them would not. 
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I also knew that different pilots had different limita- 
tions. 

This was pointed up in the mid-twenties when | 
was a test pilot at old McCook Field. At that time 
there were few facilities and little ground equipment 
to do environmental testing on new airborne de- 
vices. It was therefore necessary to test them out 
in flight, and the test pilots spent many hours flying 
around the airfield to see how a new device held up 
under the accelerations, vibrations, and changes in 
temperature and pressure experienced in flight. 

Lieutenant Alex Pearson always spent these hours 
for instance, holding 
As a result he became 


practicing precision flying 
constant speed and altitude. 
extremely proficient and could fly a better speed 
course or do a smoother saw-tooth climb than any of 
the rest of us. 

I spent the hours flying low in the vicinity of 
McCook Field and on the main air routes in and out 
I knew every high build- 
ing, tree, silo, windmill, radio tower, and high- 


memorizing the terrain. 
tension line in the area. I could therefore fly in—or 
adverse weather safely when other equally 
experienced pilots didn’t fly. 


under 
This was not because 
I was a better or more daring pilot than my col- 
Constant practice had just extended my 
limitations. 


leagues. 


The trick was to learn your limitations, gradually 
expand them, but never go beyond them. 

I thought I was being pretty sharp but the Com- 
manding Officer, learning that I frequently flew in 
that area when other pilots did not, thought differ- 
ently. Unaware of my training plan, he removed me 
from the job of chief pilot in the flying section, 
advised me that I didn’t have judgment enough to 
be a pilot, and assigned me to the airplane section 
as an aeronautical engineer. 

All of these things went through my mind as the 
I planned to fly contact all 
the way so, in order to avoid the mountains, took 
Buffalo, Syracuse, Utica, 
Schenectady, Albany, and then down the Hudson 


weather deteriorated. 


the route Rochester, 
River. 

There was no particular problem getting to Al- 
bany, but from there on the ceiling and visibility 
ccame marginal. 

Soon I had passed “‘the point of no return” and no 
longer had gasoline enough to go back to Buffalo. 

At one place I found it expedient to slow down and 
hover with the left wing of the airplane over a 
brightly lit southbound passenger train traveling 
along the east side of the Hudson River. Presently 
it went through a cut and wiped me off. This was 
too hazardous, so I left the train behind and fol- 
lowed the river bank. Considered crossing the 
river and landing in the parade ground at West 
Point, but abandoned this idea as the weather re- 
mained flyable—barely. 


Upon reaching the lights and heat of the City of 


New York, the ceiling and visibility improved 
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slightly. Flew south to the Battery hoping to be 
able to get to Mitchel Field from there, but the East 
River and the area to the south were ‘‘socked in’’ so 
I could not go on. Tried to get to Governor’s 
Island and land on the drill ground, but it was fog- 
shrouded. 

Turned north back up the Hudson intending to 
land on the Yonkers Golf Course which I knew well, 
but the fog was below the river bank so had to 
abandon that idea. 

Returned to the Battery with the intention of 
crash-landing in Battery Park. A chap ran out into 
the middle of the Park and waved me off. He 
apparently thought I mistook it for a flying field. 

It is interesting to note that the George Washing- 
ton Bridge across the Hudson at 179th Street was 
under construction at this time. There were as yet 
no suspension cables or other horizontal structure. 
Only the great vertical piers on each side of the River 
had been completed. I had passed the east pier 
three times without seeing it. 

About this time it appeared that a crash-landing 
in the River might be necessary, so I removed my 
parachute in order to be able to swim ashore. 

The water, on closer inspection, looked uninviting 
so decided on a final try—this time for Newark 
Airport—and headed across the Hudson. As soon 
as the River was crossed and the lights south of 
Jersey City appeared, it became obvious that this 
last chance was impractical. Thereupon pulled up 
into the fog and climbed through it. It was only 
about a thousand feet thick and crystal clear above. 
Intended to fly west until past the highly populated 
part of the metropolitan area and then jump. The 
gasoline gage had been fluttering on zero for some 
time. Noted about this time that my parachute 
harness was off and promptly put it on. 

About over Kenilworth, beyond Elizabeth, saw a 
revolving beacon through a hole in the fog and a 
flat-looking area adjacent to it with no lights. 
Hoping it might be an emergency field or at least an 
open area (although realizing that it might be a 
wood or lake), turned the landing lights on and 
dived through the hole and scouted the area. 
The bottom of the fog was still very low, and I tore 
the left lower wing badly on a tree top. The air- 
plane still flew, although almost completely out of 
gasoline, so I returned to the most likely spot and 
crash-landed. Took the impact by wrapping the 
left wing around a tree trunk near the ground. 
The O2U-1 was completely washed out—quite 
beyond repair—but I was not even scratched or 
bruised. 

The moral of the story is that with the NY-2 
mounting blind landing equipment and the Full 
Flight Laboratory radio station alerted at Mitchel, 
this would have been a routine cross-country flight 
with ‘‘no sweat.”’ 

The flight pointed up the importance of constant 
radio communication between the aircraft and the 


ground and the need for frequent and accurate 
weather reports obtainable by radio during flight in 
order to assure safe continuation or indicate suit- 
able alternate destinations. It also indicated the 
desirability of a special light to mark emergency 
fields for night landings. Later green lights were 
used. 


In general, the weather a pilot could fly in safely 
was a function of airplane characteristics, ground 
facilities and procedures, ground and airborne equip- 
ment and instrumentation, general pilot skill, and his 
specialized knowledge of the local aids to air naviga- 
tion, the terrain, and the weather conditions to be 
expected in the area. 


An important requirement in instrument landing 
was to have a precise measure of altitude when 
approaching the ground for a landing. 

The conventional barometric altimeters of the 
day measured, at best, to the nearest 50 or 100 feet. 
It would be handy if an altimeter were available that, 
near the ground, would measure to 10 or even 5 
feet. 


The Kollsman Instrument Company developed 
such an instrument and on August 30, 1929, I was 
very pleased to take Paul Kollsman and his new 
instrument up on its first test flight (in the second 
O2U—a more modern O2U-2 purchased after the 
crash of the O2U-1). Mr. Kollsman held the sensi- 
tive altimeter in his lap during the flight and it 
performed perfectly. Here was another important 
addition to flight instrumentation. We promptly 
obtained one and installed it in the NY-2. 


This instrument had two hands and a multiplica- 
tion factor of twenty between them. Actually the 
fast-moving hand made one complete revolution for 
each 1,000-foot change in altitude, which meant a 
movement of about 5/3. inch for 20 feet altitude. 
This was over an order of magnitude more accurate 
than earlier altimeters. 

Although the Kollsman altimeter provided a 
considerable advance in instrumentation, it still 
measured the barometric altitude or height above 
sea level. An instrument which would measure 
the height above the ground regardless of changes 
in barometric pressure would be of great value. 


With this idea in mind, several companies were 
encouraged to work on a sonic altimeter. One 
built by Sperry was tested out in the NY-2. A note 
on a frequency of about 950 cycles was directed at 
the ground from a megaphone on the bottom of the 
airplane and picked up by a detector on the airplane 
after having been reflected back from the ground. 
The elapsed time interval was measured, as in a 
marine “‘fathometer,’ and the altitude above the 
ground thus determined. The concept was theo- 
retically sound, but in its initial form the equipment 
caused considerable drag and was unduly large, 
heavy, difficult to install, and complicated. It also 
appeared that a radio altimeter measuring the phase 
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difference of radio waves reflected back from the 
ground offered more promise. 

Several radio altimeters were under development 
with Fund encouragement—so the sonic altimeter 
experiments were abandoned. 

Although major emphasis was given in the Full 
Flight Laboratory to flight tests designed to permit 
flying safely and reliably despite fog, the Fund, 
and the Laboratory, were also interested in experi- 
ments on the penetration of fog by light rays and on 
the dispersal of fog. 

Studies on fog penetration were made by Dr. 
Anderson of the University of Washington and Dr. 
Barnes of Bryn Mawr College under Fund sponsor- 
ship, and by Dr. Julius A. Stratton of M.I.T. for 
Col. E.H.L.Green. The visual, infrared, and ultra- 
violet parts of the spectrum were carefully explored. 

It was concluded that, though there was some 
difference in penetration due to wave length, there 
was little likelihood of finding any visual light which 
would penetrate thick, dense fog. Dr. Stratton’s 
further investigations gave some promise that at 
very short wavelengths—below one centimeter— 
there was a possibility that fog might be penetrated 
by dispersal. 

In addition to the scientific studies carried out, 
various inventors presented ideas for fog penetra- 
tion, all of which received careful consideration. 

Four basic methods of fog dispersal were originally 
considered : 


Mechanical 

Here the concept was to have a large propeller 
or series of propellers “churn up”’ the fog. 
Chemical 

Experiments carried out by Henry G. Houghton, 
Jr., of M.1.T., using hygroscopic materials, showed 
some promise both in the Laboratory and in full- 
scale tests which were conducted on the estate of 
Colonel Green in South Dartmouth, Mass. 


Electrical 

From the early 1920's Dr. Warren of Hartford, 
in collaboration with the Army Air Service, had 
been experimenting with cloud and fog dispersal by 
means of electrified sand particles dropped from an 
aircraft. He was occasionally successful in dis- 
persing small clouds. 

Mr. Flowers experimented with fog dispersal 
through the use of electrified water particles and 
achieved some interesting results. 


Heat 

Experimentation on, and the actual accomplish- 
ment of, fog dispersal by heat has been carried on 
until fairly recent times. 

FIDO—Fog, Intense Dispersal of. FIDO was 
effectively used in World War II on three of the 
long emergency and ten of the main landing fields 
in England. It consisted of powerful heat sources 
at intervals along the runways which heated the air 
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to a temperature above the dew point, and the 
nearby fog was thereby dispersed. 

Approximately 2,500 and fighters 
mostly RAF but including some USAF—returning 
from missions and finding their home base, and much 
of England, covered in pea soup fog were directed 
to the cleared ‘‘tunnels in the fog’’ over these FIDO 
fields and safely landed. 

Perhaps I should mention here that the Eighth 
Air Force Bomber units in England were being 
without 


bombers 


trained and equipped to operate “‘blind”’ 
FIDO—and at the end of the war two groups could 
take off and land regardless of fog. In short time, 
the entire force would have been so trained. 

In the late forties successful FIDO experiments 
were carried out at the Naval Air Base at Arcata 
on the coast of northern California. Arcata was 
chosen because it was in one of the foggiest areas on 
the entire Pacific Coast. Southwest Airways suc- 
cessfully used the Navy fog dispersal system from 
November 1947 until November 1949. 

In 1953, a FIDO installation was made at the 
Los Angeles International Airport. The cost for 
installation, modification, and test was $1,325,000. 
While the installation appeared to handle moderate 
fogs fairly well, it did not satisfactorily disperse 
dense fogs where the visibility was one eighth of a 
The dispersal difficulties seemed to 
arise from the great amount of heat required, the 
airport configuration, and the movement and in- 


mile or less. 


trusion of the fog. Further operational problems 
resulted from infrequent use of equipment and the 
high cost of installation, maintenance, and operation. 
It was also anticipated that there would be difficulty 
experienced in moving passengers to and from the 
airport and in moving aircraft to and from the cleared 
runway when the fog was really dense. The fol- 
lowing year the FIDO program was abandoned at 
Los Angeles. 

It appears that FIDO would be a successful way 
of coping with intense fog only when the air is com- 
paratively still. When there is any considerable 
movement of the air particularly if the 
movement has a component across the runway, the 
cleared air mass moves on and new fog comes in 
faster than it can be dispersed. Difficulties are also 
experienced in an absolutely dead calm because of 
conduction, particularly at runway intersections and 


and fog 


elsewhere where there are no burners. 

The Full Flight Laboratory’s first experience with 
FIDO 

Mr. Reader of Cleveland operated a gravel pit 
and utilized a large blowtorch type of heater to dry 
the gravel and sand. He observed that if there was 
an intense fog when he turned the heater on, the 
fog in the immediate area dispersed. 

Learning of the Guggenheim fog flying experi- 
ments, he wrote giving information of his expert 


of sorts—occurred in 1929. 


ence and advising that he was interested in helping 
solve the fog problem. 
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He was invited to bring one of his heaters to 
Mitchel Field where it was installed just east of the 
last hangar. 

For several months thereafter we waited in vain 
for a dense fog. Finally, on September 24, 1929, it 
Someone, I think it was Jack Dalton, 
awakened just before daylight and noted that there 
was a zero-zero fog covering the area. 
quickly called together. 

We immediately notified Mr. Reader who arrived 
shortly thereafter. Mr. Guggenheim was also noti- 
fed but he had to come from Port Washington 
and didn't arrive until later. 

The equipment was fired up, but the fog did not 
disperse except in the immediate vicinity of the blow 
torch. The experiment was a disappointing failure. 
At the time we considered the concept impractical 
and the equipment was removed. 

In retrospect, it appears that the trouble was prob- 
ably fog movement or possibly the absence of a 
mineral mass to store and reflect the heat. Also, 
for satisfactory fog dispersal over a considerable 
area, about two orders of magnitude more heat 
energy would seem to be required. 

Though we were all disappointed we were there, 
and the fog was there, so I decided to make a real 
fog flight. 

The NY-2 was pushed out of the hangar and 
warmed up. The ground radios were manned and 
the radio beacons were turned on. I taxied out to 
the middle of the field and took off. Came through 
the fog at about 500 feet and made a wide swing 
coming around into landing position. By the time 
I had landed 10 minutes after take-off, the fog had 
just started to lift. 

At about this time Mr. Guggenheim, along with 
several other people, arrived so we decided to do an 
‘official,’ under the hood, flight. I'd just made a 
real flight in the fog so wanted to go alone, but Mr. 
Guggenheim insisted that Ben Kelsey be taken along 
as safety pilot. The fog had lifted considerably by 
this time and he was afraid there might be other 
aircraft in the vicinity. 

We both got into the airplane and the hood over 
my cockpit was closed. 

The engine was again warmed up and I taxied the 
airplane out and turned into the take-off direction 
on the radio beam. Took off and flew west in a 
gradual climb. At about 1,000 feet the airplane 
was leveled off and a 180° turn was made to the left. 
This course was flown several miles and another 
180° to the left was made. The airplane was lined 
up on the leg of the radio range located on the west 
side of Mitchel Field and a gradual descent started. 
Leveled off at 200 feet above the ground and flew 
at this altitude until the fan beacon on the east; side 
of the airfield was passed. From this point the air- 
plane was flown into the ground, using the instru- 
ment landing procedure previously developed. Actu- 
ally, despite previous practice, the final approach 


came. 


Our gang was 


and landing were sloppy. This entire flight was 
made under the hood in a completely covered cock- 
pit which had been carefully sealed to keep out all 
light. The flight, from take-off to landing, lasted 
15 minutes. It was the first time an airplane had 
been taken off, flown over a set course, and landed 
by instruments alone. This was just 10 months and 
3 weeks from the first test flight of the NY-2. 


The tests carried out thus far indicated, among 
other things, that before all-weather flying would be 
practical there was a need for 


(1) Better coordination between the long-range 
aural beam and the short-range visual landing beam. 


(2) Accurate measurement of distance along the 
landing beam. 


(3) A slant glide beam at right angle to the 
vertical landing beam which might be curved to 
become tangent to the earth rather than intersect 
it in a straight line. 

(4) The possibility of automatic volume control 
on the receiving radio. 

(5) A good gyro pilot to assist the pilot in flying 
the aircraft. This could be, and later was, tied in 
with the instrument landing systems. 


(6) Much more work on ignition shielding not 
only to reduce noise but to permit flight in heavy 
rain. 

At the end of 1929, the Fund was disbanded. , It 
was felt that it had made the necessary initial con- 
tribution and that further development could better 
be carried out by other organizations. The Full 
Flight Laboratory went out of existence with the 
Fund, and the NY-2 was turned over to the U.S. 
Army Air Corps and moved to Wright Field where 
instrument landing experimentation and develop- 
ment were continued under the direction of Ist Lt. 
(later Maj. Gen.) Albert Hegenberger. 

The Bureau of Standards and the Bureau of 
Air Commerce of the Department of Commerce 
also continued their work and on March 20, 1933, 
James L. Kinney, with William La Violette as 
mechanic and Harry Diamond as passenger, flew a 
Bellanca airplane by instruments alone, from College 
Park, Md., to Newark Airport, a distance of about 
200 miles. He landed blind using the recently de- 
veloped bent landing beam. The weather was actu- 
ally below minimum. 

The Collier Trophy is awarded each year to an 
individual or group ‘“‘for the greatest achievement in 
aviation in America, with respect to improving the 
performance, efficiency, or safety of aircraft, the 
value of which has been thoroughly demonstrated 
by actual use during the preceding year.” 

In 1934 the Collier Trophy was awarded to Capt. 
Albert Hegenberger, U.S. Army Air Corps, for the 
development and demonstration of a_ successful 
blind landing system. 


On May 9, 1932, Capt. Hegenberger made the 


October 1961 + Aerospace Engineering 


67 


l 
1 
) 
g 
d 
a 
1S 
n 
m 
iT 
te 
se 
a 
to 
1e 
n- 
ns 
he 
yn. 
ty 
he 
ed 
ol- 
at 
ay 
m- 
le | 
he 
he 
inl 
lso 
of 
ind 
ith 
pit 
dry 
was 
the 
eri- 
yeri- 
ing 


first solo blind flight, alone in the plane, depending 
solely upon instruments from take-off to landing. 

Over a period of years he and his colleagues had 
greatly improved the equipment previously used 
and developed additional equipment. He devised 
a blind landing procedure which became standard 
military practice. This system was put into actual 
use in 1934. 

Flying by instruments had outgrown the early 
experimental phase. It was a practical reality, and 
aviation had entered a new era. 


High Stress Level Fatigue 


(Continued from page 19) 


Experimental Basis of the Theory 


Generally speaking, the fatigue damage from 
which the crack is formed is localized near the 
surface of the specimen.' The cracks start in broad 
slip bands, traces of which persist during electro- 
polishing treatments, and the fatigue cracks are 
found to start quite early in the life of the specimen. ! 
Grooves and ridges are formed on the surface of 
the specimen along the slip bands, and the ridges 
sometimes develop into thin tongues of metal.’ 
The mechanism of the formation of these tongues of 
metal and the holes, more familiarly called extrusions 
and intrusions, seems to be independent of tempera- 
ture in the sense that the phenomenon was observed 
at liquid air temperatures. So it appears that, 
while a thermally activated mechanism may play 
a part in fatigue failure, it is not an essential feature 
of fatigue. However, the thermally activated 
mechanisms can definitely compliment, enhance, 
and accelerate the deterioration characteristics in 
fatigue. The crack formation generally observed 
seems to occur in the intrusions on the surface. 
Indications are that this is a purely geometric 
process* in the sense that it may be understood as 
a consequence of the motion of dislocations in the 
several intersecting crystallographic planes. 

It is appropriate at this stage to introduce the 
findings of Wood.‘ He shows, on the basis of fairly 
conclusive evidence, that two essentially distinct 
types of failure have been hitherto classified as 
fatigue. One predominates at large amplitudes of 
stress and the other at small stress amplitudes 
where the curve of stress vs. number of cycles in 
fatigue tends to turn parallel to the cycle axis. 
These two ranges are termed by Wood the H and 
F ranges, respectively. The mode of deformation 
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Symbols 


A = area of plastic enclave ahead of crack 
A = area of cross section of test specimen 
= X Lm)/b, an experimental constant to 
evaluate Eq. (19) 
Ct. 
Cp) = constants 
fi = Young’s modulus 
Es = secant modulus 
E = 
K = a constant representing width of hysteresis 
loop of enclave material after it strain- 
hardens; an experimental constant to 
evaluate Eq. (19) 
Ko, 
K;, Ky = constants 
Ky = theoretical stress concentration factor in 
vicinity of a notch 
Ke =a plane stress factor associated with A in 
elastoplastic solution for area of plastic 
enclave 
= grain size 
= mean grain size 
N, = actual number of stress pulses after which 
failure may be expected in completely 
random loading consisting of Nr stress 
pulses 
N; = total number of cycles to failure in fatigue 
at a single stress o,, Eq. (21) 
N, = total number of cycles to failure at stress 
level o, 
Ng -1 = total number of stress pulses in random se- 
quential loading which can be applied be- 
fore failure occurs; subscript (S — 1) 
indicates block number 
Ny = total number of cycles assumed to be ap- 


plied during anticipated design life of a 
test specimen in random loading 

N.), = equivalent number of cycles at reference 
stress corresponding to application of Nr 
cycles in completely random loading 

N,)z = total number of cycles to failure at reference 
stress level corresponding to maximum 
stress in either program type, random, or 
acoustic loading 


Nr; = total number of cycles to failure in program 
type loading = n, X An rp 
No = cycles associated with failure at ultimate 


tensile strength, Np = 1/4 
= fatigue life at stress levels a; and o, re- 


> 
= 
| 


spectively, in discussion of two-step cumu- 
lative damage 
NV." = N2 


number of cycles at stress level o;, beyond 


which idea of cumulative damage has no 
meaning in two-step cumulative damage 

Nie = total number of cycles to failure in cumula- 
tive damage when stress g; is applied first 
and then ge is applied, 0; < o 

Noy = total number of cycles when order of ap- 
plication of stress o; and go» is reversed, 
< G2 

AN = sum of cycles of applied stress for dormant 
period and active incremental growth of 
crack AN = AN, + AN, 


in the H or hardening mechanism range is sub- 
stantially the same as in static deformation. The 
process of deformation is essentially one of coarse 
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dormant period in incremental crack growtl 

active period in incremental crack growth 

mean number of cycles at which failure may 
be expected in random loading consisting 
of possible application of Np cycles 

maximum number of cycles after which fail- 
ure may be expected 

equivalent number of cycles at reference 
stress corresponding to Nmean 

residual static strength of test specimen 
after subjecting it to fatigue for m cycles, 
Eq. (23) 

ratio of minimum stress amplitude oa,’ to 
maximum stress amplitude o, in tension 

surface energy density in Griffith theory of 
cracks 

sum of fractional lives in two-step cumula- 
tive damage when order of stress applica- 
tion is 01, (6; < 

same as Sj. when order of application of 
stress is o2, < a2) 

initial area of hysteresis loop of plastic en- 
clave ahead of crack before it hardens 

final area of hysteresis loop after enclave 


hardens 
width of plastic enclave; Wp ~ A‘/? 
Burger’s vector, a characteristic lattice 


spacing 


= crack length at any stage 


Il 


length of a Frank-Read source 

length of a crack which is just critical for 
given applied tensile load 

mean length of crack when variations of 
grain size are taken into account 

crack length at the end of No cycles of ap- 
plied stress, assumed here to be equal to 
Griffith crack length for failure at o,,;, 

crack lengths in discussion of cumulative 
damage after application of cycles of 
stress a; and respectively 

incremental length of crack when original 
crack length is (J + Al) 

arbitrary number of cycles at any stress 
level 

number of dormant periods in fatigue, mp = 
N/AN, 

number of blocks to failure in program type 
of loading 

number of cycles of stress applied at stress 
oq in cumulative damage 

distribution function of grain size 

distribution function of Frank-Read source 
lengths 

equivalent number of cycles at stress levels 
o; OF o2 in two-step cumulative damage 
when n cycles are applied first at o2 or oi, 
respectively 

number of cycles at stress o; and o2 applied 
in cumulative damage in familiar cumu- 
lative damage notation 

total number of cycles in a block in program 


type of loading Anrp = pay An, 
q=1 
number of cycles at o, in program loading 
(An,) 
r/q 
q=1 


An 


erg 


Fult 


To 


Toct 


= equivalent number of cycles at reference 


| 


= Rayleigh 


= tensile stress appropriate for 7, 


= reference 


stress level o, corresponding to Ang cycles 
at a stress level o, 


= total equivalent number of cycles at refer- 


ence stress level corresponding to stress 
cycles at different frequencies in acoustic 
fatigue 


= actual time to failure of panel in acoustic 


fatigue 


= time to failure at reference stress level and 


reference frequency 
q 
erg 


q=1 


= equivalent time at reference stress level and 


frequency corresponding to Amerg cycles 


functional constants in discussion of cumu- 
lative damage in two-step loading 


= number of dislocations piled up in a grain 


due to operation of a Frank-Read source 


= shear modulus 


= normalization factors of distribution func- 


tions used in discussion of fatigue under 
random loading 

maximum tensile stress at leading-edge dis- 
location of a pile-up of dislocations 


= internal tensile stress in direction of applied 


stress for which 7; is appropriate shear 
stress; also interpreted as endurance limit 
or natural elastic limit 

peak amplitude stress for rms 
stress at w, 


= an arbitrary stress and its range in discussion 


of multiple-step and random loading 


= rms stress peak in acoustic fatigue at fre- 


quency w, 


oct Na tension 


test 


| 


= maximum applied tensile stress in a fatigue 


cycle for which 7, is appropriate resolved 
shear stress 


= minimum stress amplitude in a fatigue cycle 


= stresses in discussion of two-step cumula- 


tive damage 


= stress far away from local discontinuities 

= residual static strength in psi; o,, = P,,/A, 
= ultimate tensile strength 

= endurance limit in fatigue when test stress 


ratio is R 


= mean stress in acoustic fatigue at frequency 


Wg 


= endurance limit in fatigue when mean stress 


1S 


= internal shear stress that opposes dislocation 


motion 


= resolved shear stress in slip plane in slip 


direction when a tensile stress o, is applied 


= limiting octahedral stress in elastoplastic 


theory of cracks 


= frequency of a mode of vibration in acoustic 


fatigue 


frequency in acoustic fatigue 


analysis 


slip and appears to occur primarily through the 
operation of 


Wood calls the fatigue failure in the high stress 


Frank-Read sources.5 Accordingly, 


range a “delayed static fracture.” 


| 


It is delayed 


because the hardening process, which appears to 
precede the nucleation and propagation of cracks, 
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occurs in a less efficient manner in this range of 
fatigue stress than it does in normal static stressing. 

On the other hand, in the F range, conditions are 
significantly different. The deformation is primar- 
ily due to fine slip and involves little or no harden- 
ing. Failure of the material appears to be due to 
a gradual deterioration of structure which manifests 
itself as the “‘persistent slip bands,” and the fissures 
that lead to eventual failure seem to start in these 
regions. This process is presumably the result 
of the to-and-fro motion of a large number of 
individual dislocations rather than the avalanche 
release of dislocations characteristic of the Frank- 
Read sources that are apparently responsible for 
coarse slip. Wood calls the failure in this region 
pure fatigue since it appears to be an important 
feature of low stress level fatigue and is not a sig- 
nificant aspect of plastic deformation leading to 
failure in static testing. 

Because of the distinctly different modes of de- 
formation, it is not reasonable to expect that a 
quantitative theory of fatigue may be developed 
which will be valid in a logical manner in both 
ranges. If one simple theory appears to do so, it 
must be considered a coincidence. It may be 
observed in passing that even if a crack develops 
in the F range, when it grows, the local stresses in 
the immediate vicinity will be relatively high be- 
cause of stress concentration. The result will be 
that at some stage during the propagation, a crack 


Growth steps in high mens Growth steps in medium loads 


100 microinches 55 microinches 


Growth steps in low loads | | 
17 microinthes 


Fig. 4. Electron micrographs of progressive fracture on microscale. 


which started in the F range will behave more as 
if it were in the H range. Consequently, the region 
of validity of a quantitative theory of high stress 
level fatigue may be expected to extend somewhat 
beyond the strictly H range. It will also be shown 
that by choosing two of the unknown constants 
judiciously, the range of practical validity of a 
theory of fatigue developed for the H range may be 
extended artificially to the F range also. 

There is evidence to show that crack propagation 
in fatigue occurs intermittently.6.7 Christensen, by 
focusing attention on the cross sections of broken 
test specimens in fatigue, was able to show that 
the rings characteristically observed in fatigue were 
in fact an aggregate of a large number of finer rings 
(Figs. 2, 3, and 4). It is suggested that the growth 
time involved by these rings constitutes a ‘“‘dormant 
period’”’ in fatigue. Christensen’s pictures show 
that the dormant period is substantially larger 
than the active period, and that the incremental 
increase in length of the crack is much smaller 
during the former. He finds that, typically, the 
incremental growth of the crack in the dormant 
period, which is dependent upon the stress, is of 
the order of 20 to 100 microns and that correspond- 
ing crack growths in the active period are several 
orders of magnitude larger than these values. The 
number of cycles during which the crack is dormant 
is typically about 50 to 500 cycles and the active 
period approximately 1 cycle. It seems that, 
macroscopically smooth as the crack growth may 
appear to be, it is in fact growing in a stepwise 
manner (Fig. 2). This observation is of funda- 
mental importance in the theory to be presented 
here. 

The discussion that follows is restricted to fatigue 
failure of materials in the H range. Since plastic 
deformation in this region proceeds primarily 
through coarse slip which occurs because of activa- 
tion of the Frank-Read sources, any characteristic 
stresses associated with these sources may be 
expected to play a prominent role in the theory. 
In view of certain assumptions made, the theory 
as developed is not expected to be valid below 50 
cycles on the high stress range and, strictly speak- 
ing, beyond about 5 million cycles on the low stress 
side. Since the range is sufficiently large, it is 
hoped that the theory will be of some use in explain- 
ing the various aspects of fatigue in a systematic 
manner from a unified viewpoint. 


Basic Assumptions 


The material is considered to be a polycrystalline 
aggregate susceptible to strain hardening with a 
known statistical distribution of grain size. The 
mean grain size is L,, and the distribution is assumed 
to be log-normal in accordance with Beck’s observa- 
tions reported by Lunchick.* While there may 
be a number of cracks generated by the fatigue 
stresses, it is assumed that the final failure is due 
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to the growth of a ‘dominant crack.’’ In accord- 
ance with some of the existing findings, this dominant 
crack in the H range is assumed to grow intermit- 
tently and consists of active and dormant periods. 
The theory concerns itself with developing quan- 
titative expressions for the growth of this dominant 
crack in terms of the applied stress, the stress range, 
and the number of cycles of the applied stress. As 
a corollary, this will result in most of the other 
quantitative expressions that we desire in order to 
understand the various aspects of fatigue. 

While these assumptions are based substantially 
on physical observations, we shall make additional 
assumptions based on broad and plausible con- 
jectures. It is assumed that the dormant period 
in the crack growth is associated with the strain 
hardening of a plastic region or enclave that is 
presumed to be formed ahead of the crack. Com- 
pletion of the hardening of the enclave is assumed 
to be followed by a rapid growth of the crack 
through an incremental length represented approxi- 
mately by the width of the hardened enclave. It 
has been found that the hardening process in 
reversed stressing is less efficient than that occurring 
during static stressing.‘ In static stressing the 
process of strain hardening is accomplished by pro- 
ducing an network of dislocations 
whose stability is assured by the applied stress. 
Under repeated loading in fatigue when the stress is 
reversed, presumably some of the dislocations 
generated during the increasing part of the stress 
cycle go back to their sources and get absorbed. 
In the subsequent cycle, more are generated, and 
from them additional ones are prevented from 
going back to their sources. This marginal increase 
in the number of dislocations that interact and 
become fixed is apparently a decreasing function 
of the number of cycles. When the hardening is 
completed, a stable, intersecting array of disloca- 
tions is set up in the sense that their stability is 
not dependent upon the applied stress. It is 
postulated that this motion of the dislocations 
leads to interactions, presumably leading to sub- 
microscopic cracks.’ It is therefore suggested 
that growth of the dominant crack at the end of the 
hardening stage is accomplished through the process 
of these submicroscopic cracks joining together 
with the dominant crack to extend the crack through 
the region of the hardened enclave. This may be 
termed “advanced submicroscopic crack initiation 
ahead of the dominant crack in fatigue.” It ap- 
pears, therefore, that inefficient as it may be in the 
rate of hardening, large amplitude stress reversals 
play a fundamental part in the fatigue mechanism. 
Since the number of submicroscopic cracks pro- 
duced seems to be dependent upon the extent and 
number of dislocation interactions and since this 
in turn is dependent upon the stress range, range of 
stress appears to be an essential parameter in 
fatigue. 


intersecting 


These conclusions aré in a sense an 
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extension of Wood’s hypothesis of fatigue in the 
H range. This conjectured type of crack growth 
in fatigue has in fact recently been observed.’ 

It has already been stated that since the Frank- 
Read sources appear to play a prominent part in 
deformation in the H range, any characteristic 
stresses associated with them may be important 
in the theory. Operation of a Frank-Read source 
causes a pile-up of dislocations in a slip plane, 
and such a pile-up has been shown to result in 
large tensile stresses. It has also been shown 
that when the number of dislocations in the pile-up 
the order of 1,000, the coalescence of the 
leading-edge dislocations may take place leading to 
an incipient crack. This suggestion was made 
originally by Zener." It can be shown that the 
number of dislocations in a pile-up is given approxi- 
mately by the relation® 


is of 


n = [2Ln(t. — (1) 
where » is number of dislocations in a pile-up, 
L» is grain size, t, is resolved shear stress® in slip 
plane in direction of slip, 7; is a mean internal stress’ 
that opposes motion of dislocations, is shear 
modulus, and 6 is a characteristic lattice spacing 
called Burger’s vector.® 

The tensile stress o,, produced at the leading- 
edge dislocation of the pile-up (Fig. 1) may be 
shown to be of the order 
b 

The classical theoretical 
strength!? of a material is of the order E/2a where 
E is Young’s modulus. 


where K, is a constant. 


Assuming that o, and «; 
are the tensile stresses for which 7, and 7; are the 
appropriate shear stresses, we may nondimensionally 
express the tensile stress oy as 


m ~ Ke 


(3) 
with Ky again aconstant. It appears that the value 
of the mean internal stress o; is of vital importance 
in any problem of plastic deformation, and so it is 
in fatigue. If the fatigue failure is due to the 
nucleation of cracks produced by the generation and 
motion of dislocations, since these respond only 
to the stress difference (¢, — o;), we assert that 
the applied stress by itself is of no importance, 
but it is the difference between the applied stress 
and the mean internal stress that is important. 
Also, since it is conceived that the propagation of 
the dominant crack is due to the joining together 
of the cracks generated within the enclave and the 
dominant crack, the response is to (¢, o;)? and 
not simply to (¢, — o;). The exponent 2 is essen- 
tially a factor that gives added weight to higher 
stresses. It is therefore proposed that the right- 


ha 
qu 
Or 
wl 
ot 
th 
th 
ha 
re 
ca 
in 
co 
| cr: 
G 
as 
to 
at 
T 
cr 
in 
ul 
al 
al 
te 
f¢ 
W 
dl 
iI 
— P 
I 
a 
n 
t 
( 
t 
( 
1 


hand side of Eq. (3) must also find a place in a 
quantitative theory of fatigue. 

In order to derive an expression for the rate of 
crack propagation, we need two important relations. 
One is an estimate of the width of the plastic enclave 
which gives the incremental crack growth, and the 
other is an expression for the rate of hardening of 
the enclave. These will enable us to determine 
the rate of crack propagation. If, in addition, we 
have an estimate of the initial crack length and a 
relation for the limiting crack length when the 
catastrophic failure of the specimen occurs, we are 
in a position to define the fatigue failure expression 
completely. The estimate of initial length of the 
crack and its limiting length are provided by the 
Grifith theory of cracks. Specifically, we shall 
assume that unless the crack length is at least equal 
to the equivalent Griffith crack length for failure 
at ultimate tensile strength, it will not propagate. 
This would in essence mean that the dominant 
crack will be that which would have caused failure 
in a simple tension test. This assumption is 
untenable in low stress level fatigue. We shall 
also assume that the limiting length of the crack is 
also governed by the Griffith theory. Thus, if the 
tensile stress amplitude in fatigue is o,, we write 
for failure 


hes = Cult? cr, (4) 


where /,, is critical length of fatigue crack for stress 
o,, and /, is critical Griffith crack length for failure 
at ultimate tensile strength o,,,. This equation is 
an alternate statement of the Griffith expression’? 


= KsVES/ler (5) 


where Kg is a constant and S is the equivalent 
surface energy.'® Regarding rate of strain harden- 
ing, in the absence of more pertinent information 
we shall assume that the hardening obeys a rate 
process and that it can conveniently be measured 
by collapse of the hysteresis loop of the material in 
the enclave. This assumption is not inviolate 
and can be substituted by a more accurate, experi- 
mentally observed expression. 


Quantitative Development of the Theory 


With these assumptions, we are now in a position 
to formulate the theory quantitatively. In a 
recent paper on analysis of stress distribution in 
the vicinity of a crack based upon some calculations 
of Eaton and Stimpson,'* Williams estimates 
the extent of plastic enclave produced ahead of the 
crack in a material that has a well-defined limiting 
octahedral stress, the assumption being that as 
soon as the stress exceeds this value, the distribution 
changes from elastic to plastic. The expression 
that he obtains for the area of the plastic enclave 
is of the form 


A = (6) 


where A is area of plastic enclave, / is crack length, 
og. is applied stress far away from crack, 7,.; is 
limiting octahedral stress, and AK, is an appropriate 
plane stress constant. While shape of the enclave 
is a function of applied stress,‘ it is a satisfactory 
approximation for the present purpose to assume 
that the width of the enclave W, is equal to (A)'”?, 
Hence, for a crack in a sheet specimen, we write 


W, = (7) 


where oc, is appropriate tensile stress which produces 
limiting octahedral stress 7,,;. For three-dimen- 
sional specimens a similar relation may be expected 
to hold good except that the constant of multi- 
plication will be different. 

The disadvantage of Eq. (7) is that it does not 
take into account the strain hardening charac- 
teristics of the material. Since the stress de- 
pendence in Eqs. (7) and (3) is the same and since 
Eq. (3) describes a realistic state of affairs in a 
material capable of strain hardening, it is proposed 
that the width of the plastic enclave for materials 
susceptible to strain hardening be written intuitively 
in the form 


(8) 


where C, is now a constant to be determined experi- 
mentally at this stage of development of theory. 
This intuitive generalization is a fundamental 
assumption of the theory. If the test specimen is 
notched with a theoretical stress concentration 
factor Ky, it will be necessary to modify Eq. (8) 
so that o, will represent actual macroscopic stress 
in the vicinity of the notch. In such a case, we 
shall use the results of Hardrath and Ohman* to 
write o, as 


- Es 
= {1 + (Kr — I) (9) 


where E, is secant modulus of stress-stra’n curve 
at stress o,, E. is ordinary Young’s modulus, Kr 
is theoretical stress concentration factor, 
is stress far away from the notch. 

When we proceed to introduce the range factor 
[(o, — o,')/o;|? into the expression for the width 
of the plastic enclave, we can identify enclave 
width with the incremental crack length A/. We 
may therefore write 


b CO; 


The exponent two and the nondimensional form 
of the range factor are introduced intuitively and 
justified by test results. When the crack extends 
through this incremental length, the new incre- 
mental length when the crack propagates again 
will be given by 


and 
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Fig. 5. Schematic representation of collapse of 
hysteresis loop during strain hardening of plastic enclave 
ahead of the crack. 


Ah = C,(l + Al) a (11) 
b E 


The significance of this expression is clear. It 
states that at this stage, while we do not know 
when the crack suddenly becomes active, we do 
know the approximate distance through which it 
moves when it becomes active. 

Before we can derive the rate of crack propagation, 
we need a relation for the rate of hardening of the 
plastic enclave, for unless the enclave hardens, 
we cannot (strictly speaking) identify it with the 
incremental crack length. At the present stage 
of development of the dislocation theory, it is 
extremely difficult to determine this from purely 
theoretical considerations. Therefore, as indicated 
previously, we shall assume that a rate process 
law governs the hardening. We shall take the rate 
of decrease of the area of the hysteresis loop—that 
is, its collapse—as a measure of this hardening. It 
may be noted that there will always be a residual 
hysteresis loop, for—no matter how small—there 
will always be some free dislocations that will 
execute a to-and-fro motion and thus dissipate 
energy, apart from other dissipating mechanisms 
such as grain boundary relaxation, stress induced 
ordering, etc. We therefore invoke a relation of 
the type 


log. (W/Wy) = KsAN, (12) 


where A, is number of cycles it takes for collapse 
of a loop of area W to one of area Wy (Fig. 5). It 
is evident that 


AN, = Cz log, Ki 
Jo — Go 


(13) 
where Ky, (2, and AK are constants. A, width of 
the collapsed loop, is dependent upon the final 
dislocation structure in the hardened enclave. If 
AN, is number of cycles for hardening of the enclave 
or what amounts to the same, the dormant period, 
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and AN, is number of cycles for the active period, 
we have 


AN, > AN, (14) 


In fact, Christensen finds that they are of the order 
250 to 1. Therefore, the total number of cycles 
to failure V may be written to a good approximation 


N = NpAN, (15) 


where ”p is number of dormant periods (Fig. 6). 
Incremental propagation rate A//AN (AN = AN, + 
AN2) may therefore be approximated to Al/AN,. 
Rate of crack propagation may therefore be written 
from Eqs. (10) and (13) as 


(* 
Al E 


dN AN; (* = 6) 
> 
K 


where C, is again a constant. The dependence of 
rate of crack propagation on crack length is too 
well known for any additional comment. 
Frost and 


Liu,” 
Dugdale,'*"* and several others have 
observed this relationship and have given theoretical 
interpretation of the result. The main difference 
between Eq. (16) and its counterparts suggested 
by others is the proposed additional dependence on 
grain size, stress and its range, and stress o;._ Inte- 
grating Eq. (16), we obtain 


(17) 


where /, is crack length appropriate for N,. We 
shall interpret, as mentioned already, that /, is the 
equivalent Griffith crack length for failure at the 
ultimate static tensile strength. Therefore, JN. 
|/4 and, as a matter of convenience, we shall neglect 


0 


it compared to the large number of cycles to fatigue 
important 


failure. An consequence of this as- 
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Fig. 6. Schematic representation of crack propagation in high 
stress level fatigue. 
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sumption regarding /, is that, until the crack gener- 
ated by fatigue is at least equal to /,, no significant 
general damage has occurred in the material. The 
macroscopic residual static strength will then be 
governed by the pre-existing Griffith cracks or the 
cracks generated during determination of the tensile 
strength, rather than by the fatigue cracks. It is 
entirely possible that the number of cycles it takes 
for the fatigue cracks to grow to the equivalent 
Griffith crack length can be substantially large 
in the low stress level fatigue (F range). Hence, 
it is here that we explicitly invoke the concept of 
high stress level fatigue. In the low stress level 
fatigue, this approximation is not justified since 


the “crack initiation stage’? may be expected to 


be relatively large compared to the total cycles to 
failure. Therefore, we have to take into account 
N, and introduce a more appropriate value of /,. 
Otherwise the expression would predict cycles to 
failure that are less than what the tests would 
indicate. Since the grain size is not frequently 
measured, we may introduce a new constant 


C = CL,,/b (18) 


and write 


Eq. (17) or (19) may be considered the funda- 
mental equation of this theory of fatigue. In 
view of Eq. (4), it is a simple matter to determine 
the cycles to catastrophic failure. The condition 
clearly is 


E 
log, = loge = C= x 


so that 


where N; is number of cycles to failure at stress 
level ¢, and its appropriate range. 


Comments Regarding the Value of o; 


Before determining the constants C and K and 
examining the explanation of the several observed 
features of fatigue, some comments are in order 
regarding o;. In the dislocation theory of metals, 
7; is assumed to be the stress that must be over- 


come before dislocations can move. This stress 
is supposed to arise from various sources, such as 
the existing network of dislocations, solute atours, 
etc. For certain aspects of the dislocation theory, 
the additional assumption is made that the material 
may be looked upon as consisting of potential hills 
and valleys’ and that the maximum slope on all 
these hills is the same and has the value o;. The 
dislocations with their associated stress fields may 
then be looked upon as climbing 1p and down these 
potential hills and valleys. When external stresses 
are applied, they drive the dislocations up and down 
these potential hills and valleys. Since the material 
does not exhibit plastic deformation until the dis- 
locations move and since the external stress has to 
exceed a; before dislocations can move, it is natural 
to interpret o,; as the proportional limit of the 
material. But in fatigue this interpretation causes 
difficulty for two reasons. First, determination of 
the proportional limit for a material is an elusive 
thing, and its value depends upon the accuracy 
and technique of measurement. Second, it is 
well known that if a test specimen is subjected to 
a fatigue stress amplitude somewhat less than the 
proportional limit commonly measured in engineer- 
ing practice, the specimen is liable to fatigue failure. 
We may therefore return to Bairstow’s original 
hypothesis” and assume that the real or natural 
elastic limit of a material is the same as the endur- 
ance limit exhibited by the material under com- 
pletely reversed fatigue stressing. With this inter- 
pretation of the o;, the contradiction of having 
fatigue failures at stresses below the normal pro- 
portional limits is removed. This interpretation 
has the natural advantage that variations of the 
endurance limit under different conditions, such as 
understressing and overstressing, can be discussed 
in terms of the variation of the internal stresses 
due to changes in dislocation structure. A simple 
example of this is the effect of cold working. It is 
fairly easy to show that the effect of cold work is 
to produce an increase in the dislocation density 
resulting in an increase in the value of o;. The 
result of this on Eq. (21) will be to increase the 
value of N as compared to that of an annealed 
material for the same applied stress conditions. 
The constants C and K have to be determined 
from test results at this stage of the theory’s de- 
velopment. Thus, given the test specimens of 
any material, we perform fatigue tests on a statistical 
basis choosing stress levels that extend all the way 
from ultimate tensile strength down to endurance 
limit. Taking two points in the intermediate range 
of stress, we can determine the constants C and K. 
By a judicious choice of the points, the region of 
validity of this equation may be extended well 
beyond the H range. This artificial extension of 
the range may well destroy the physical significance 
of the constants C and K in the sense that, instead 
of K implying strain hardening, it may imply strain 
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softening. But the indications are that the practical 
utility of the expressions is not destroyed by such 
a process. 

Test results of Illg?! for 2024-T3 have been used 
to determine the constants C and K in Eq. (21). 
The chosen test points are indicated by the big 
squares in Fig. 7 with the computed curve super- 
posed on the test points. The same constants 
are used to compute the curves of Eq. (19), the 
results of which are presented in Fig. 8. Ina sense, 
the curves of Fig. 8 may be considered the funda- 
mental curves of the proposed theory of fatigue 
since from them stems the discussion of almost 
all other aspects of fatigue. 

At this stage it is desirable to sum up briefly 
what we have done so far. Fatigue failure of a 
specimen is attributed to the growth of a dominant 
crack. This growth is assumed to occur in jerks 
and is affected by the hardening characteristics of 
the material and the range of applied stress. On 
this basis, an expression for the growth of the 
dominant crack as a function of the number of 
cycles was derived. Failure is assumed to occur 
when the growing crack just becomes critical for 
catastrophic failure at the applied stress. This 
condition determines the number of cycles to failure. 
The fundamental equations of the theory, Eqs. 
(19) and (21), contain essentially four constants of 
which two are the ultimate tensile strength and 
the endurance limit. The other two are the con- 
stants C and K which have to be determined from 
statistically reliable fatigue tests. It will be found 
that, if the statistical mean test points chosen to 
determine the constants are spaced suitably, the 
region of practical validity of the equation is ex- 
tended well beyond the H range. 


Dependence of Fatigue on Different Variables 


Shape of Fatigue Curve 


The expression for the shape of the fatigue curve 
is given by Eq. (21), and Fig. 7 shows a typical 
result where the test points for determining the 
constants have been obtained from Illg.2!_ While 
the theoretical curve appears to fit the test results 
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Fig. 7. Predicted curve plotted against test results. (Material: 
2024-T3) 
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in a satisfactory manner, in the high stress level 
range the predictions will be somewhat spurious. 
The actual slope of the crack propagation curve 
dl/dN has been approximated to d//dN, by virtue 
of the physically justifiable assumption AN. << AN. 
But, in fact, as we approach closer and closer to the 
oun, AN, will become smaller and smaller and 
comparable in magnitude to AN»: Hence, this 
assumption will no longer be justified, and the 
consequence will be that a higher propagation rate 
will have been assumed with the result that the 
predicted life to failure may be less than what may 
be expected. For this reason the upper limit of 
stress validity of the theory has been restricted 
arbitrarily to about 50 cycles. Therefore, the pre- 
dicted results in this region must be used with some 
caution. The low stress side of the curve and the 
artificial extension of the region of practical validity 
have already been discussed. 


Existence or Nonexistence of Endurance Limit 


This property is closely related in principle to 
the value of the mean internal stress o;. There 
are two alternate possibilities for the existence of 
a well-defined endurance limit. Since it is generally 
accepted that the fatigue process generates its own 
cracks, a well-defined endurance limit would imply 
either that the cracks’ formation is prevented or 
that, if they form, their propagation is inhibited. 
Since formation of the cracks is due to dislocation 
movements and interactions, a natural way of 
obtaining a well-defined endurance limit is to prevent 
the motion of dislocations. This requires the 
pinning of existing dislocations to prevent them 
from moving and also the pinning of any new dis- 
locations generated in isolated places as the result of 
local stress concentrations. This anchoring can 
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be effected fairly easily in interstitial solid solutions 
(mild steel, for example) through the formation of 
Cottrell consequently, one 
should expect a well-defined endurance limit in 
such materials. This fact is well known. It also 
follows that, since the formation of these atmos- 
pheres is a diffusion controlled process and hence 
time dependent, the endurance limit of such ma- 
terials should be frequency dependent. It also 
follows that, if the interstitial solute atoms of such 
materials are removed by a suitable treatment, the 
endurance limit must disappear. 
fact.**> The second alternative—preventing propa- 
gation of existing cracks—is somewhat more difficult. 
One possibility for this is an indefinite capacity for 
plastic deformation—that is, an ideally plastic 
material. The enclave in such a material will not 
strain harden, and hence the crack cannot propa- 
gate. This alternative appears to be one of theo- 
retical interest only. What this amounts to is 
that, apart from interstitial solid solutions, one 
may not in principle expect a well-defined endurance 
limit. An approximate endurance limit may be 
defined on the basis of nonfailure at a large number 
of cycles, say of the order of 500 million. 


atmospheres,** and, 


This is a known 


Residual Static Strength 


The determination of the residual static strength 
after an arbitrary amount of fatigue is straight- 
forward when the curves of Fig. 8 are computed. 
It is clear from Eq. (4) that the scale of crack 
length can be converted to a corresponding scale 
of the critical Griffith stress. Thus, we may convert 
the scale of log ///, to that of the residual static 
strength by the relation 


Ors = (22) 
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where ¢;,; is residual static stress at which specimen 
would fail. Residual static strength in pounds 
will be given by 


= (23) 


where A, is area of cross section of test specimen. 
Graphical construction for the residual static 
strength appropriate for o vs. N curve of Fig. 7 is 
indicated in Fig. 9. When size of the crack is 
significant compared to initial area of cross section 
of the test specimen, a correction factor will be 
necessary in these expressions. Analytical repre- 
sentation for the residual static strength aiter 
application of m cycles of stress at stress level ¢ 
and range (0, — o,') is given by 


1/3 
0d = Cult ( ] ) = Cult exp x 
1 
Og. 


A parametric representation of Eq. (24) for 
different values of o, is given in Fig. 10. Some 
inherent limitations of this method of representation 
of the residual static strength are pointed out in 
the original report** on which this paper is based. 
The number 1/4 is reintroduced here in order to 
satisfy the condition that when o, = o,, and 


n 1/4, Ors = Curt 


Effect of Frequency 

The effect of frequency is a_ straightforward 
consequence of the theory and has already been 
treated elsewhere.” It is not yet possible to treat 
it in a quantitative manner, but qualitatively it 
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can be looked upon as due to the possibility of 
thermal activation of some partially activated 
Frank-Read sources in the enclave. Briefly, the 
principle of approach may be outlined as follows. 
The stress concentration at the leading-edge dis- 
location of a pile-up of dislocations is dependent 
upon the number of dislocations piled up in the slip 
band. When the external stress o, is applied, of 
the Frank-Read sources with a spectral distribution 
of lengths n(/,)d/,, those with lengths defined by 
the condition 


1. 2 pb/(o, — a3) (25) 


where /, is source length, become immediately 
activated. Sources with lengths smaller than this 
value are only partially activated in the sense that 
they assume shapes less than the critical semi- 
circular shape. Under suitable thermal fluctua- 
tions, these partially activated sources can become 
momentarily active, send out an additional dis- 
location loop into the slip band, and hence cause 
more local stress concentration. It is presumed 
that this increase will result in a decrease in fatigue 
life. The process is similar to the exhaustion 
theories of creep, and the measured activation 
energies (34,000 cal/gm mole for aluminum) seem 
to support a creep-type mechanism since they are 
of the same order as volume diffusion energies. 
The frequency effect comes into the picture by 
virtue of the fact that if the frequency is high, 
not enough time is available for successful thermal 
fluctuations, with the result that the number of 
marginally active sources becomes less and produces 
fewer sources of stress concentration and pre- 
sumably, therefore, a higher fatigue life. Some 
conclusions arrived at on the basis of this approach” 
are briefly summarized : 

(1) The effect appears to be marginal. But 
there are experimental observations which show an 
increase in life of about 800 percent for high-purity 
aluminum” for a three-decade increase in frequency. 
This would apparently mean that other, more 


simply activated mechanisms may also be operative 

(2) The effect is more pronounced at higher stress 
levels than at lower stress levels—that is, what 
one may call the endurance limit is not sensibly 
effected by the frequency. This would imply 
that the possible frequency dependence of the 
endurance limit of the interstitial solid solutions 
(mild steel, etc.) may be insignificant. 

(3) Cold worked materials show more pronounced 
frequency effect than annealed materials. 

(4) The effect of frequency for a given temperature 
is significant only within a certain frequency range. 
Beyond this range, changes in frequency do not 
result in significant additional changes in fatigue 
life.” 


Effect of Grain Size 


This effect appears through variations in Lp. 
Larger grain size implies larger crack propagation 
rates, Eq. (16), and hence shorter lives. Coarse 
grained materials may therefore be expected to 
have shorter lives. This fact is well known.” 


Effect of Cold Work 


To predict quantitatively the effect of cold work 
from purely theoretical considerations is somewhat 
difficult since it is sensitively dependent upon the 
dislocation structure. Approximate calculations** 
based upon dislocation theory predict more signifi- 
cant effect of cold work than is normally observed. 
The effect occurs for two possible reasons. One 
is the increase in the density of dislocations’ due 
to the cold work, and this in turn may be shown? 
to increase the mean internal stress o;. It is 
evident that the result of this will be to reduce the 
net effective stress (¢, — o;) responsible for fatigue 
damage, and hence the fatigue life would be higher. 
The second effect is due to the introduction of 
compressive stress fields on the surface where cracks 
generally occur, thus inhibiting their growth. 


Effect of Temperature 


The dependence on temperature does not appear 
explicitly in the theory, but it is fairly easy to see 
how it occurs. In the first place, an increase in 
temperature means a decrease in the modulus 
FE and hence a higher rate of crack propagation. 
However, the second and probably more important 
effect of temperature is in its capacity to reduce the 
mean internal stress level. The substantially stable 
array of dislocations that exists at the room tempera- 
ture will be upset at the elevated temperatures 
where creep-type mechanisms can operate. Ther- 
mally activated dislocation climb’ and _ possible 
mutual annihilation of these activated dislocations 
can reduce the value of o; with the result that the 
crack propagation rates will again become larger. 
If the temperatures are higher than the normal 
recrystallization temperatures of the material, the 
temperature effect occurs through the grain growth 
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during testing leading to larger values of the grain 
and, hence, again a decrease in the fatigue life. 
While qualitatively the dependence of fatigue 
through the variation of E, o;, and L,, is fairly easy 
to understand, we are not yet in a position to take 
it into account quantitatively in a satisfactory 
manner, since the dislocation theories pertaining to 
these aspects are not too well developed. 


Effect of Changing Test Stress Ratio 


The effect of changing the test ratio may be 
investigated in a straightforward manner. Intro- 
ducing the test stress ratio R = o,'/o,, we can re- 
write the expression for cycles to failure in the 
alternate form 


N 2 log. (oy Fo) loge [(o — Cin) /K| 


Evin 


(26) 


where oj, is “endurance limit’ appropriate for 
given test stress ratio. Thus, for example, when 
R = —1, the minimum stress is in compression and 
is exactly equal in magnitude to the maximum 
tensile stress. Test ratio RK = 1 implies no fatigue 
stress. If the form of the equation has general 
validity, it should be possible to compute the 
constants C and K for one value of R (say, R = —1), 
and the same constants must be valid for different 
values of R for the same material, provided appro- 
priate values of the endurance limit corresponding 
to the different test stress ratios are used. No 
exhaustive or statistically reliable test results are 
available to explore this idea thoroughly, but 
Grover ef al.** present test results in which the 
range was varied from KR = —1 to R = +0.80. 
Figs. 11 and 12 show the comparison between the 
theory and their test results which may be considered 
fair. More exhaustive tests are necessary to con- 
firm these different aspects. 


80 Aerospace Engineering +* October 1961 


Cumulative Damage 


The generally observed results of cumulative 
damage are among perhaps the most direct conse- 
quences of the theory. The Miner theory of 
cumulative damage suggests that, when different 
stress levels are applied for a different number of 
cycles, total life to failure is governed by the condi- 
tion that 


q 


t=1 (27) 


q=1 Ng 


where ”, is number of cycles applied at stress level 
g,, and N, is total number of cycles to failure at 
this constant stress level. It has been evident 
from the beginning that this theory is inadequate 
in the sense that the experimentally observed sum 
is rarely equal to 1 as the theory predicts. Further- 
more, it is well known that if two stress levels o, 
and o(o; < o2) are used in a two-step cumulative 
damage test, the sum of the partial life fractions 
is very much dependent upon the order in which 
the stresses are applied. If the low stress is applied 
first, the sum is generally found to be less than one, 
and if the high stress is applied first, the sum, is 
found to be greater than one. This sensitivity to 
the order of stressing is not explainable on the basis 
of Miner’s theory. More recently, several theories 
of an essentially empirical nature have been pro- 
posed in order to explain this phenomenon. Among 
them are those of Freudenthal and Heller,*! Corten 
and Dolan,*? Liu and Corten,** and Newmark.* 
These theories suggest, in general, that the damage 
is not simply proportional to u,/N, but to some 
power of this fraction. All these theories have 
been able to explain some of the different aspects 
of cumulative damage with varying degrees of 
success. 

It appears, however, that a cumulative principle 
sensitive to the order of stressing is a natural 
consequence of the present theory. The theoretical 
model has in it the natural means of extension of 
the theory to program-type loading, random loading, 
and—as a formal extension—the theory of acoustic 
fatigue. In view of the importance of the cumula- 
tive damage aspects of fatigue to the aircraft 
industry, the theory will be developed here in some 
detail. 


Two-Step Cumulative Damage 


Consider the application of two stresses o; and 
o2 with the corresponding stress ranges — ) 
and (o2 — o»’). Assume o; < oe. With this 
information and the statistically reliable results 
of a basic curve of stress vs. number of cycles, we 
can construct curves of Eq. (19) for these two 
stresses and the appropriate ranges (Fig. 13). Let 
us now consider the application of m cycles at stress 
level o; and its appropriate range. As a conse- 
quence, the initial crack of length / will grow to 
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length /; as indicated by the vertical line with the 
open arrow. The crack of the same length can be 
grown by mz cycles at stress level o2 and its appro- 
It is easy to see from Eq. (19) that 
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priate range. 


where the notations o;; and oo; indicate that when 

o1| | —oy'| and | o2| | then ¥ 9; 
—that is, if there exists a mean stress, the endurance 
limit is not necessarily the same, as indicated in 
the discussion of the effect of test stress ratio. We 
may, therefore, write 


Insofar as the test specimen is concerned, it is of 
little consequence how the damage is caused. An 
irreversible damage exists in the specimen, and the 
relationship between the number of cycles of stress 
at different stress levels and their appropriate 
ranges for the same amount of damage is given 
by Eq. (29). In a sense, this is the fundamental 
equation for the discussion of different aspects of 
cumulative damage. If now we change over to 
stress o2 and its appropriate range and subject the 
test specimen to fatigue failure at this stress level, 
the specimen will fail when the crack reaches the 
appropriate critical length for failure at this level. 
The progress of damage therefore follows the line 
associated with stress 2 (indicated by open arrows, 
Fig. 13) and terminates at point No. Therefore, 
the additional number of cycles to failure sustained 
by the specimen is (N2 — me), and the total number 
of cycles to failure is 


Nw = n + (No — m) (30) 


The sum of the familiar damage fractions is 
given by 


It may be noted that in the normal practice 
n/N, is plotted against m2/ of the present notation. 
In the normal notation this is the same as plotting 
m’/N, against 1 — (m’/Ne) where m’ and nm,’ are 
the partial lives at the two stress levels. Sub- 


stituting appropriate values of N, and N2, it may 
be seen that 


N, 


~ Stress 9, 


——/  Range(o,-9, ) 
| 
| 
4 
2 Stress 
Fig. 13. Schematic $ 1, } { A 
. / 
representation for dis- = | 
cussion of two-step 
cumulative damage. a | / 
| / 
A 
| q 
| 
| // 
Me Log % 


Crock length 


N, N, 2 1 1 (2 
Og, 
g g K 
(32 
E i 
= am (33) 


lo (=) lo (2 
Se Se K 


It is clear from the structure of the two equations 
that a; and a: are predetermined constants once 
the stresses (and their ranges) are given. In similar 
manner, we can see that if the order of application 
of stresses is reversed and o2 is applied first for m 
cycles, we obtain 


Ny = n+ (N — m) (34) 


where 7”; is equivalent number of cycles to failure 
at stress level o; to cause the same damage as that 
due to 7 cycles at stress level oo. 
partial life fractions are 


The corresponding 
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It is immediately seen that the results of this 
theory are different from those of Miner’s theory. 
In Miner’s theory, a; = 6; and a2 = 62. The sum 
of the partial life fractions when stress a; is applied 


first and then oo is Sj = ayn + (1 — aon). There- 
fore, 
nC E 
So = x 
y = 
Og, 
K 
(log, — log, 
(37) 


uit 
log, ( ) log, ( ) 
02 


Similarly, = 1 + — 62). Therefore, 


(log, — log, o;) 


Ou u 
log, ( log, 
02 


It is easy to see that since 0; < o2, Spy < 1 and 
So, > 1, a result that is consistent with the generally 
observed values. 


(38) 


As already mentioned, a” is 
plotted against and 6.” against to describe 
cumulative damage. The slopes of the two lines 
thus obtained are given by 


ay (log, — log, o1) 
= = (39) 


dx (log, our — log, 


_ Bx (loge eure — loge 
21 


(40) 
dx Bi (log. — log, a1) 


where the notation (dy/dx)y implies the slope of 
the damage line when stress o; is applied first and 
then oe, and vice versa for (dy/dx)x. It is seen 
from these relations that when the low stress is 
applied first and then the high stress, the damage 
line falls above the Miner representation (45° 
line), and the opposite occurs when the order of 
application of stress is reversed. We may also note 
a symmetry relation between the two lines repre- 
sented by the two Eqs. (39) and (40), by virtue 
of the fact that 


(dy/dx)2(dy/dx) = 1 (41) 


This symmetry of the damage curves has been 
observed.* 

An extremely important aspect of cumulative 
damage may be pointed out at this stage. Suppose 
that low stress o; is applied with the expectation 
that at some stage the stress will be changed to 
higher stress o2 (Fig. 13). When the number ot 
cycles at stress level o; reaches Nj’ as indicated in 
the figure, the length of thé crack will have become 
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just critical for stress level 2. Therefore, while in 
single stress fatigue it is meaningful to talk in terms 
of fatigue cycles beyond N,’, in two-step application 
it is not meaningful to do so. The value of N,’ is 
given by the condition 


log log, = 


so that 
Og, og, 
Ni = (43) 
E Oli 
[he maximum value beyond which u has no mean- 


ing in two-step cumulative damage is given by the 
condition 


T? 


n= (44) 


Substituting this value, the limiting points of the 
damage lines in Fig. 14 are found to be 


(log, gui — log, o2) 
(log. — log, o1) 
Qo) l (45) 
5 
= 1 = Xo 
(log, Gur — loge a2) 
DoMmar = = 
(loge — loge o1) 


It is clear that, given the two stresses o, and o 
and the ultimate tensile strength, we can draw the 
damage lines from the information contained in 
Eqs. (39), (40), and (45). Thus, for every pair of 
stresses, we obtain a pair of damage lines. As an 
example, such a damage line fan is drawn in Fig. 14 
for the same numerical data for which curves in 
In the scheme of 
things presented here, it is not meaningful to join 


Figs. 7 and 8 are appropriate. 


these pairs of lines to the point (1, 1) as is commonly 
done In fact, it is difficult to see why these curves 


Fig. 14. Cumv- 
lative damage 
curves for different 
pairs of stresses. 
(Curves computed 
for the same speci- 
mens as Fig. 7.) 
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should converge on this point. As an example, 
test results from Grover et al.” and the predicted 
damage lines for a pair of stresses are given in Fig. 15. 

Generalization to cumulative damage involving 
three or more stress levels is straightforward. It 
is seen that if the stresses are given by a1, a2, 93, 

Am, Anz, Ang are the number of cycles 
of stress applied at each one of these stress levels, 
and if o, is the last stress to be applied (denoted 
as a reference stress), we have the number of cycles 
sustained until failure occurs at the last stress given 
by 

qd 


An, = N,— > (An)q (46) 


q=1 


where 


9 ‘\2 2 
Or — On; Cy Ge Ogi 


and An, is number of cycles sustained at stress 
level o, applied last after applying Am, Anz, An; ... 
cycles at stress levels o1, 02, 03, etc., and N, is cycles 
to failure at stress level o,. 


Program-Type Loading 


In program-type loading, the common practice is 
to apply a sequence of cycles of stress at various 
stress levels and repeat the sequence until failure 
occurs. The problem then is one of determining 
the number of times the sequence can be applied 
before failure occurs. 

Let us, for example, assume that stresses oj, 9», 
03, ... O are applied for Am, Amz, Anz, ... Ang 
cycles. An, is arbitrary; the only limitation is 
that An, < N,, where N, is number of cycles to 
failure at stress level o,. As before, let o, be the 
reference stress. It can be any one of stresses oi, 
2, 03, the only limitation is that o, > the 
endurance limit. During the application of the 
sequence once, the extent of equivalent life con- 
sumed at reference stress o, is given by 


(An), = (48) 
where (An,), is given by Eq. (47). Failure would 
then occur when the crack in the specimen reaches 
a length just critical for the largest stress of the 
sequence. If o, is the largest stress of the sequence, 
the number of cycles (N,), needed at the reference 
stress to reach the critical length is given by 


Cult Or — On 
on, ( K ) 
(N,), = = 


E Cri 


Therefore, the number of times the sequence can 
be applied before failure occurs is 


ns = (N,),/(An.); (50) 


It will be noted that (Am,), is essentially the 
weighted sum of the stress pulses in the sequence. 
Consequently, the result ms is sensitive to the 
number of cycles at each stress level. If the size 
of the block of the program—that is, Am, + An: + 
An; ... An, = Any,—the total number of actual 
cycles to failure is 


Nrp = nsAnrp (51) 


For example, let us consider the tests performed 
by Liu and Corten*® in which o, = —¢,’, so that 
og = oO, = o; The program-type loading used 
by Liu and Corten consists of applying two stress 
levels in blocks of 10,000 cycles in reversed bending, 
varying the ratio of cycles corresponding to the two 
stress levels. Eq. (47) therefore reduces to 


(An,), = An, (f= x 


The value obtained for o; on a statistical basis 
was 20,000 psi, and from among the different test 
levels used by them we shall choose the high stress 
as 50,000 psi, and the low stress 45,000, 35,000, and 
25,000 psi. The percentage of life at high stress 
was taken as 9.95, 3.97, and 0.95 percent. Let us, 
for example, consider that o, = 50,000 psi and 
a2 = 45,000 psi, and let us choose o; as the reference 
stress o,. Therefore, 


50,000-20,000/7 \ 5 


50,000—20,000 
log, 


45 ,000—20,000 
log. 


x 
OL. 
K 


x 
K 


Since the logarithmic terms are so close to each 
other, considerable simplification results without 
serious loss of accuracy if we neglect them. Thus, 


= Am + An X 0.563 


An,, = An, + Ane X 0.563 


The life to failure at the reference stress level of 
50,000 psi on a statistical basis is 18,770 cycles 
based on 72 tests. Therefore, the number of blocks 
to failure is 
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Fig. 15. Response 
to order of appli- 
cation of stress in 
cumulative dam- 


age. 
Predicted curves 
T and test results 
| from NACA TN 2324, 
| poge 32 
| Mati. SAE 4130 
02 0.4 o6 10 
n, 
18,770 18,770 
ns = — = 


Aner An, + Ane X 0.563 


Since the size of each block is 10,000 cycles, the 
actual life to failure is 
18,770 


An, + Any X 0.563 
In one test An; = 995 and An, = 9,005. Sub- 
stituting, we obtain the life to failure in this program 

of loading as 
18,770 
Nrp = 10° X ———————_—— =. 30,800 cycles 
995 + 9,005 & 0.563 

The actually observed mean statistical life for 
tests was 30,550 cycles, thereby giving a percentage 
error of 0.82 percent. Comparing the logarithms 
of life, the error is approximately 0.1 percent. 
Some other computed results are presented in Table 
1 and Fig. 16. It will be noted, therefore, that 
where statistically reliable test results are available, 
comparison between the theory and experiment 
appears to be good. It is desirable to obtain exten- 
sive additional experimental confirmation of this 
and the other aspects of the present theory of fatigue. 


Fatigue Damage in Random Loading 


The case of random loading is a simple extension 
of the preceding equations. We shall assume that 
a frequency distribution function of An, vs. a, is 
given. The problem then is one of determining 
the number of cycles to failure for such a loading. 
The extent of randomness may be one of two types. 
In one case, while the loading is random, there is 
regularity in the loading in the sense that after 
the application of the number of cycles represented 
by the area under the distribution curve, the load 
pulses repeat themselves again in a random manner. 
The second case will be the problem of predicting 
the probability of failure during the anticipated 
design life of a specimen. We shall, for example, 
imagine that the specimen will be subjected to 
Nr stress pulses in a random manner during its 
anticipated life, and thatthe distribution function 
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is known in terms of the maximum, mean, and 
minimum stresses and the standard deviation, 
The problem then is one of ascertaining the proba- 
bility of failure at any stage. 

Considering the first case, let us assume, for 
example, that the stress pulses obey a normal 
distribution function, 


1 Anr -1/,(% 


where a,, is mean stress of the sequence, s is standard 
deviation, Anr is total number of cycles in this 
random sequential loading, and p; is a normalization 
factor introduced by virtue of the fact that the 
maximum and minimum stresses do not go to 
infinity. The equivalent number of cycles corre- 
sponding to An, at reference a, is 


An, = (53) 


An; 


An; — x 


/ 9 
| : do, (54) 


Or Gri Op Ogi 


The logarithmic terms are neglected in view of 
the relatively small change that they cause in the 
final results and the considerable simplification 
Therefore, 
the number of blocks of random sequential pulses 
Anr which may be applied before failure occurs is 


that results from this approximation. 


= (N,)1/An., (55) 


where (N,)z is given by Eq. (49), with the assump- 
tion o, = Omar. Because of the randomness of the 
application of stress pulses, the failure in the ns 
block has a probability of failure whose cumulative 
probability increases from (0 at the beginning of the 
block to 1 at the number of cycles corresponding to 
failure at the minimum stress of the block. It 
must be noted that this probability of failure is 
quite different from the inherently statistical nature 
of the fatigue problem and arises by virtue of the 
randomness of the load application. As before, the 
actual number of random stress pulses which may 
be applied before any failure occurs is 


(N,) 


= Anr = Anr X ns_1 (56) 


ler 

During application of the next Amz cycles, failure 
would occur as indicated in Eq. (56) and the cumu- 
lative probability of failure will apparently be a 
function of the cumulative distribution function of 
the applied stress pulses and the instantaneous 
crack length. This will be discussed in detail 
elsewhere. 

The second problem may be treated as a formal 
extension of this case. We shall now interpret 
An, not as the number of pulses applied in a random 
manner in a block, but as the total number of 
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cycles Vr that are presumed to be applied during 
the anticipated design life of the specimen. The 
minimum number of cycles at which failure may be 
expected will correspond to the critical crack length 
for maximum anticipated stress ¢, = Omar. There- 
fore, the formal number of blocks after which failure 
may be expected is 


NS, min = (N,) L, (57) 


where (N,), is given by Eq. (54) with Nr sub- 
stituted for Any and the appropriate distribution 
function. The actual number of stress pulses after 
which failure may be expected is therefore 


Na = N T NS; min (58) 


It is clear that if ms, min < 1, the failure would 
occur during the anticipated life of the specimen. 
Correspondingly, if s, min > 1, no failure may be 
expected. The mean anticipated life, following the 
same reasoning, is 


N = N ( N,) mean 
Nim 
(Ne), 


Tult Or — Ori 
log, 
= Omean K 


(59) 


where 


The maximum life that may be expected before 
failure occurs is 


Nr [(N maz/(Ne)r] (61) 


where it will be noted that 


N = 
4VYa, mar 
Or — On; 
2 
— o,'\2 
Or Ori Or Or 
E Ori 


Any theory of acoustic fatigue has three essential 
aspects. The first is an adequate quantitative 
representation of the acoustic field. The second is 


( N,) mar 


Acoustic Fatigue 


the determination of the stress response of a panel 
to this given acoustic input. The third is the esti- 
mation of the time to failure of the panel for the 
stress response so determined. The emphasis in 
this section is on the third aspect; it will be treated 
as a formal extension of the present theory of fatigue. 
This is possible since it appears that acoustic fatigue 
is essentially the response of the structure to random 
excitation. 

The stress response of the panel to an acoustic 
input is given generally in either of the equivalent 
alternate forms, its power spectrum or an rms stress 
as a function of the frequency. An additional 
assumption is frequently made that the amplitude 
peaks associated with an rms stress peak obey 
the Rayleigh distribution function given by 


P = (63) 


where &, is rms stress peak occurring at frequency 
w,, and P is probability of occurrence of amplitude 
peak o,”. It is apparent that there is an additional 
complexity in this problem because of the different 
frequencies at which fatigue stress cycles are simul- 
taneously imposed. Because of this, it will be found 
necessary to estimate equivalent times for the 
damage at a reference frequency and stress corre- 
sponding to different times for the same damage 
at the different frequencies and stresses. 

Let us consider a frequency w,, and let the rms 
value of the stress at this frequency be 6,. In 
unit time, the probable number of stress pulses 
An, related with each stress amplitude ¢,” is given 
by 


An, = — exp} (64) 
p2 Gq 2 
where po is a normalization constant that takes 
into account the fact that the maximum and mini- 
mum anticipated stresses are finite. There are 
two additional aspects worth noting in acoustic 
fatigue. The first is that there is generally a mean 
stress over which is superposed the vibratory stresses. 
The second is that the vibrations are taking place 
over a range of frequencies (Fig. 17), and so pre- 


Table 1. Comparison Between Theory and Test Results for Program-Type Loading 


Percent of Predicted Actual Log of Log of 
High Low Life at High Life to Life to Percent Predicted Actual Percent 
Stress Stress Stress Failure Failure Variation Life Life Variation 
50,000 45,000 9.95 30,800 30, 550 0.82 4.4886 4.4851 0 
50,000 45,000 3.97 32,000 33,600 4.72 4.5051 4.5262 0.246 
50,000 45,000 0.95 32,600 31,220 4.42 4.4997 4.4945 0 
50,000 35,000 9.95 85, 100 7,600 2.86 4.9299 4.9424 0.25 
50,000 35,000 3.97 118,500 118,100 0.34 5.0737 5.0722 0 
50,000 35,000 0.95 143 , 500 132,500 8.3 5.1568 5.1221 0.68 
50,000 25,000 9.95 162,500 163 , 700 0.75 5.2108 5.214 
50,000 25,000 3.97 404,000 362,000 1.7 5.6064 5.5583 0.865 
50,000 25,000 0.97 1,135,000 1,045,000 8.6 6.0550 6.0191 0.6 
Mean variation 4.7 Mean variation 0.29 


Material 7075-T6 aluminum alloy. 
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Fig. 16. Comparison between statistical mean values 
in program-type loading. Test results from NASA D-647. 


sumably one has to take into account the effect of 
frequency on fatigue. It has already been men- 
tioned that in high-purity aluminum one may 
observe about a decade change in fatigue life for 
three-decade change in frequency. Apparently, 
for many alloys, the change is much less. Typically, 
in the case of acoustic fatigue, the frequency range 
over which significant rms stress peaks occur is 
rarely more than one decade. This would probably 
correspond to a maximum of 200-300 percent 
change in life, well within the normal scatter band 
of fatigue. So to obtain a reasonable solution to 
the problem of acoustic fatigue, the assumption 
may be made that the effect of frequency as such 
may be neglected. We can now write the governing 
equation for acoustic fatigue (neglecting the log- 
arithmic terms) 


q mqt q mq mri 
(An,), = An, 


(65) 
” ” 


Omr Omar 

where ¢,,; iS endurance limit at frequency of 
vibration w, corresponding to mean stress Om, 
Eq. (65) is exactly the same as Eq. (47) except that 
it is rewritten in terms of mean stress instead of 
stress range. The total equivalent number of 
cycles per unit time at the reference stress corre- 
sponding to stresses ¢, may therefore be written as 


q 


— — Omr) mai: sil 
The limit ¢,,,; is chosen on the basis of the as- 
sumption that any stresses below the endurance 
limit have no effect on the fatigue life. This as- 
sumption does not appear to be rigorously true. 
The stress o,” is an arbitrary upper limit of the 
Rayleigh distribution function. If the reference 
frequency is w,, the time it takes to impose An,,, 
cycles is clearly given by 
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Terqg = Atterg/ (67) 


It must be noted that there are rms stress ampli- 
tude peaks at different frequencies w, and these 
have to be added. The total equivalent time at 
the reference frequency and the reference stress 
level is 


= berg (OS) 


In order to obtain actual time to failure, we have 
to divide time to failure ¢, at the reference stress 
level and the reference frequency by ¢,,. Thus, 
actual time to failure is given by 


(69) 


The time ¢, may be seen to be equal to 


og, , 10g. 
(N,) L 1 OL K 


” 


\ Eomr 


| 


(70) 


Physically, ¢,, is the sum of the equivalent times 
at the reference frequency for the damage to which 
the material is being subjected simultaneously at 
different frequencies and stress levels in unit standard 
time. Essentially ¢,, may be looked upon as a 
coordinate transformation for time. The trans- 
formation is from the frequencies w, and stress 
go,” to the reference frequency w, and reference 
Eqs. (66-68) may be looked upon as 
appropriate coordinate transformations. 


” 
. 


Simplifying Assumptions and Resulting Equations 


Solution to the problem of acoustic fatigue is 
essentially given by Eq. (69). A comment is in 
order regarding the choice of the reference frequency 
and reference stress level. Since the basic stress 
vs. number of cycles curve of fatigue is to be de- 
termined in any case for the panel in order to de- 
termine the constants C and K, and since this is 
most easily done by exciting the panel in its funda- 
mental mode, it simplifies the analysis considerably 
to choose the frequency of the fundamental mode 
as the reference frequency. The reference stress 
may correspondingly be chosen as the largest stress 
to which the panel may be subjected. The mean 
stress of vibration which is generally small may 
be assumed to be approximately the same at the 
various frequencies. This would also imply the 
endurance limits to be the same at the different 
frequencies—that is, 


Cn Cop Car = = Om 


In addition, it will be noted that ¢, is the 
time to failure at the reference frequency and 
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the largest stress o,” which may be expected 
on the test panel. This is a term that can be 
determined fairly accurately from tests on a statis- 
tical basis. In such a case, only the denominator 
of Eq. (69) needs to be computed. We may there- 
fore write Eqs. (66-69) in the simplified form 


q ” WV 

> 1 1 /o,"\*\ 

jo "— ¢ a,’ — oa 5 
qd mi q ” 79) 
” ” \ do, (72) 


One implicit assumption made in the formal 
extension of the present theory of fatigue to the 
problem of acoustic fatigue must be mentioned. 
The basic fatigue theory as developed is strictly 
valid only for high stress levels. In acoustic 
fatigue, the damage is due more to a large number 
of cycles of low stress in short time than to a small 
number of cycles of large stress. Hence, some 
caution is needed in the extension of the present 
theory to acoustic fatigue. It will be noted that 
this item of caution is introduced by using the 
experimentally determined value of ¢, and using 
the general form of the equations to predict times 
to failure in acoustic fatigue. Extensive experi- 
mental verification is needed to substantiate this 
approach. 


Statistical Aspects 

The preceding discussion of the several aspects of 
fatigue was closely related to the value of o;, the 
internal stress that is presumed to oppose the 
motion of dislocations. In this theory the value of 
this internal stress was identified with the endurance 
limit in fatigue by invoking Bairstow’s hypothesis 
which suggests that the endurance limit of a material 
is the same as its natural elastic limit. Apart 
from determining the endurance limit from extensive 
tests, the determination of the actual value of a, is 
an extremely difficult task at this stage. This is 
the main disadvantage of the present theory and, 
in a sense, reflects the bridge between the micro- 
scopic and macroscopic approaches. 

The difficulty of accurate determination of o; 
brings us directly into the statistical aspects of 
fatigue, of which there are three items. First, we 
have assumed that the value o; (maximum slope of 
the potential hills in the material at the level of the 
lattice structure) is a constant. Actually, it is 
clear that this is not true. Therefore, in certain 
regions the plastic enclave will be larger and pre- 
sumably the deformation more severe than in other 
regions. In other words, the width of the enclave, 
and hence the incremental crack length, the rate 
of strain hardening, and the rate of crack propaga- 
tion are indeed statistical variables and not the 
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Fig. 17a (left). Schematic representation of typical rms stress 
response in acoustic fatigue. Fig. 17b (right). Idealized stress 
response curve used in discussion of acoustic fatigue. 


nonstatistical quantities we have assumed in 
developing the theory. 

The second aspect is concerned with grain size. 
While in the beginning we assumed that the grain 
size had a statistical distribution with a mean size 
L,» and a logarithmic normal distribution, nowhere 
did we take this specific property into account, 
although it is not difficult to do so in principle. If 
grain size distribution is n(L)dL, then one can 
define a mean rate of crack propagation 


dln L E 

iN n(L) dL (73) 


Og. K 

and use this value instead of Eq. (16) in subsequent 
developments of the theory. In principle, we should 
do the same thing with o, also. But at the present 
stage of the dislocation theory, we know nothing 
about its distribution, with the result that there 
will always be an element of statistical uncertainty 
in the fatigue results. This will, in principle, be 
reflected in the theory by a spread of the value 
a; This necessarily implies that the familiarly 
defined endurance limit for metals is also a statistical 
quantity, a fact that has already been shown experi- 
mentally by Epremian and Mehl.” 

The third aspect is concerned with the cracks 
themselves. The analysis presented here concerned 
itself with the progress of the “dominant crack’”’ 
only, and presumably the propagation of this crack 
proceeds from the opening up of other smaller cracks 
in its path. In reality, one should expect the 
possibility of having more than one major crack. 
In sheet metal structures, this may not be so serious 
as in other types of three-dimensional structure 
where it may play a prominent role. 

These three aspects therefore seem to indicate 
that for a considerable time to come, the problem 
of fatigue will remain essentially statistical in nature. 
One has to take recourse to extensive tests to obtain 
statistically reliable data. It appears that the main 
contribution of the theory is that, when such 
statistically reliable test data concerning the basic 
stress vs. number of cycles curve in fatigue are 
given, most of the other aspects of fatigue may be 
predicted in a systematic manner without any 
arbitrary assumptions other than those made or 
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implied in development of the basic theory. The 
systematic approach in this method has become 
possible by virtue of the equivalence relations 
obtained for determining the number of cycles of 
stress at any arbitrary reference stress level and 
range to cause the same damage as the number 
of cycles of stress at any stress and its appropriate 
range (and, in acoustic fatigue, its appropriate 
frequency). This relationship—or, what amounts 
to the same, the coordinate transformations from 
the o,, — o,'), An, plane to a,, (o, — a,’), 
(An,), plane—are provided by Eq. (47). This 
equation enables one to treat the problem of fatigue 
for any arbitrary intermixing of stresses, stress 
ranges (and, by implication, temperatures and 
frequencies) as if the whole damage is occurring at 
one stress amplitude, temperature, and frequency, 
and discuss the problem at the reference state 
i.e., the reference coordinate system. 

It was mentioned earlier that any theory is only 
as good as the assumptions implied in its develop- 
ment. It is therefore highly desirable to call 
attention to the limitations implied by the assump- 
tions in the present theory. The first is the fact 
that strain softening, which appears to be an im- 
portant aspect of fatigue under certain conditions, 
is not physically taken into account and, hence, 
cannot be explained by this theory. The second is 
that there is some evidence to indicate that, if a 
stress lower than the endurance limit is included in 
a program, it appears to make some difference in 
the test results. This would certainly mean that 
the form of Eq. (17) must respond in some manner 
to this effect. Presumably, it may be explained 
on the basis that the value of o; changes in some 
manner in such a test, but there is as yet no way 
of taking it into account. The third is that there 
is also some evidence to indicate that the order of 
application of the loads in program-type loading 
alters the fatigue life in a more significant manner 
than the present theory would indicate. The 
reported test results do not appear to be statistically 
reliable, and this aspect must be investigated 
further. 

The analysis of cumulative damage in random 
loading presented here is a first-order approximation 
only. More detailed investigations are under way 
at the present. A formal quantitative theory of 
temperature and frequency which appears to explain 
many of the known aspects is substantially com- 
pleted. 


Conclusions 


An analytical theory of fatigue for high stress 
levels has been formulated based on a fairly simple 
model derived from a synthesis of the macroscopic 
elastoplastic fracture theory and concepts derived 
from the dislocation theory of metals. The theory 
(1) predicts the shape of the stress vs. number of 
cycles curve in fatigue; (2) after determining two 
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of the unknown constants from test results, makes 


quantitative predictions regarding the residual 


static strength, changes resulting’ from change of 
test stress ratio (and mean stress), cumulative 
damage due either to multiple step or to random 
loading; and (3) as a formal extension, gives a 
Once the 
constants have been determined, the theory asserts 
that no additional assumptions of a fundamental 
nature are necessary to predict or explain the 
sensitivity of cumulative damage to the order of 
The basic 
statistical nature of the fatigue problem has been 
attributed to two material variables of which one 
is grain size and the other, probably more important, 


quantitative theory of acoustic fatigue. 


stressing and the stress magnitude. 


a mean internal stress denoted as oa; which is essen- 
tially a stress that opposes the motion of dislocations. 
This stress, characterized at present by its mean 
value, depends on the prior stress history and is 
expected to vary throughout the material with its 
variation in terms of mean value and its distribution 
not known. It is suggested that, because of this, 
the fatigue problem will remain statistical in nature. 
In the present theory this stress o; has been identi- 
fied with the endurance limit under completely 
reversed loading. Confirmation of the theory 
through tests performed by different people has been 
satisfactory, though considerably more experimental 
verification is needed. When statistically reliable 
basic test data are available, quantitative pre- 
dictions from the theory for program-type loading 
have been found to be excellent, the average error 
of prediction being about 5 percent, a value well 
within the requirements of most designers and 
engineers. The general form of the equations of 
this theory of fatigue is such that the treatment of 
the various aspects of the fatigue problem from a 
unified viewpoint appears to be feasible. 
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Man—Machine Dynamics of Nose-Wheel Steering (continued from page 21) 


Thus the force output required of the operator is 
seen to be a complex function of a number of quanti- 
ties, including several variables. 

Up to this point we have made no assumptions 
concerning the primacy of one or another perceptual 
input, so that the results are perfectly general. If 
we now assume that the variable y is the pilot’s 
primary visual cue in steering and F his primary out- 
put, then the function F/y represents a simple linear 
approximation to his transfer function. 
(13), we obtain 


From Eq. 


F/y = — D)/D) (14) 


This statement represents the assumption that 
the perceptual variable x (lateral distance from 
centerline) does not play an important part in the 
steering task. It does not, however, tell us any- 
thing about the variable D other than that it must 
be taken into account by the operator. 

Some observations on steering behavior are in 
order at this point. Most drivers are probably 
aware of the fact that they tend to look farther ahead 
at high speeds than at low speeds. This is a logical 
outgrowth of the fact that, given a particular control 
task, a human operator is characterized by a more 
or less constant response time so that his anticipation 
requirements dictate that he “‘look ahead”’ at least a 
distance equal to response time divided by speed 
(this is then the distance traveled during the interval 
of one response time). It may be possible to quan- 
tify the relationship experimentally by means of 
precise eye movement studies to determine fixation 
point as a function of vehicle speed. 


RUNWAY 
CENTERLINE DEFINITIONS 
a x= DEVIATION NORMAL TO 


RUNWAY CENTERLINE, FT. 


= DIRECTION OF MOTION OF 
AIRCRAFT WITH RESPECT 
TO RUNWAY CENTERLINE, 
RADIANS 


D= INSTANTANEOUS DISTANCE 
TO RUNWAY CENTERLINE 
ALONG DIRECTION OF 
MOTION, FEET 


NOTE: x= yD FORW= 10° 


(ALL ANGLES AND DISTANCES 
ARE EXAGGERATED FOR 
CLARITY ) 


Fig. 2. Geometry of visual perception problem associated with 
steering an aircraft. 
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In line with these observations and in the absence 
of definitive data, then, let us make two further ex- 
tremely important assumptions: (1) that the pilot 
looks a ‘‘distance ahead,”’ d, which is directly propor- 
tional to speed—viz., d = rv, where 7 is his response 

and (2) that he continually at- 
tempts to maintain the condition D = d by adjust- 
ing w via the medium of the nose-wheel control. The 
latter assumption can best be illustrated by assum- 
ing that the pilot extends an imaginary line of 


time in seconds; 


length d = rv ahead of the aircraft and pointed in 
the direction of motion, and attempts to keep the 
far end of this imaginary line exactly on the center- 
line of the runway by his steering behavior. We 
may then introduce the resulting relationship repre- 
senting these two assumptions: 


D = ww (15) 


Now D = 0 since in the present analysis v is also a 
constant (from a previous assumption). With 
these relationships in mind we can rewrite Eq. (14) 
as follows: 


F/y = W/kRD = W/rkv (16) 


|The significance of Eq. (16) will be discussed later. | 

Next, let us suppose that x—1.e., lateral deviation 

from centerline—is the pilot’s primary visual input 
instead of y. From Eq. (16) we see that 

F = (p/k)(W/D) (17) 

But y = x/D, from Eq. (7). 
Eq. (17) we obtain 


Substituting this in 


F = (x/k)(W/D?) (18) 


As in Eq. (15), we assume again that D = 7. 
Substituting this in Eq. (18) and solving for F/x, we 
obtain the desired result: 


F/x = W/krv? (19) 


It is instructive to assess the relative meanings of 
Eqs. (16) and (19) as possible human transfer func- 
tions. First, we see that both are dependent only 
upon the three constants W (wheel base), v (speed), 
and 7 (response time). Second, we can determine 
whether stable steering system operation is possible 
in either case by solving both Eqs. (16) and (19) for 
F, then substituting successively in Eqs. (12) and 
(6 


From Eq. (16), F = Wy/krv. Substituting this 


in Eqs. (12) and (6) we obtain ¢ = vy/r. But y = 
x/D, from Eq. (7). Therefore, 
* = (20) 


Again applying the assumption of Eq. (15) and 
substituting for D in Eq. (20), we obtain 
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Fig. 3. Human transfer gain required to compensate for different 
constant vehicle speeds under differing assumptions of perceptual input. 


The well-known solution to Eq. (21) is 
me (22) 

This means that if the man exhibits the transfer 
function defined in Eq. (14) and uses y as his pri- 
mary perceptual input, the aircraft will always ap- 
proach the runway centerline in an exponential man- 
ner for any given displacement x) (the constant 1/7 
represents the decay rate of the resulting exponen- 
tial). 

Now let us look at the second possibility. Eq. 
(19) represents the case where the man uses x as his 
primary perceptual input. Substituting Eqs. (19) 
and (12) in Eq. (6), we obtain directly Eq. (21) 
again. So it is evident that both methods of per- 
ceiving apparently lead to stable operation in the 
nose-wheel steering task. How then do they differ? 
To illustrate the difference graphically, both gain 
functions are plotted in Fig. 2 for r = 0.5 sec and 
W = 22 ft. 

It is clearly evident that the human’s transfer 
gain must vary over a wide range to handle different 
vehicle speeds—in the one case, his gain must be 
inversely proportional to speed, while in the other 
case it must be inversely proportional to the square 
of speed. Pilot’s transfer gain in either case is seen 
to be independent of time and/or frequency and may 
be characterized quite simply as a velocity-modulated 
gain. 

It is quite possible that, in general, the operator 
may include components of both transfer functions. 
For instance, his output might be characterized as 
follows: 


(23) 
where c, and ¢, represent the relative proportions of 
each function being utilized by the operator at any 
given time. Indeed, it is possible that these relative 
proportions may vary with learning, with time, or 
even with ground speed. 

However, it appears more reasonable to assume 
that, since transfer gain variation with the square of 
speed is the more difficult task, requiring greatest 
sensitivity and largest variation, the operator will 
eventually learn or will prefer to steer the vehicle 
using only aircraft heading y and distance to center- 
line D. The considerable weight of opinion in favor 
of this assumption is summarized by Gibson! thus: 
“The behavior involved in steering an automobile, 
for instance, has usually been misunderstood. It is 
less a matter of aligning the car with the road than it 
is a matter of keeping the focus of expansion* in the 
direction one must go.”’ (Italics mine.) 

That a changing gain as a function of speed is re- 
quired of the pilot is not a new observation; this 
requirement has been recognized, at least implicitly, 
in connection with the steering of all types of ground 
vehicles (see the previous discussion of ‘‘distance 
ahead’’). Nonetheless, the extent of the variation 
in required gain as indicated in Fig. 3 is rather 
startling. It derives from the fact that nose-wheel 
steering efectiveness must be maintained throughout 
a range of speeds up to that at which the increasing 
aerodynamic effectiveness of the rudder can provide 
adequate lateral control of the aircraft. It is likely 
to be at this speed that nose-wheel lift-oif is per- 
formed, relieving the pilot or copilot of the extreme 
nose-wheel steering sensitivity required because of 
high speed and permitting the lower sensitivity of 
rudder steering to become effective. In eifect, nose- 
wheel lift-oif represents a sudden change in vehicle 
steering gain of a sort indicated in Fig. 4. This 
change serves to simplify the pilot’s steering task at a 
crucial point in the ground roll. 

In contrast, at the low end of the speed range, 
very high nose-wheel gain is required for turn-oif, 
taxiing, parking, and initial alignment on runway 
during start of take-off ground roll. 

Taken together, these factors indicate a rather 
complex steering task. As take-oif speeds have 
increased with the advent of jet aircraft, the inade- 
quacies of linear steering systems with respect to 
the pilot’s steering task have become increasingly 

* The focus of expansion ia the visual field always repre 
sents the direction of motion of the observer. In steering on 
a runway, relative motion of the aircraft with respect to the 
runway centerline is simply the angle y formed by the inter- 
section of the runway centerline with an extended line from 
the pilot to his focus of expansion, as shown in Fig. 2. In the 
absence of wind, skid, or creep, this line is usually coincident 
with the centerline of the aircraft as shown in Fig. 1. Never- 
theless, the foregoing analysis is perfectly valid in the pres- 
ence of a crosswind. 
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serious and, hence, apparent. If, for example, a 
nose-wheel steering gain is selected for best operation 
at 60 mph, the gain will be found much too low for 
low-speed taxiing and parking. In an attempt to 
overcome this particular difficulty, it has been pro- 
posed that a two-valued nonlinear steering gain be 
incorporated into steering systems as exemplified in 
Fig. 5. The higher gain would be obtained with the 
large excursions of the rudder bar normally encoun- 
tered in taxiing and parking. 


Proposed Solution 


Of the various solutions proposed, none has gone 
far enough toward a complete solution valid for all 
conditions of the ground-roll phase of aircraft oper- 
ation. It is proposed, then, that such a solution be 
attempted, based upon known facts and principles 
of human performance in a servoloop. 

It appears desirable from a human factors stand- 
point for an operator in a servoloop to maintain a 
high, fixed gain and to act as a simple amplifier. 
Experimental evidence indicates that such opera- 
tions will result in least variability and greatest ac- 
curacy in performance.?~> 

Translating this statement into terms of the 
present problem, it means that a given rudder pedal 
(or steering wheel) force should always result in the 
same rate of change of aircraft direction of motion 
y, regardless of speed [refer back to Eq. (4)]. In 
order that this may occur, the steering system must 
be so designed that the nose-wheel angle 6 is made 
directly proportional to the product of applied rudder 
pedal force F and wheel base W and inversely pro- 
portional to ground speed v—viz., 


6 = kF(W/v) (24) 


where & is the constant of proportionality as in Eq. 
(12). Eq. (24) specifies the operation of a velocity- 
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Fig. 4. Variation in required pilot's transfer gain as a function of 
speed, showing effect df nose-wheel lift-off. 
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Fig. 5. Steering gains in current use. 


modulated steering system. To see whether this 
mechanization has the desired effect, first substitute 


Eq. (24) in Eq. (11). This results in 


F = (y/k) [(v — D)/D] (25) 


The contrast between Eq. (25) and Eq. (13) is 
evident in that the dependency upon wheel base and 
inverse velocity has been eliminated. 

Next substitute Eq. (15) in Eq. (25) to obtain 


F = — (70/v)] (26) 


For the assumption of constant speed, Eq. (26) 
reduces to 


F = y/rk (27) 


This equation says that, with the foregoing 
mechanization of the steering system, the force out- 
put required is simply proportional to the perceived 
angle y regardless of speed. 

Since constant speed has been assumed in the de- 
velopment to this point, it is instructive to examine 
the new requirements which the proposed system of 
velocity-modulated steering places upon the human 
operator or pilot under the more usual conditions of 
varying speed. More specifically, the case of par- 
ticular interest is that of take-off, wherein relatively 
constant acceleration occurs so that d = a constant. 
For aircraft of interest, this value is of the order of 
8.8 ft/sec’, representing a take-off run of 3,900 ft 
over 30 sec. The value of 7, or human response 
time, can be taken as approximately 1/2 sec. Sub- 
stituting these values in Eq. (26), we obtain 


F = (p/0.5k)[1 — (4.4/v)] (28) 


The quantity rkF/y is plotted against [1 — 
(4.4/v)| in Fig. 6 to indicate the required variation 
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in pilot’s gain with constant take-off acceleration 
using the proposed velocity-modulated steering 
system On the same graph is plotted the solution 
for the constant steering gain system. The success 
of the proposed system in reducing pilot transfer 
gain variation over a wide range of ground speed is 
evident. It does not require increasing pilot steer- 
ing sensitivity with speed as does the normal con- 
stant gain system. Over the important speed range 
of 30-80 mph the pilot gain requirement is essenti- 
ally constant. Its specific level is fixed by the ap- 
propriate selection of k. 

The type of mechanization required for such a 
steering system is indicated schematically in Fig. 7. 
The portion of the diagram within the dotted lines 
already exists on most aircraft of interest, the only 
difference being that a constant voltage is applied 
across the Request Potentiometer at A”-B”. In- 
troduction of velocity-modulated steering in these 
aircraft requires only that this applied voltage be 
made an inverse linear function of ground speed, 
instead of a constant. The scheme shown in Fig. 7 
is only one of a number of possible schemes, the 
particular choice of which depends upon the type of 
steering system in an aircraft, the selection of a re- 
liable speed sensor, and the desirability of introducing 
fail-safe features. References 6 and 7 describe other 
mechanization techniques which will achieve the 
same goal after suitable modification. (It is of in- 
terest to note that those aircraft presently equipped 
with anti-skid braking devices have inherent in 
them a ground-speed sensor.) 

Referring back to Fig. 7, it is seen that the steering 
signal voltage is directly proportional to pilot’s 
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Fig. 6. Comparison of human transfer gain requirements for 
standard vs. proposed steering systems under conditions of constant 
longitudinal acceleration. 
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Fig. 7. Schematic diagram of velocity-modulated steering system. 


force F (which produces rudder bar or steering wheel 
deflection), aircraft wheel base W, and a constant 
of proportionality k; and inversely proportional to 
speed v; all as required by Eq. (24). Such a 
scheme would indeed result in high gain at low speeds 
for turn-off, taxiing and parking, and initial ground 
roll during take-off, and would result in low gain at 
high speeds where overcontrol might be critical. 

The most attractive result, however, for those 
aircraft with rudder pedal control of nose-wheel 
steering is the fact that rudder pedal travel can be 
large at high speeds, thus permitting full use of the 
rudder as soon as it becomes aerodynamically effec- 
tive with no sharp rudder transition. 

An additional important advantage is that, 
through proper choice of the gain constant k, the 
normal side-load on the aircraft due to nose-wheel 
steering can be operationally limited to some nom- 
inal amount (as required by some current aircraft 
specifications) while at the same time allowing the 
pilot to override this limit in certain situations where 
a ground loop might be the lesser of two evils. 

Finally, the proposed velocity-modulated steering 
system concept retains the safety feature of requiring 
that the pilot look ahead a distance proportional to 
speed, without imposing on him the additional 
burden of varying his steering transfer gain with 
speed as required by present systems. 
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Flight Mechanics Reliability 
existent. A much higher level of effort is needed in 
the collection, correlation, and dissemination of 
information on components and subsystems. The 
initiative here can be taken by the Air Force, 
particularly in improving the organization and 
unification of the program to place emphasis on all 
aspects. However, the willing cooperation of the 
industry is most essential. 

The solution of the problem of design of reliable 
aerospace vehicles demands careful and complete 
analyses at every level of industry, from the smallest 
item up to components, subsystems, and systems. 
After these careful analyses, judgments can be made 
as to re-examination of the requirement, develop- 
ment of ‘“‘super’’ components, redundancy, testing, 
and inspection rigidity. Fundamentally this job 
can be done only by industry. The Air Force can 
exhort, establish requirements and goals, participate 
in compromises, etc., but it is not close enough, 
early enough, to the design to be able to do the job. 
Neither can the job be done well by the specialized 
reliability groups in industry. These groups have 
a utility, but the key to the situation is the indoctri- 
nation of all levels associated with the design and an 
almost fanaticism on the part of group leaders to 
accept nothing less in reliability than what is 
needed regardless of what it means in discarding 
old ways of doing things. Only when this indoc- 
trination is accomplished will we be on the way to 
solution of the problem we face. 
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Relating Industrial Research to Company Goals 


including whether emphasis will be placed on basic 
or applied research; over what time spans are the 
results of research expected to have an economic 
impact in the company; and how follow-on develop- 
ment is to be supported. Successful research leads 
to ideas, and ideas lead to product development. 
The extent to which the research team participates 
in the development will determine to a large degree 
the stability of the research organization. If re- 
search scientists are encouraged to continue with 
development, research will be stopped. On the 
other hand, if development derives little benefit 
from the research background, the development 
project may fail. Policy must be established upon 
which continuity from research to development is 
maintained with maximum advantage to company 
goals 

Under the topic of facilities, decisions are needed 
with regard to questions such as whether emphasis 
will be given to facilities for research or for develop- 
ment; geographic accessibility of the laboratory 
from the home company, employees’ homes, local 
universities; who will operate the facilities; and does 
the long-range financing of the facilities provide for 
all important needs? 

Organization of the research function is important 
because it determines to a large extent the environ- 
ment in which scientists and engineers doing re- 
search must work.* Organization determines not 
only the immediate working environment but also 
the extent to which the research director may control 
facilities operation, support services, procurement, 
etc. Decisions on these details are required. 

Finance involves topics such as time span of re- 
search budgets (one year or more); pricing policy 
of research contracts; circumstances under which 
the company sets aside the requirement for profits in 
research contracts; 


impact of Armed Services 


procurement regulations on research financing; 
and accounting practices in determining research 
costs 

Relating the research mission to company goals 
requires careful consideration of all details just 
enumerated. Reconciliation of the issues involved 
in the company goals with the factors involved in 
the conduct of research requires top management 
decisions.‘ Although the problems are complex, 
skilled managers are solving equally difficult prob- 
lems in other areas.° 

To SUMMARIZE—a problem exists in matching re- 
search with business. The problem is important 
and complex, but there is no reason to believe that 


it cannot be adequately solved. A recognition of 
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the principal details involved in company goals and 
research, an appreciation for new management 
techniques involving complicated details, and a 
willingness for all—including top management— 
to participate in the approach to the problem are the 
essential ingredients in its solution. 
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(b) Transistors rejected via temperature cycling 
tests. 

(c) ‘“Gettered”’ transistors. 

Tests results were as follows: 

(a) Transistors accepted via temperature cycling 
tests showed approximately 4 percent temperature- 
sensitive units. 

(b) Transistors rejected via temperature cycling 
tests showed approximately 80 percent tempera- 
ture-sensitive units. 

(c) ““Gettered”’ transistors did not show tempera- 
ture sensitivity. 


This test answered several questions: 
(a) The temperature cycle screening was removing 
the majority of temperature sensitive transistors and 
removing all that would cause malfunctions at time 
of screening. 

(b) ‘‘Gettered’’ transistors did not show tempera- 
ture sensitivity at the time of the test. 

It asked at least one more: Would time cause 
temperature sensitivity degradation of temperature 
cycle screened transistors? 

The test groups were subjected to 96 hours at 
100°C (equivalent to one year’s storage at 25°C). 
The temperature plateau test was repeated, the 
groups were subjected to a second 96-hour period 
at 100°C, and the temperature plateau test was 
again performed. 

It was found that temperature sensitivity degrada- 
tion increased in transistors that were temperature 
cycle accepted. An approximately linear increase 
in the number of transistors exhibiting temperature 
sensitivity with time was noted. The severity 
of sensitivities also increased with time. A definite 
upward shift was noted in temperature sensitivity 
point with respect to temperature. ‘‘Gettered”’ 
transistors showed no evidence of degradation or 
temperature sensitivity. 

A program of remedial action was instituted on the 
basis of test results. Ona short-term basis, a tem- 
perature plateau test with plateaus at —35, —25, 
—15, +5, +15, and +25°C was established to 
screen existing stocks. Steps were taken to reduce 
transistor high-temperature experience in storage, 
manufacturing, and _ test. 
action is a change to transistors of ‘“‘gettered’’ con- 


Long-term corrective 


struction. 

The deficiency noted would not affect equipment 
operating at average or high ambient conditions 
except as an undesirable shortening of transistor life. 
However, this phenomenon should be carefully con- 
sidered in the design of equipments where low-tem- 
perature operation is desired. 
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advertising appears in this issue 
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you know what they are doing 
to help you solve research and 
design problems in_ missile, 
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projects. 
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plies useful to your professional 
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Engineering is now underway on the strange look- 
6,7 | ing bird pictured here. We call it the “Humming- 
bird,” and we are developing it for the U. S. Army. 


It is to utilize the jet ejector augmentation principle 
39 for VTOL aerial vehicles to fly straight up, straight 
down, hover in mid-air, and convert to conventional 
over | forward flying at speeds in excess of 500 miles per 
hour. Just one example of what is going on here 
2 | at our Engineering Center where dreams come true. 
Like to join us? 


AIRCRAFT ENGINEERS 


Openings in the fields of: Aircraft Design Engineering * En- 
gineering Drawings Checking * Aircraft Structures * Flutter 
and Vibration * Sonic Fatigue * Aircraft Specifications En- 
gineering ¢ Reliability Engineering * Operations Research 


¢ Aircraft Research Engineering. 


Write to: Hugh L. Gordon, Professional Employment Megr., 


Lockheed-Georgia Co., 834 W. Peachtree St., Atlanta 8, Ga. 
Dept. Z-77 


All qualified applicants will receive consideration for em- 
ployment without regard to race, creed, color or national 


origin. 


|THE ENGINEERING CENTER 


“\LOCKHEED -GEORGIA COMPANY 


A DIVISION OF LOCKHEED AIRCRAFT CORPORATION 
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If you are an 

advanced degree physicist, 

scientist or electronic 

management engineer who seeks to 
shape air and space vehicles to 
their unique environmental 
problems rather than stretch existing 
designs for marginal success, 

you are needed at McDonnell. 


Write R. F. Kaletta, Professional Placement, 
McDonnell Aircraft, St. Louis 66, Missouri 


AN EQUAL OPPORTUNITY EMPLOYER 


lified Sagicants will receive consideration for 


yment without regard to race. creed, color or national orig 


\ 
<j 
The | 
% 
4 


‘ 
© 
i 
2 
iy 
: 
i 
Ag 


